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APPENDIX A 

Parametric Modeling Of 
Discounted Benefit-Sacrifice 
Streams 


1 


Bene fit- cost flow f( + ) 


INTRODUCTION 


This Appendix provides an introduction to modeling of benefit- 
sacrifice streams by parametric approach: benefit and sacrifice 

streams are shaped in terms of mathematical functions with an appropri- 
ate number of shaping parameters that are interpreted in seme relevant 
manner. Each parameter can be treated also as a "random variable" 
whereby the benefit-cost streams can be subjected into effects of a 
multidimensional uncertainty. Although this Appendix treats primarily 
internal rate of return models , the extensions can be made to any 
other discounted streams . 


SOME BASIC DEFINITIONS AND CONCEPTS OF BENEFIT -SACRIFICE STREAMS 





Figure A-l 


Figure A-l illustrates a typical net benefit-cost or benefit-sacrifice 
stream given by a function f(t). By definition, when the net f(t) is 
negative, the flow is called net sacr ifice flow or "investment flow", 
-i{t); and when it is nonnegative , it is called net returns benefit 
flow, +n(t). i(+) is then a positive flow running out from the pocket , 

while n(+) is a positive flow running into the pocket. 

Further, 00 

K = / n(+)dt 
t=0 


2 


is the total volume of sacrifice or total invested value in the program 
while 

03 

N ~ j n(+)dt 
t=0 

is the total net returns benefit volume. 


The internal rate of return from the program, r, is of the form 


r - f 1 

Tpb \ 


Where Tpb is the so-called pay-back period possible only if N>K in 
Figure A-l. f f a monotonic increasing function of the benefit to 

sacrifice volume ratio N/K. The payback period Tpb for N>K is defined 
as follows: 

Tpb 

K = / n(+)dt 


In Figure A-2, also the break-even point Tv e is indicated. 



The discounted net benefit-net sacrifice flow g(t) is related to 
the net benefit-net sacrifice flow f(t) as follows: 

gCt) = fit) e“ rt 

r here is the rate of discounting and has dimension 1/year or inverse 
year. Figure A-3 illustrates a flow f(t) and its discounted flow gCt). 

The internal rate of return on investment, r, of the total flow 
f(t) is that rate of discounting for which 


/ f(t) e rt dt = / g(t) dt = 0 
0 0 

i.e. is that value of r for which the discounted total value of volume 
invested equals to the discounted total net returns volume . Let K' 
be this discounted total dollar value invested, and N' be the total 
value volume net returns. Then the internal rate of return on invest- 
ment is that value of discount rate r for which 

J** - YVh - 00 — r*t 

K* = / i(t)e dt = / n(t) e dt = N' 

0 0 

One notes then that there is this one particular rate of discounting 
that satisfies the above equality on the nose. 



Figure A-3. 
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For a POSITIVE rate of return on investment, r, one must have initially 
N>K in order that N 1 = K'. If N = K, then r =0, and if N<K, then r must 
be made negative in order to make N T = K* . The effect of how n(t) and 
i(t) are distributed over time will become clarified by numerous sub- 
sequent examples. 

The method of continuous discounting has been chosen here rather 
than the conventionally accepted method of discrete discounting. There 
are some very good reasons for this. First, continuous discounting turns 
out to be mathematically simpler than the discrete discounting method. 
Second, it is much easier to find the right and unique value of the rate 
of return on Investment for a particular problem In the case of continu- 
ous discounting than in the case of discrete discounting. Third, it is 
much easier to introduce parametric cash flow models subject to continu- 
ous discounting than those subject to discrete discounting in order to 
find the rates of returns in terms of theee parameters in closed forms. 
Fourth, it is rather simple to discount discrete cash flow models by 
continuous discounting so that basically nothing is lost even in dis- 
crete cases by using continuous discounting techniques. Fifth, appro- 
priate dimensionality of various financial and other business parameters 
becomes properly accounted for by the use of continuous discounting 
techniques. The general overall advantages of continuous cash flow 
and discounting techniques over the conventional discrete methods of 
discounting becomes obvious as one introduces several examples in the 
subsequent treatment. It is useful to review a simple case of discrete 
discounting and then show how one can introduce from it the case of 
continuous discounting. 


A REVIEW OF CONCEPTUAL RELATIONSHIPS AND DIFFERENCES BETWEEN THE DIS- 
CRETE PERIODIC DISCOUNTING METHODS AND THE CONTINUOUS DISCOUNTING METHODS 


A commonly used conventional method of discounting is the discrete 
one where the event of discounting occurs at definite intervals of time. 
Typically the event of discounting occurs once per annum, semi-annually, 
or quarterly. 

Consider a program for which it is desired to invest K dollars in 
a productive asset at the beginning of the first period (k = 1). Let 
be the lump sum dollar net return for the k 1 " period and let J be the 
residue value of the asset in dollars at the end of its life at the with 
period’s end. Then the usual "present value equation" for this cas^ is 
as follows: 


(1 + r Q ) k + (1 + r 0 ) m 
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It is very important to note that Nj_» N 2 , N^» N m are all sums 

of net dollars obtained over the respective periods. That is, the dimen- 
sionality of is dollars. Thus the above sequence of dollar sums is 
not yet properly a "cash flow" although it is so called at times. Further, 
one should note that r„ here is a dimensionless discounting factor (i.e. 
not a discounting RATE) that can be eventually related to a discount- 


ing RATE. If K,J,N ,N, 


are all given then the problem is 


ii XJ » -L. X jU jin 9 • • • 9 4 * • 3 Ifl Lilt J-Cj 

to find an appropriately UNIQUE value r Q satisfying the above equation. 

By a fundamental theorem of algebra one will find m roots for the poly- 
nomial in r Q from the above equation. Many of these roots could be 
negative or complex, and many of them could be positive. A tedious 
problem facing an analyst is to pick up a unique real valued root that 
corresponds to the RELEVANT fraction of return on investment r n for 
his problem. It will be seen that this tedious aspect of the discounting 
becomes greatly simplified by shifting to the continuous discounting 
methods . 


The relationship between the fraction or return on investment r n 
and the rate of return on investment r {note, r must have a dimension 
of INVERSE TIME) can be handled by introducing a standard unit period 
of time (say, one year) 0 whereby 

v = r Q /0 or r Q = r 0 . 

Thus, since r Q is a dimensionless fraction, r now has the dimensionality 
of inverse years or 1/years. The percentage rate of return on investment 
is just 100 r. In terms of the standard period 0 and the RATE of return 
on investment r one can rewrite the present value equation in the follow- 
ing form: 


N /(I + 0 r) 


+ J/(l + 6 r) 1 


Here then the discounting is done once per period 0 (say, once a year). 

Assume now that one wants to discount p times per period 0 where 
p is a positive real number greater than one. Thus the NEW discounting 
period is NO MORE 0 but is now 0/p. Further, the old index k = 1,2,... 
..,m-l, m must be replaced by a new index k(p) which depends on p, k(p) 
1,2, 3, 4, pm. Note also that the net return dollar sums over the 

old periods of length 0, N^,N 2 ,N m are now to be replaced by 
new dollar sums of about only 1/p times the old values over the re- 
spective time periods. However, one still has the following equality: 


k(p)=l 
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Note , the old k^ interval of time was divided to p parts so that over 
it the old would be replaced roughly b new values N^/p = N, , Thus 
one now has a present value equation where the discounting is aoSe not 
once per period 6 but p times over this period, or once per a new period 

0/p: 


pm 

K = N k(p) / [1 = (9/p)r] kCp) + J/Cl + (0/p)r 

k(p)=l 


pm 



i 

k(p)=l 


CCN k( /(0/p)](0/p)]/[(l + 


ft p/ 9r 
ir) ] 


(0r/p)k(p) p/0r (0r/p)mp 

+ J/CU tip) 1 


It is now relatively easy to see what happens if one lets p become 
large or 0/p, respectively, small. As p approaches 00 , one can then in- 
troduce the following notations and limit values : 

1. 0/p -*■ dt, an infinitesimal increment of time; 

2. [6/p] k(p) + t, the continuous time variable; 

3. N k( n)/(®/p) ^NCt)/dt = n(t) the true CASH FLOW with the correct 

’ r ’ dimensionality of DOLLARS/YEAR; n(t) is then the 

expected net returns cash flow. 

4. Further one obtains the following well known limit: 

lira [1 + (0/p)r]P /6r = e 

p 


Thus, as p goes to infinity, one obtains in the limit' the following 
continuously discounted present value equation: 


K = / n(t) e” r1: dt + J e 

t=0 


If K, n(t), J and T are given then one Is supposed to find the rate of 
discounting r that satisfies the above equation. This r is then called 
the rate of return on Investment, For this equation recall now 

K is the sacrified volume generating productive asset. 

T is the planning time horizon or the planned life of the asset; 
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J is the expected residue value of the asset at the end t=T 
of its life* 

n(t) is the expected net returns benefit flow in dollars per 
annum from the venture utilizing the productive asset; 

r is the rate of return expected from the venture , in 1/years ; 

t is the continuous time variable in units of years. 

One can also define the pay-back period T , for this case as follows 
7. is that period of time for which 


^"Tpb 

K - J = / n(t)dt. 

t=0 


It should be noted that if one cannot find a finite Tp^ for a given K-J 
and n(t ) , then one should not get involved in such a venture. It should 
be intuitively clear that a venture should have a finite pay-back period 
in order that it would have a positive rate of return on investment. 

Insofar as future events are concerned, one does not deal with 
certainties but rather with expectations subject to uncertainties. Thus 
T,K,J,n(t), T and r as well as T ^ are expected values in most cases 
when contemplated for future actions. It should also be clear that all 
such expected values are subject to the relevancy of technological, 
financial, production, etc. information necessary for putting together 
a meaningful benefit-sacrifice flow picture for the contemplated venture. 
Typically such pieces of information must be integrated together to a 
total meaningful picture or pattern caricaturing a venture in a reason- 
able manner. Once a relevant model for the venture is obtained, one can 
proceed to find whether a finite pay-back period would exist. If it 
does , then it is appropriate to find the magnitude for the rate of re- 
turn on sacrifice for this contemplated venture. 


REVIEW OF MATHEMATICAL CONCEPTS USEFUL FOR BENEFIT-SACRIFICE FLOW MODELS 


The subsequent treatments of several benefit-sacrifice (B-S) flow 
models can be greatly faciliated and simplified by utilizing concepts 
of applied mathematics . For the cases of continuously discounted flows 
techniques of Laplace transforms and related mathematical concepts are 
useful. 


UNIT STEP FUNCTION u(t-T) is defined as follows: 

u(t-T) = 1 for all t>T 

= 0 for all other values of t. 
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An example of a linear combination of unit step functions is the follow- 
ing "pulse function": 

u(0 - u(t-T) = 1 for 0 t T 

= 0 for all other values of t. 

This kind of a pulse function allows one to cut off an appropriate portion 
of any time function. For example, 

T 03 

/ n(t) e -rt dt = / Cu(t) - u(t-T)] n(t) e “^dt 
0 0 


DELTA FUNCTION OR IMPULSE FUNCTION 6 (t-T) is defined as follows: 
5 (t-T) = 0 fol all t i T 

CO 

/ 6 (t-T) dt = 1 


An immediately useful application for this function is the case where K 
dollars are laid ou" instantly for the purchase of an asset. In this 
case the investment outlay cash flow would be as follows : 

- i(t) = - K6(t) . 

This would be then the investment outlay cash flow in the "present value" 
case, as will be pointed out later on. 


THE LAPLACE TRANSFORM OF A FUNCTION f(t) is defined as follows: 

+*» 

L[f(t)] = F(r) = / f(t) e -rt dt 
0 

The Laplace transforms for many different given time functions f(t) 
are tabulated in tables. 

Some useful examples of Laplace transforms are the following ones: 

L[i 0 u(t)] = i 0 /r L[K5(t)] = K 

— vi T wv 

L[n 0 u(t-T)] = (n c /r) e ' L[N<5(t-T)] = N e 

For any function f(t) that is shifted tc the right by the amount T of 
time one has : 

L[u(t-T)f(t-T)] = e~ rT L[f(t)] = e“ rT F(r) . 

Further discussion on how to use Laplace transform techniques will be 
done in the conjunction of several examples which are introduced sub- 
sequently. 
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PERIODIC FUNCTION: The function f(t) is said to be periodic if for 

n - 0 ,1,2 ,3,4 ,. . . it is true that f(t) = f(t-nT). Then T is called 
the period of this periodic function. 

If one has the integral of the form 

y _ rt 

/ f(t)e dt 

0 

then it can be rewritten into the following form: 

T 

/ f(t)e' rt dt = /”[u(t)-u(t-T)^f(t)e" rt dt = L[(u(t)-u(t-T))f(t)] 

0 0 

CO 

~ F(r) - / u(t-T)f Ct)e _rt dt . 

0 

If f(t) is either constant or periodic with period T* then 


os 

/ u(t-T)f(t)e -rt dt = e -r '^F(r) and so for this particular case 

0 

T t J f(t)e _rt dt 

/ f(t)e dt = F(r)[l - e -r "] or F(r) = - — 

0 [1 - e _rT ] 

Further, note if rT»l then one can use the approximation where 

T _ 

/ f(t) e“ rc dt ~ F(r). This can, indeed, simplify treatment of several 
0 cash flow discounting problems. 


THE INTERNAL RATE OF RETURN ON SACRIFICE 


The intrinsic or internal rate of return on sacrifice for a B-S 
flow associated with the respective venture is a commonly accepted and 
used measure for expected success. Consider the total B-C flow f(t) 
whose negative portions represent investment outlay (and thus with 
negative sign) and whose positive portions represent the expected net 
returns flow (hence with positive sign). If i(t) is the investment cash 
flow and n(t) is the net returns cash flow, then 

f(t) = n(t) - i(t) . 

In general then the condition for obtaining the internal rate of return 
on investment is as follows: 
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CO 


/ f(t) e rt dt = / CnCt) - i(t)] e'^dt 

0 0 

= / n(t) e -rt dt - / i{t) e“ rt dt = 0 


/ i(t) e“ rt dt = / n(t) e" rt dt 

0 0 

or L[i(t)] = l(r) = L[n(t)] = N(r). 


In this case it has been assumed that the residue value of the asset is 
zero, i.e. J = 0. This assumption is also reasonable if rT»l for the 
case at hand. If the residue value cannot be ignored, then 

f(t) = n(t) + J6(t-T) - i(t) 


and 


I(r) = N(r) + e r ^J . 

The "present value" case is one where i (t) = K6(t) whereby then 
l(r) = K. In such a case one obtains the condition shown below: 

K = N(r) + e -rT J . 

Here the outlay occurs at the time t = 0 in a one lump sum of a 
magnitude K so that it is a delta function with a strength K. 

In an interim summary, then, the basic equations of interest are 
as follows : 

In a general case 

-rT 

I(r) = N(r) + e J 

where I(r) = L[i(t)] and N(r) = L[n(t)], J is the residue value or 
benefit of the asset at the end of its life occurring at t=T where T 
is the life of the asset or planning time horizon, and r is the expect- 
ed internal rate of return on the investment. 


EXAMPLES OF DISTRIBUTED BENEFIT-LUMPSUM SACRIFICE MODELS 


The distributed benefit-lump sum sacrifice streams represent eases 
where the sacrifice lump sum is the present value of the discounted 
benefit stream, and the rate of discounting is the internal rate of 
return. 
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EXAMPLE NO. 1 


A program requires an asset of K dollars at the time t=0. Thus the 
investment cash outflow is just - i(t) = - K5(t). The net returns from 
the venture occurs in a lump sum of N dollars at the time t=T so that 
the net returns cash inflow is + n(t) = + NS(t-T). The overall cash flow 
of this venture is then f(t) = n(t) - i(t) = NS(t-T) - K6(t), For the 
internal rate of return on investment one has then the following require- 
ment: 


LCf(t)] = LCn(t)] - LCi(t)] 

= N e - K = 0. 

Solving for r yields the following expression for the rate of return on 
investment r: 

r = (1/T) In (N/K) . 

Figure A-4 illustrates the cash flow plot for this case. 


f(t) 

/N 


N<S(t-T) 

* 


t=TT 


K6(t) 

I 


W 


>t 


FIGURE A -4 


For this model K (capital outlay) and N (net returns volume) are 
financial parameters while T is a non-financial time parameter. It is 
noted that the rate of return on investment is inversely proportional 
to T and directly proportional to the logarithm of the financial ratio 
N/K. Note then how strongly r varies with T and how mildly it varies 
with N/K. This example illustrates simply how important the element 
of time is in ventures. One can readily obtain the total differential 
for r in order to see how sensitive r is to the changes in T and N/K: 


Ar = -E(ln(N/K))/T] 2 AT + Cl/T(N/K)3A(N/K) 


EXAMPLE NO. 2 


Let i(t) = K6(t) as before but assume now n(t) = n Q for 0<t“. 
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Then 


L[i(t)] = K = L[n(t)] = n 0 /i> . 


Thus r = n /K . 

o 

Note K is in units of dollars while is in dollars per year. Thus r 
has the correct dimension of a rate, i.e. inverse year. If n)t) = n Q 
for 0<t<T and zero elsewhere , then 

K = (n Q /r)[l - e“ rT ] . 

If rT>3 5 then exp(-rT) is much smaller than unity whereby r^n 0 /K is quite 
a good approximation. Note that in this case N - n Q T is now the net re- 
turns dollar volume so that one has the following expression: 

"rT, 

rT/(l - e ) = N/K 
or 

K = [(1 - e -rT )/rl] N 

The term in the square brackets is a dimensionless "time price of benefit 
value or money" 5 i.e. the worth of N when "purchased" by K. Note that 
for high rate of return r N is cheap and for low rate of return r it 
becomes dear. For high r little K is needed to purchase a unit of N; 
for low r much K is needed to purchase a unit of N. Each venture has its 
N or net returns value volume that must be "purchased" by investing a 
value volume K. Thus one has a clear relationship between the time price 
of money and the rate of return on investment. However , other time para- 
meters besides the rate of return on investment must enter into this 
expression of the time price of value. This will becomes more evident 
in subsequent discussions. It is for this and other good reasons that 
parametric flow models are extremely useful in venture analysis. Again 
one notes N/K enters in as a "financial" parameter or a bene fit /sacrifice 
ratio while T and r are time parameters. 


EXAMPLE NO, 3 


Consider a venture facing a competition with a constant force of 
mortality , h(t) = a, and an exponertial survival probability, exp(-at). 
Assume no renovation force. The expected net returns cash flow for the 
venture is assumed to be proportional to the probability of its survival. 
It is appropriate to note the mean life of the venture is iust 1/a in 
this case. Then the following net returns cash flow is assumed: 

n(t) = n Q e 

N = n c /a = n D T pl 

Here T pl = 1/a is the expected life of the venture in the competitive 
market. The capital outlay is K for this venture. Then one obtains the 
following result for the rate of return on investment: 

r = n 0 /K - a ~ [(N/K)-l]a . 
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Thus the internal rate of return on investment is reduced linearly by 
the force of product mortality. That is, the force of mortality to an 
added discount rate on the benefit stream. Note further: 

K = Ca/Cr+a)] N . 

Thus the time price of M is the higher a gets in comparison to r, and 
approaches unity in the ultimate limit. One notes again that a and r are 
the two time parameters associated with the expression for the time price 
of N. But a is now the frequency called the force of venture mortality. 


EXAMPLE NO. 4 


The previous model can be expanded to a bit more interesting one by 
introducing some additional relevant parameters. Figure A- 5 illustrates 
the basic setting. In this case the net returns flow need not rise 
instantly to its maximum value but rather gradually if this be relevant 
to the venture at hand. Thus the cash flow model has a t=T max point. 

It also has the force of product mortality, a. The other parameters 
incl ude a characteristic time T Q and the financial parameters K and N 
or their equivalents. Thus one has gained some additional degrees of 
freedom to portray more complex cases than the ones discussed previously. 
With these added degrees of freedom or descriptive' parameters one can 
shape net return flows with growth, maturation and decline character- 
istics found so often in everyday ventures. Yet one wishes to obtain 
closed form expressions that simplify greatly any computational efforts 
while retaining at the same time an easy pattern recognition for the 
total behavior of the venture. 
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The maximum expected net returns cash flow occurs at t=T max . 

The expected net returns B-S flow model that could portray or carica- 
ture such phenomena is as follows: 

n(t) = n 0 (t/T 0 ) k e" a t = [k! n 0 /T Q k l(t k /k ! ) e~ a * 

The constant n c , in value units per annum, scales the general amplitude 
of the whole flow cycle. The exponential term with the force of product 
mortality a will assure eventual decline of the whole net benefit flow. 
The characteristic time parameter T and the dimensionless parameter k 
can be used to caricature the growth portion while k and a can be used 
to portray the maximum cash flow time Assuming J=0 dollars for 

the residual value of the productive assets, one has the following re~ 
suit after appropriate discounting operation: 

K = [k! n 0 /T k ][l/(r + a)] k+1 

r = [k! n 0 /KT k ] 1/(k+l) - a 

One can easily verify that the maximum of the expected net returns 
benefit flow occurs when 


t = 


T 

max 


k/a , 


Further 9 the net returns volume is obtained as follows: 


co co -at 

N = / n(t) dt = (n /T k ) / t k e dt 

0 0 
k+1 

= k! V T o a 

One has then, also, the following maximum value for the net returns 
benefit flow: 

n max “ D</eaT 0 ] k n Q = (k/e)^CaN/k ! ] * 

This flow model can be then expressed in the following alternative 
but equivalent forms : 

n(t) = <a k+1 /k!) N t k e“ at - (ea/k) k n max t k e‘ at . 

The internal rate of return on investment could be expressed 
in the form involving the benefit-cost ration N/K, the force of mortal- 
ity a and the constant k: 

r = a C(N/K) 1/(1<+;L) - 11 . 

As an illustration, consider a case where the expected venture life is 
five years or, thus, the force of product mortality is a = 0.2 failures 
per year. Assume the net returns dollar volume N = 1,000,000 dollars. 
Let the capital outlay at the time t=Q be K- 200,000. dollars. Let then 


15 


too DATA 4*0*40*1*0* 7-5K4 

200 un A=cncf.)^{co.OTio/i > i^CAt i )+r:Ar'<c-o*;i) 4 ^> 
210 LET B*<lKA) + <CG-2t35/2) + CXT2> + rApC<-0. !!)-*,*) 
220 LET* C-C 1E6) 4)/6)*CKr 3)*^ < A) 

230 LET D=OE6>*C<0.e!T5)/£4) + <At4>*EXr<C-U*fc>*A) 
KUNKEV 


PL0TT0 20 £ S3 ACTS WE Do 09/03/69 


MULTIPLE PLOT OF THE FUNCTIONS A* B* C* D 


X-VALUE 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
n 
12 
13 
1 A 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 * 


H0KI2 INCREMENT - IB75 

0 7S000 

£oo&sS*»i»*oo3ai>5*JoiojJOi>Su3*o.Tae»t> 6*i c » g a ■ 


For Ap-r » 0.247 or 24.7 % 


A 

C B 

DC B A 

° D C c 0 For Ji 6 r=0.142 or H.2 % A 


B 


B 


b 

. . B 


A 

A 


D C 
D C 

D C 

DC BA 

D C BA 

D C A 

D CAB 

D CA B 

D A C B 

DA CB 

AD B 

A DBG 

A Db c For c r ^Q 0 099l;or 9.91 % 

A BD C 

A B DC 

A B C 

A B CD 

A B CD 

A B C D 

A B C D 

A H ~ C D 

A B c D For D, r =0.0759 or 7.59 % 

A B CD 


A B CD 

A B CD 

A b CD 
A B CD 
A B C D * 

A B C D 
ABC D 

AB C D 
AB C D 
ABC D 
ABC D 
AE C D 
Ab C D 

16 


FIGURE A-6 
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75000 


k vary from 1 to 4. First one wishes tc plot the net returns cash flow. 
The table below facilitates scaling a multiple plot of the four cases 
for comparison: 


k 

T 

max 

n max 


n(t ) 


1 

5 

7.35E+4 

A = 

(1E+6 )*( (0.2+2 )/l )*(T+1)*EXP(~ 

0.2*T) 

2 

10 

b . 40E+4 

B = 

(lE+6)*((0.2+3)/2)*(T+2)*EXP(- 

0.2*T) 

3 

15 

4.51E+4 

C = 

(lE-h6)*((0.2+4)/6)*(T+3)*EXP(- 

0.2*T) 

4 

20 

3.89E+4 

D = 

(1E+6 )*( (0 .2+5 )/24)*(T+4)*EXP( 

-0.2*T) 


One can now scale and provide appropriate plotting programming at some 
ease. The time scale should run from t=0 to t=40 years in steps of one 
year. For the cash flow the range should be from zero to 7.5E+4. 

The time price of money N in this case is [l+(r/a)]~^ + ^\ Figure 
A-6 illustrates the general shape of the cash flows as k goes from 1 to 4. 
The respective rates of returns on investment, r, are shown for each case. 


EXAMPLE NO. 5 


In the discrete periodic discounting process one has as many roots 
for r as there are discounting periods. This is no more so with the 
continuous discounting. Rather the number of roots obtained for r now 
depend on the general shapes of cash flows and is typically much smaller 
a number than in the discrete case. Further, the choice of an appropri- 
ate root for r is greatly simplified by considering at a relative ease a 
number of relevant limiting conditions that allow one to disregard all 
r except the relevant one. 

Consider a venture with seasonal variations. Management considers 
to purchase such a venture at the time t=0 for K dollars. The net re- 
turns flow is expected to have the following periodic form: 

n(t) = n Q [1 - cos w t] 

w = 2ir/T 

w is the angular frequency of the seasonal variation in the flow and T 
is the respective time period for the repetitive fluctuations in the 
net returns cash flow. Typically T is one year. Using the previous 
principles of discounting one obtains the following result: 

K = L[n(t)D = (n 0 w 2 )/[r(r 2 +w 2 )]. 


This is a cubic equation in r. If K, n Q , and w are, as they should be, 
positive real numbers , then this cubic equation has only one real root . 
The other two roots form a complex conjugate pair and are not used for 
any relevant interpretation for the rate of return on investment. The 
real root turns out to be as follows: 


1/3 C 1/3 

= (n Q w 2 /2K) ? Cl + A + (l/27)(2Kw/n 0 ) 2 : 




+ Cl - /l + ( 1/ 2 7 ) ( 2Kw/n n ) 2 ] 




One notes this egression behaves correctly when w*“ as in this case one 
has an average cash flow of a constant amplitude n Q for which in this 
limit r+n 0 /K. Also, as w-K> then n(t)->0 for whicu r-H3. 


In order to illustrate how the rate of return on investment, r, now 
depends on the seasonal periodicity w, consider a particular numerical 
example. Let K = 400,000 dollars for the purchase of an asset. Let 
n Q = 100,000 dollars per annum for the amplitude of the oscillation of 
the net returns flow. Then 


2/3 


r = (1/2) w ) Cl + /l + (64/27) w 2 ] 


1/3 


+ Cl - /l + (64/27) w ] 


1/3, 


Using time sharing computer terminal one readily finds the following 
results: 


w 

r 

0.00 

0.000000 

0.05 

0.075799 

0.10 

0.111475 

0.15 

0.136633 

0.20 

0.155677 

0.25 

0.1705 82 

0.30 

0.182482 

0.35 

0.1S2116 

0.40 

0.200000 

0.45 

0.206509 

0.50 

. 0.211927 


w r 


0,55 

0.216468 

0.60 

0.220301 

0.65 

0.223556 

0.70 

0.226337 

0.75 

0.22S727 

0.80 

0.230792 

0.85 

0.232586 

0.90 

0.234151 

0.95 

0.235524 

1.00 

0.236733 

1.05 

0.237803 


The ultimate limit for r is 0.25 as w increases toward infinity. Thus 
here the most sensitive range for r is when w goes from zero to unity 
years . 

Figure A-7 illustrates the behavior of the B-S flow for this kind 
of a periodic case. 
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FIGURE A- 7 
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EXAMPLE NO. 6 


If one assumes a case where the B-S flow has the form 

f(t) = n 'cos wt - K5(t) - 
o 

then by the previous discounting procedures one finds for r the following 
expressions: 


r = (n 0 /2K)[l - A - (2Kw/n Q ) 2 'J . 

which. r would be correct? If w-H) one should get a li" 1 "^ing case where 
n(t) = n Q for all t>0. Then + sign is appropriate sii. -e for w=0 one 
obtains the correct limiting case where r = n Q /K. Therefore, one can 
pidc up uniquely an appropriate rate of return on investment. The above 
case is somewhat academic since there will be no real valued r unless 

(2Kw/n )<1 . 

O 

Note that in this example the flow has periodic negative portions that 
could be then interpreted as investments for survival of the venture. 
This interpretation was discussed previously. 


EXAMPLE NO. 7 

A potentially usable flow model is as follows: 
f(t) = n Q e at cos wt - K6(t) . 

Here a is again the force of venture mortality. If w = 0 then one 
should obtain for r the expression 


r = (n /K) - a 
o 

and if a = 0 then one should obtain the results of the previous example. 
The condition after discounting is as follows: 

K = Cn 0 (r+a)]/[(r+a) 2 + w 2 ] 

This is a quadratic equation in r yielding two roots. The correct rel- 
evant root satisfying the above mentioned limiting conditions is as 
follows : 

r = [(n Q /2K) - al + /(n Q /2K) 2 - w 2 

For a real valued r it is necessary that (n Q /2K) 2 >w2 e Figure A-8 il- 
lustrates the general shape for this kind of a flow with oscillations 
and negative portions, which, again should be interpreted as sacrifices 
or investments for survival. 



It is appropriate now to illustrate how to treat the discounting 
of ’’sample data” flows or discrete periodic flows by using the principles 
of continuous discounting. The first such example is illustrated in 
Figure A- 9 

ANS(t) ANA(t-T) AN6(t-2T) AN6(t-3T)AN6(t-4T) AN6(t-5T) AN6(t-6T) 


T { 

1 1 

» i 

' 

> 4 

► i 

> 

^ 

j— ■ ■ — — - — — - — — ^ 

1 T 2T 3T 4T 5T 6T 

4 


-K6(t ) 


FIGURE A- 9 


In the above illustration lump sum net returns occur uniformly periodically 
as shown. They are impulse functions AN6(t - nT), n = 1,2,3, ... . Then, 
using the previous principles of discounting, one obtains the following 
condition: 

00 

K - l AN e“ mT = AN/(l - e“ rT ) . 
n=0 

Let AN/T = n Q . If 0<rT<<l, and so e -r ^=l - rT, then r = n 0 /K. For the 
above case one can solve for r yielding the following result: 

r = (l/T ) In [1/(1 -AN/K)] . 


EXAMPLE NO. 9 


If the net returns dollar pulses are attenuated by a factor g n , 
o< g<l , then these pulses have the form 


Then 


and 


AN g n 6(t - nT) , n = 0,1, 2, 3, 

0<ge _rT <l . 

- CO 

K = l m g 1 e^ 7 = AN/(l-ge~ rT ) 
n=0 

r = (l/T) In [g/(l - AN/K)] . 


EXAMPLE NO. 10 


One can introcuce the force of venture mortality, a, and the prob- 
ability of product survival exp (-at) as follows: 

AN e“ anT 6(t - nT) ; n = 1,2, 3, 4, 

For such a train of net returns pulses one obtains the following results 
of discounting: 

Cl = l AN e‘ (r+a)nT = N/(l-e“ (r+a)T ) 

) n=0 

(r = (l/T) In [1/(1 - AN/K)] - a 

If now AN«K, and if AN/T = n then 

o 

r = (n Q /K) - a 

approximating thus the case with the continuous flow n(t) = n Q e _ai: . At 
this point it should be noted that the continuous flows can be approxi- 
mated by sample data discrete flows or vice versa. 
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EXAMPLE NO. 11 


Consider now the net returns pulses of the following form: 

AN n g n 6(t - nT); n = 1,2, 3, 4,. 

In this case 

K = £ AN n g 11 e -rn T = (AN g e -I *^)/(l - g e“ r> ^)^ 

n=0 

0<ge~ r> '^ < 1 . 

Solving for r yields the following expression: 

r = (l/T) In Eg/(1 - (AN/2K) ( v'l + 4K/AN - 1)] . 

If 0<AN«K, then for g=l this sample data discrete case approximates 
a continuous case with the net returns cash flow of the form 

n(t) = (n Q /T) t . 

BENEFIT-SACRIFICE STREAM MODELING: A GENERALIZATION 


The previous B-S flow models for ventures did not take into account 
the fact that in many cases the assets for a venture cannot be obtained 
at an instant of time but rather the investment activity must occur over 
extended periods of time. The general B-S flow has the form 

f(t) = nCt) - i(t) + jS(t-T) 

as pointed out previously. Figure A- 2 provides an illustration of a rather 
typical B-S flow cycle. Its nature has been discussed in some detail. 

It is the purpose of this section to present several B-S flow models of 
this more general type. 


EXAMPLE NO. 12 


Figure A-io illustrates a case where a constant outflow of invest- 
ment lasts for a period of time T whereafter a constant net returns flow 
will take place over another time period T. 
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Then the discounting relationship ior the total cash flow is as 
follows : 


/ f(t) e rt dt 

0 


T 2T 

/ (-i o )e -rt dt + / n Q e -rt dt = 0 


This relationship yields the following equivalent condition: 

(e" 1 ^") 2 -2[(n + i Q )/2n 0 ] e -rT + (i Q /n 0 ) = 0 . 

Thus 

e rT = (2n 0 )/[(n 0 + i) - (n 3 - i c )] . 

In a limiting case where rT«l one must have e r ^~l+rT. If the - sign is 
used one gets 

r - (l/T)[(n 0 -i 0 )/i 0 !l 

whereas if + sign is used this would imply that r=0 which is not of any 
relevant Interest. Therefore - sign is the appropriate one. Further, 


r = (l/T) In (n 0 /i o ) = (1/T) In (N/K) 


It is noted again how the rate of return on investment is inversely 
proportional to the investment duration or lag time T while it is only 
logarithmically proportional to the relevant B-S ratio. For a given 
N/K ratio it is desirable to organize the program in such a fashion that 
T is minimized. This again is related to basic realization of the pro- 
gram’s goals in a minimum time. It is, indeed, possible to illustrate 
that the maximization of the rate of return on investment may require 
a well defined restriction for the financial ration N/K, i.e. the real 
performance for a program is of such an essential nature that the B-S 
ration must be matched for it in order to maximize the rate of return 
on investment. This demonstration of the principle can be done as 
follows: A short duration T or a large inverse duration 1/T cannot be 

obtained at no cost in a typical modern economy. Let s be some non-negative 
real number. Assume that 1/T is proportional to K raised to this number 


s, i.e. 

l/TsdT ; s>0 . 


Then 


r«K S In (N/K) . 


If N is estimated or given, then the only manipulated factor is K, Then 
r can be maximized with respect to K. The necessary condition is as 
follows: . 

dr/dk«sK s " 1 ln(N/K) + K s (-K x ) * o 
or 

N/K 


- e 1/s 


If s is a positive number much less than unity, the B-S ratio N/K for 
the maximum return on investment can be high. However, if s becomes 
larger than unity, then this benefit-sacrifice ratio must be reduced 
for the realization of a maximum rate of return on investment. Indeed, 
it is illustrative to plot the function K s ln(N/K) for a given N and few 
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different s. Choose N = 1,000,000 dollars. Let s take values 1/2, 1 
and 2. Figure A- 11 illustrates the results for these cases indicating 
the limited role of benefit-sacrifice ratios for best results. 

One might think then that in this illustration the cost of obtain- 
ing the shortened time T for an increased rate of return on investment 
is overriding the inportance of the B-S ratio N/K and that the control 
of such a cost of reduced time T is not a problem of, say, financial 
management as much as it is the problem of the management responsible 
for realisation of actual performance for the goals of the program. 


EXAMPLE NO. 13 


If the problem remains the same as the previous one except in that 
now the net returns flow is n(t) = n Q for all t>T, then one finds that 
the rate of return on investment is given by the following expression: 

r = (1/T) ln[l + (n 0 /i 0 )] . 

Again the time period T is the duration of the investment activity or 
the investment lag time. Thus K = i Q T. Figure A-12 illustrates this 
case ♦ 



EXAMPLE NO. 14 


Considerable mathematical simplifications can be obtained by con- 
structing B-S models losing appropriate superposition of relatively 
simple component functions which, nevertheless, could be characterized 
by relevant parameters for the discription of a program. Consider, for 
example, a case where the force of venture mortality is not a factor but 
where the lag in realizing productive assets by a process of investment 
is controlled by the force or rate of renovation or innovation. If such 
a rate is low, one would expect long lag times. What kind of a model 
jj could be constructed to portray this situation with a reasonable 


£ 

£ 
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simplicity? In the most simple case one might assume the force of in- 
novation or rate of innovation is just constant, say, b. The probability 
that the job will not be realized within a given time would be then an 
exponentially decaying one in time, i.e. exp(-bt). This concept will be 
now used to develop a B-S flow model for an innovation limited venture 
facing no product mortality. 

It is assumed that the venture is characterized by a one way proc- 
ess from a state of no business and building up to the state of steady 
business. Figure A-13 illustrates the model for the cash flow of this 
venture . 


f(t) 



For the given total cash flow the rate of return on investment is 
simply as follows: 

r = b 

Note then that this rate of return on investment is directly proportional 
to the rate of innovation b. If b=0 then r=0. The brake even period T Q 
(or zero crossing time) occurs when f(t) = 0, i.e. is given by the fol- 
lowing expression: 

T q = (1/b) In [1 + (i 0 /n Q )] . 

Thus the time T Q to the brake even condition (f(t)=0) is inversely pro- 
portional to the rate of innovation . Since in this case 1/b is the mean 
time of realizing the productive asset by the investment process, then 
T Q is directly proportional to this mean realization time of the pro- 
ductive asset for the venture . The total dollar volume of the investment 
is given as follows: 

T T 

K = / °[-f(t)]dt = / °C(n +i )e _bt - nldt = C(n 0 +i 0 )/b](l-e" bT °)-n 0 T 0 . 
0 0 ° 0 
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Noting the expression for T 0 one has e~^° = n 0 /(n 0 +i 0 ). Thus 
K = i Q /b - n 0 T Q = (l/b)[i Q - n D In (l + i Q /n 0 )] 

= (i G /b)[l - (n 0 /i Q ) In [1 + (i 0 /n 0 )]] 

As an illustration, if i 0 /n 0 =10, then In 11 = 2.4. Thus k=0.76(i o /b). 

Thus far, an earlier model illustrated a venture with no force of 
innovation acquiring an asset at the time t=0, worth K dollars. It was 
subjected to a constant force of mortality, a. Its net B-S flow was 
exponentially decaying and the consequent rate of return on investment 
was expressed as follows: 


r - (n Q /K) - a 

In the last example there were no such force of mortality present. The 
only limitation was the rate of innovation or force of innovation. For 
this case one had the following expression for r: 

r = (n 0 /i Q ) b . 

It would be, indeed, interesting to develop a model where a productive 
asset is built up in an innovation limited way while the output of the 
venture would be subjected to a force of mortality without subsequent 
renovation processes. 


EXAMPLE NO. 15 


Consider a venture which starts in a state of being conceived as an 
opportunity for investment in productive assets. The process of invest- 
ment would then take this venture from its state of being conceived to 
the state where it has realized the planned productive assets. Mean- 
while there is a force of mortality to which force this venture is 
subjected. 

A relatively simple B-S flow model for this kind of situation can 
be constructed by a superposition of two exponential functions , one 
portraying the investment process with a constant force or rate of 
innovation, and the other portraying the process of product mortality 
with a constant force of mortality. Figure A-14 illustrates this case. 
For a meaningful venture "rat race" it is assumed that the rate of 
innovation b is magnificantly greater than the rate or force of pro- 
duct mortality a. 
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0 



In this case one interprets the investment and net returns flows as 
follows: 

i(t) = -f(t) = (n o +i o )e -i:,t -n 0 e~ at for 0£t<T o 
= 0 elsewhere. 

n(t) = +f(t) = n e"kt_ (n_+i )e"'^ )T for t>T 0 
= 0 elsewhere. 

For this case the zero crossing time or time to the break even con 
dition occurs when f(t) = 0 which condition yields the following ex- 
pression for T : 

T q = [1/ (b-a) ] In [1 + (i 0 /n Q )3 

The rate of return on investment in this case is obtained by the method 
discussed previously, and is as follows: 


r = (n 0 /i 0 )[b - a] - a 

With this result one obtains the following interpretation for (b - a): 
b - a is the net rate of innovation in an environment of a continuously 
existing force of mortality a. It is not enough that a venture provides 
just some innovation rate in realizing its new productive existence. 

It must provide a rate of innovation in net over the existing rate of 
mortality of its products and methods of production and marketing, etc. 

For this model the total investment dollar volume K is as follows: 

T T 

K = / °i(t) dt = / °[(n 0 +i 0 )e~ bt - 'i 0 e -at ]dt = [(n 0 +i 0 )/b](l-e“ bT °) 

0 0 


- (n 0 /a)(l-e~ aT °). 


The total net returns dollar volume N is as follows: 


H = / n(t ) dt = / [n Q e~ -(n 0 +i 0 )e~ bt ]dt = (n 0 /a)e“ a °- [(n Q +i 0 )/b]e bTo . 


Noting the expression for T one finds (^“ a ^° = [l+(i 0 /no)] a/ ^ b 


e' bT ° = [l+(i 0 /n 0 )]“ b/(b_a) 


Therefore, K and N as well as N/K can be computed readily if i 0 ,n 0 , a and 
b are specified, a four parameter model. There are two financial para- 
meters and two "performance” parameters. 


The descriptive flexibility and power of parametric flow models can 
be increased by increasing the number of usable relevant parameters* 
Although the above may be of interest in caricaturing some ventures, 
it may be far too simple for others. Therefore it may be desirable to 
introduce more complex models. 


One way to increase the descriptive power of simple models is to 
add into a model a sufficient number of independent parameters or "degrees 
of freedom" which allow one to match such a model to various empirical 
or conceptual situations with appropriate levers for a reasonably good 
fit so that a model, indeed, serves as an adequate MAP for its intended 
purposes of representation* 


EXAMPLE NO* 16 


As an exarple of a multiple parameter flow model, consider a par- 
ticular example where 

i(t) = i 0 (t/T) k "V bt 
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and 


n(t) = 


M q (t/0) m_1 e' at -(t/T) R_ V bt 


used in a particular case with b_>2a. If one chooses 0=T=1 unit of time 
(years), then 

k = / i(t) dt = Cr(k)/T k_ 1 b k ] i 0 
0 

co 

N = / n(t) dt = [r( J n)/6 m - 1 a m : n 


The discounting equation is as follows: L[i(t)!l = L[n(t)] ; thus 

(r+b) k /(r+a) m _= [i + (i 0 /n c )] [r(k)/r(m)3 [o m_1 /T k ” 1 :] = [r(k)/r(m)3 

[0 Tn ~ 1 /T k>1 ] £L + [(r(m)T k_1 b k /r(k)e tn_1 a m ) - 1](K/N)^ = a . 

A subclass of such potentially useful flow models is obtained if m = 2k. 
In such a case one has the following relationship: 

(rtb) k /(r+a) 2k = A or (r+b) = A 1//k (r+a) 2 . 

This quadratic equation has two roots one of which is relevant since it 
satisfies the necessary condition r>0 if N/K>1: 


(1/2) A~ 1/k [1 + A + a (b - a) A 1 ^] - a . 


For this model one still has the following relationships so long 

b>2a: . . 

i = [ (k-1 ) /ebT] K_1 i_ at t = T- = (k-l)/b 
max ° 1 

n max = ^ m - 1 Veae] m-1 n o at t = T n = (m-l)/a . 

If one lets T = 0 = 1 and m = 2k, then the zero crossing or break even 
time T 0 is as follows: 

[k/(b-a)] In T 0 t T 0 = [l/(b-a)3 In [1 + (i 0 /n 0 )L 

Consider a particular case for which T=0=l;b=l,a= 1/4, 
k = 3, m = 2k - 6, K = 1,000,000 $, N = 4,000,000 $. For this case 

i 0 = 500,000 $/year 

i = 275,000 $/year at T. = 2 years 

n Q = 8.15 $/year 

n mak 186,000 $/year at T fi = 20 years 

The zero crossing time occurs when 4 In T Q + T Q = 14.65 years. This 
would yield T Q around 7 jears. Figure A -15 is a plot for this case. The 
points C correspond to f(t), A to i(t), and B to n(t). Thus f(t) 
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corresponds very closely to i(t) and n(t) except in a relatively small 
neighborhood of the zero crossing. 


Note k need not be an integer. However if it is an integer while 
m = 2k, N/K = 4, a = 1/4 inverse years and b = 1 inverse year, then 


[(k-l)!/(2k-l)n )l + [ (2k-l) !4 2k /(k-l) ! - l](l/4) 


A- 


4 2k-l 


This approximation is the better the higher k becomes. Then for this 
example one has the following expression for the rate of return on in- 
vestment: 

r = (l/2)(l/4^ 2J< - 1)/k )i:i + 1 + 3*4^k-l)/k] _ (1/4) 


In this case one obtains the following time parameters, and the values 
for the rate of return on investment: 


k 

%r 

T i 

T n 

1 

32.57 

0 

4 

2 

12.50 

1 

12 

3 

7.69 

2 

20 

4 

5.54 

3 

28 

5 

4.33 

4 

36 

6 

3.56 

5 

44 

7 

3.02 

6 

52 

8 

2.62 

7 

60 

9 

'2.31 

8 

68 

10 

2.07 

9 

76 


The behavior of this cash flow model, as k increases, is illustrated in 
Figure A-16. 

f(t) 



FIGURE a-16 
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It is seen that the general features of this model can Indeed provide a 
relevant portrayal of "typical'’ features of cash flow cycles. 


AN INTRODUCTION TO RISK ANALYSIS UTILIZING PARAMETRIC BENEFIT-SACRIFICE 
FLOW MODELS 


An extension of the parametric B-S flow models is into the areas of 
risk analysis. Each parameter in the cash flow model can be treated as 
a random variable to which one can assign some appropriate probability 
density function. In particular, one could consider a set of parameters 
as a set of mutually independent random variables if this seems desirable. 
The general idea is to find the mean and variance of a general performance 
measure such as the rate of return on investment or pay-back period in 
terms of the statistical measures of the relevant parameters of the cash 
flow. In order to illustrate a variety of cases , four examples will be 
given. 


EXAMPLE NO. 17 


In the Example No. 3 one was considering a persent value model where 
the net returns flow was as follows: 


n(t) = h Q e“ a ^ . 

For this case the rate of return on investment with a total dollar vol- 
ume outlay of K dollars was as follows: 

r = (n 0 /K) - a = [(N/K)-l3 a. 

Let us assume a decision maker has come up with information to his 
satisfaction allowing time to assign from the past and present accepta- 
ble experiences a probability density distribution for the random variables 
n 0 and a. For an i litis trati on, assume K is given with certainty by the 
situation facing a decision maker. Assume, further, n^ and a are two 
mutually statistically independent variables. The decision maker assigns 
the relevant probability density distributions g(n Q ) and s(a) to the ran- 
dom variables n Q and a, respectively 

/ g(n Q )^0 for all n Q >0 s(a)>0 for all a>0 

/ / g(n 0 )dn 0 =1 / /s(a)ds = 1 

L° 

The expected value of r is then as follows: 



CO CO 

E(r) = / / [(n Q /i 0 ) - a]g(n 0 )s(a)dn Q da 

0 0 
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OO CO 

= (1/K) J n Q g(n )dn 0 - / as (a)da - U/K)E(n 0 ) - E(a) . 

0 0 

For the above case one has the following relationships for the expected 
value E(r) and variance Var(r) of r: 

r = (n 0 /K) - a 

E(r) = (l/K)E(n 0 ) - E(a) 

Var(r) = o 2 = (1/K) 2 cr 2 + c 2 

r n Q a 


Let the random variable r lie within the interval E(r) + k. or where 
k^>0 is a real number chosen so that the probability P of r being within 
this interval is a desired one. Then one can define the following ranges 
for r: 


The range of pessimistic rates of returns on investment: 

r<rp g = E(r) - kptfr, kp chosen for pessimistic limit 
The range of neutral rates of return on investment: 

r pe = " k p a r E(r> ) + k o a r = r op 

The range of optimistic rates of return on investment: 


r>r Q p = E(r) + k 0 o r , k o chosen for optimistic limit 

There could be several ways of picking kp and k . The way a decision 
maker wishes to choose these values depends on the particular criteria 
he wishes to emphasize. For example, if k Q + kp = 1, both being non- 
negative , then 

r<Tp e with the probability kp(l-P) 

Tp e <_ r <_ r Q p with the probability P 
r>r Q p with the probability k 0 (l-F). 

With this interpretation, k is the conditional probability that r is in 
the pessimistic range given P it is not within the middle range. Also k Q 
is the conditional probability r is in the optimistic range given it is 
not in the middle or normal range. 


EXAMPLE NO. 18 


In the previous example No. 15 one had a special case for a B-C flow 
characterized by a constant "force of product mortality" a and a constant 
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"rate of innovation" b. It was found in this very special case that the 
rate of return on investment was as follows: 

_ -i' 

r = (n Q /i o )(b - a) - a . 

Consider now a risk analysis where i c is a given deterministic parameter 
while n Q , b , and a are considered as mutually independent random variables. 
With the calculus of means and variances of a set of mutually independent 
random variables one can derive the following basic expressions: 


E(r) = (l/i 0 )E(n 0 )[E(b)-E(a)] - E(a) 

0 2 = E(r 2 ) - [E(r)] 2 

= (l/i n ) 2 ^E(nQ)[E(b 2 )+E(a 2 )]-E(n 0 ) 2 [E(b) 2 +E(a) 2 3-2E(b)E(a)a 2 

+C(2/i o )E(n 0 )+l]a 2 ’. 


Again one can choose k Q and kp which are non-negative and whose sum is 
unity such that one obtains the desired pessimistic, "normal" and 
optimistic ranges for r. 


EXAMPLE NO. 19 


Consider again Example No. 13. In this case the total investment 
dollar volume was K = i c T. In order to get back this money one needs 
to collect n D dollars per year for T-^ years so that i Q T = n 0 T^. Thus 
T-^ = (i 0 /n 0 )T and the pay-back time would be then as follows: 

T pb ■ T + T 1 * t 1 + (n 0 /i o )3 . 

Assume now n G and are random variables while i Q is a given deter- 
ministic parameter. Assume also n Q and are mutually independent. 
Then for the pay-back period one obtains the following mean and vari- 
ance: 


E(T pb ) ’ <i/V E(ll o> E < T i> + E(T 1 ) 

<4 = (l/io)[c 2 o o| 1 + EC^)^ 2 ^ E(n 0 ) 2 a 2 1 l + 


EXAMP T E NO. 20 


In Example No. 3 one found the time price relationship between K 
and N to be as follows: 


K = Ca/(r+a)] N. 
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Suppose a decision maker sets r to be some desired value set, for example 
by capital markets. Eut be feels uncertain about a and N, and there- 
fore he feels uncertain about the required investment dollar volume which 
is clearly a function of the two random variables a and N. He would con- 
ceive some probability density functions for a and N, respectively. 
Assuming a and N are mutually independent random variables one obtains 
the following egressions: 

E(K) = E[a/(r+a)] E(N) 

c K 2 = E(K 2 ) - CE(K)3 2 

= E[a 2 /(r+a) 2 ] E(N 2 ) - [E(a/(a+r)f[E(N)] 2 . 

The ranges of interest may then be as follows: 

The optimistic range for the total capital outlay: 

K < E(K) - k 0 o K 

The "normal 1 ' range for the total capital outlay: 

E(K) - k o 0j^K<_ E(K) + k p c K 
The pessimistic range for the total capital outlay: 

K > E(K) + k a v . 

P K 


EXAMPLES OF MAXIMIZATION OF A DECISION MAKER'S UTILITY UNDER A PATTERN 
OF CONJOINED SIMULTANEOUS BENEFIT-SACRIFICE STREAMS 


The earlier classical utility maximization problems fell in the 
category where the one-dimensional scalar utility measure was considered 
as a response to a pattern of conjoined simultaneous stimuli each of 
which was necessary among others for the generation of a positive re- 
action. It is the purpose of this section to introduce few examples of 
this type. 


• SAMPLE NO. 21 


Consider a "decision maker" who is completely convinced that he 
needs to execute simultaneously two distinctly different programs in 
order to be happy at all. For the program No. 1 he expects a net re- 
turn of benefits N^ at a time price p^. For the program No. 2 he expects 
a net return of benefits N 2 at a time price pg • Both N]_ and 
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N 2 are measured in same units as is the total resources K for the re- 
quired sacrifice. Let Nj_ m and N 2m be the respective saturation values 
for and N 2 , respectively, to which this "decision maker” will be 
sensitive at the maximum. Also, let N]_ 0 and N 20 be the respective norms 
for N, and N 2 . Assume that for a "typical case” Ni 0 <<N im 311(1 N 2 o <<N 2 m* 
Then the utility response relationship to these stimuli is assumed to 
be of the following form: 

ll ll 

U/U 0 = C(N 1 /N lo )/[l + (%/N^)]] 1 C(N 2 /N 2o )/[1 + N 2 /N 2m )]] 2 

Here h. and h 2 are two non-negative dimensionless real numbers character- 
izing the power lav: coefficients. This relationship satisfies the follow- 
ing additional requirements assuming all and N 2 are non-negative real 
measures: 

1. The ordering of the utility U is preserved no matter how the 
positive values U Q , N lm , N lo «N lm , N 2m , N 2o «N 2ni , hi and h2 are 
chosen. Thus this function has the ordinal property in the 
Pareto sense and, respectively, and appropriate structure for 
psycho-physiological scaling for a stimulus-response relation- 
ship . 

2. U/U 0 is a monotonic increasing function of, both, N-^ and N 2 . 

That is, the relative utility will not ever decrease if either 

or N 2 increase, given any positive values of U c , N^ m , N^ o ,N2 ms 
N 2o , h^ and h 2 « 

3. U/U Q has continuous first and second partial derivatives for 
all non-negative values of and/or N^. The total second 
differential will exist for all such non-negative values of 
N. and N 2 in whatsoever direction one wishes to obtain it. 

Thus the necessary properties for the classical marginal 
utility theory are provided. 

4. The saturation property built into this model is equivalent 

to the principle of an eventually diminishing marginal utility. 

The total sacrifice volume K relates to the net benefits and N 2 re- 
ceived from the two programs as follows: 

K = p^ + p 2 N 2 

p 1 and p 2 are the dimensionless time prices of and Nj, respectively. 

The general objective of the optimization process is to maximise U/U c 
under this constraint of available resources. 

This simple case depending on the two conjoined benefit variables 
N. and N 2 and on a sacrifice constraint that the total allowable re- 
sources in the same units as N-^ and N 2 amounts to K can be now general- 
ized to a finite number of benefit variables N-^, N 2 , .... N n as follows: 
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Consider a total sacrifice K measured in any appropriate units to 
be allocated for the realization of a conjoined set of n simultaneous 
programs each with its respective total net benefit N,, Ng* . ..., N n 
where all these benefits are measured in the same units as K. Let p^, 
p 2 , ... , pn be the respective time prices of the above benefits. 

Then the problem is to maximize 
n 

U/U Q = it C(N 1 /N 1q )/[1 + 
i=l 

under the constraint 

K = l p.N. . 

1 = 1 

This kind of a problem is best solved by the method of Lagrange’s 
multipliers. Introduce first the modified function 

n 

V(N l5 N 2 ,...,N n ) = U/U 0 -A[ l p.N. - K] 

i=l 

The necessary first order conditions for maximum value are as follows: 
aV/BNi = - A Pi = 0 ; i = l, 2, 3, ..., n ; 

n 

3V/3A = l p.N. - K = 0 . 

1=1 

Actually the second order conditions are satisfied autonrtically by the 
way the function U/U D was constructed, but they can be assured formally by 
investigating the sign of the bordered Hessian determinants will alter- 
nate properly. It follows then from the first order conditions that 

(l/ppHau/aN^ = ci/p 2 )H3U/3N 2 ] = = (i/p n )Csu/aN n ]. 

One may now consider a particular illustration how some individual 
might optimize the allocation of sacrifice K between two conjoined 
simultaneous programs characterized by the following two sacrifice- 
benefit flows: 
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Assume now saturation effects are not important. Thus, consider the 
following utility function for and N 2 : 

U/U 0 = CN 1 /N lo 3 h EN 2 /N 2o ] 1 - h . 

The rates of returns on sacrifices and the respective time price relation- 
ships for the two programs to be executed conjointly are, respectively, 
as follows: 

- £ = a 1 [(N 1 /K 1 )-l] or = [a 1 /(r 1 +a 1 )l N-j_ . 

Pi = a l/( r i +a l)* 

r 2 = a 2 [(N 2 /K 2 ) 1/2 -l] or I< 2 = [a 2 /(r 2 +a 2 )] 2 N 2 . 

O 

P 2 = [a 2 /(r 2 +a 2 )3 . 

Then one has further the following constraint: 

K = + K 2 = p^N^ + p 2 W 2 = Ea 1 /(r 1 +a i )]N 1 +[a 2 /Cr 2 +a 2 )] 2 N2 . 

The first order condition yields the following relationship for the 
maximum relative utility: 

(l/ Pl )[ 311/9^] = (l/p 2 )[8U/3N 2 3 or 

E(r 1 +a 1 )/a 1 ] h Ni' h “ 1 N 2 1-h = [(r 2 +a 2 )/a 2 3 2 (l-h) N^ 1 or 

h K 2 = (1 - h) K ± or K x = h K and K 2 = (l-h) K 

This represents the optimal splitting of resources between the two con- 
joined programs. Since 

N 2 = hK/p x and N 2 = (l-h)K/p 2 , 

the maximum utility would have the following relative value: 

U/U q = [N 1q N 2 ^ 1 [(r 1 +a 1 )/a 1 3 h [(r 2 +a 2 ] 2 *' 1 . 

One notes the higher r^ and r„ becomes and the smaller a^ and a 2 become 
the higher the relative utility, and a 2 represent venture failure 
rates. If h = 1/2, then 

U/U G = (N lo N 2o ) 1 ^ 2 [(r 1 +a 1 )/a 1 ] 1/,2 [(r 2 +a 2 )/a 2 ] K . 


tr~i 

■' i ,‘i 

a 

1 

1 

I 


CONOLUSIONS 

The several examples presented here were done in order to illustrate 
how parametric modeling of benefit-sacrifice flows works. Each para- 
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meter introduced into the flow model has a relevant interpretation, and 
can be treated also as a variable. The number of parameters is also a 
number of degrees of freedom in shaping benefit and sacrifice streams 
to fit data or provide an adequately complete portrayal of a situation 
involving several performance, time, and financial factors. Benefit- 
Sacrifice streams relate also to product institutional life cycle proc- 
esses not specifically discussed here. Some references on life-cylce 
modeling appear in Appendix M. 


APPENDIX B 


Definitions And Interest Factor Notation Used In 
Engineering Economic Calculations 


TABLE 1 


DEFINITIONS AND SYMBOLS USED FOR PARAMETERS 
(Suggested Standards) 

No. Definition of Parameter Symbol 


1. Effective interest rate per interest period. i 

2. Nominal interest rate per year. r 

3. Number of compounding periods. N 

4. Number of compounding periods per year. M 

5. Present sum of money. The letter "P" implies 

present, (or equivalent present value) P 

6. Future sum of money. The letter "F" implies 

future, (or equivalent future value) F 

7. End-of-period cash flows (or equivalent end-of- 
period values) in a uniform series continuing 
•For * specified number of periods. The letter 

’’A" implies annual or annuity. A 

8. Uniform period-by-period increase or decrease 

in cash flows (or equivalent values); the 
arithmetic gradient. G 

9 Amount of money (or equivalent value) flowing 

continuously and uniformly during a given ___ _ 

period. P or F 

10. Amount of money (or equivalent value) flowing 

continuously and uniformly during each and 
every period continuing for a specific number 
of periods. A 

Appendix B was obtained from [13C] 
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TABLE 2 



Name of Factor 


MNEMONIC/FUNCTIONAL FORMS 
OF 

COMPOUND INTEREST FACTORS 
(Suggested Standards) 


Mnemonic 

Format 


Functional 

Format 


Group I. All cash flows discrete; end-of-period compounding 




1. Compound Amount Factor 
(Single Payment) 

2. Present Worth Factor 
(Single Payment) 

3. Sinking Fund Factor 

4. Capital Recovery Factor 

5. Compound Amount Factor 
(Uniform Series) 

6. Present Worth Factor 
(Uniform Series) 

7. Arithmetic Gradient Conversion 
(to Uniform Series) 

8. Arithmetic Gradient Conversion 
(to Present Value) 


(CA-i%-N) (F/P,i%,N) 

(PW-i%-N) (P/F , i% f N) 
(SF-i%-N) (A/F,i%,N) 

(CR-i%-N) (A/P,i%,N) 

(SCA-i%-N) (F/A,i% f N) 

(SPW-i%-N) (P/A,i%,N) 

Factor 

(GUS-i%^*N) (A/G , i% ,N) 

Factor 

(GPW-i%-N) (P/G,i%,N) 


Group II. All cash flows discrete: continuous compoundinc 


9 . 

U 

■lE 

wt) 



Continuous Compounding 
Compound Amount Factor 
(Single Payment) 

Continuous Compounding 
Present Worth Factor 
(Single Payment) 

Continuous Compounding 
Sinking Fund Factor 

Continuous Compounding 
Capital Recovery Factor 

Continuous Compounding 
Compound Amount Factor 
(Uniform Series) 

Continuous Compounding 
Present Worth Factor 
(Uniform Series) 


(CCA-r % -N ) ( F/P , r % , N ) 

(CPW-rS-N) (P/F,r%,N) 
(CSF-r%-N) (A/F,r%,N) 
(CCR-r%-N) (A/P , r % , N) 

(CSCA-r%-N) (F/A,r%,N) 

(CSPW-r%-N) (P/A,r%,N) 
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Ref. 

No . Name of Factor 


TABLE 2 Continued 


Mnemonic 

Format 


Functional 

Format 


Group III. Continuous, uniform cash flows: continuous compounding 

(payments during one period only) 

15. Continuous Compounding 

Present Worth Factor 

(single, continuous payment) (CPW-i%-N) (P/F,i%,N) 

16. Continuous Compounding 

Compound Amount Factor __ 

(single, continuous payment) (CCA-i%-N) (F/P,i%,N) 

Group IV. Continuous, uniform cash flows: continuous compounding 

(payments during a continuous series of periods ) '• 

17. Continuous Compounding 

Sinking Fund Factor 

(continuous, uniform payments) (CSF-i%-N) (A/F,i%,N) 

18. Continuous Compounding 

Capital Recovery Factor __ 

(continuous, uniform payments (CCR-i%-N) (A/P,i%,N) 

19. Continuous Compounding 

Present Worth Factor _____ „ 

(continuous, uniform payments) (CSCA-i%-N) (F/A,i%,N) 

20. Continuous Compounding „ 

Present Worth Factor 

(continuous, uniform payments) (CSPW-i%-N) (P/A,i% ,N) 
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TABLE 3 


DIAGRAMS,, ALGEBRAIC FORMS, AND USES FOR COMPOUND INTEREST FACTORS 
(Explanatory Supplement to Table 2) 

Group I. All cash flows discrete: end-of-period compounding 
Cash Flow diagram for factors 1 through 6 (and 9 through 14) s 


No, 

©" 

© 

© 

© 

© 

© 


Cash 




1 

2 

3 

• ■ • 

N-l 

N 

1 J 

r 

* \ 

3 



Name of Factor 


Algebraic Form 


Compound Amount Factor 
(Single Payment) 

Present Worth Factor 
(Single Payment) 

Sinking Fund Factor 


Capital Recovery Factor 


Compound Amount Factor 
(Uniform Series) 


(l+i) N 

(l+i)“ N 

i 

(l+i) 1 

i(l+i) M 
(l+i) N ~1 

i 


Present Worth Factor 
(Uniform Series) 

flow diagram for factors Q) 


and 


(l+i) — 1 
i (l+i) ^ 




Use when; 

Given P, to find F 
Given F, to find P 
Given F, to find A 

Given P, to find A 

Given A, to find F 
Given A, to find P 



Arithmetic Gradient Conver- 
sion Factor (to uniform series 

Arithmetic Gradient Conver- 
sion. Factor (to present value) fjl. 

Li 



N I Given G , to find A 
(l+i) — 1 _ 


— V-1 

(i+i) -lj Li (l+i) 




Given G*to 
find P 
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Group II. 
No. 


All cash flows discrete: continuous compounding 


© 

0 . 


Name of Factor 

Continuous Compounding 
Compound Amount Factor 
(Single Payment) 

Continuous Compounding 
Present Worth Factor 
(Single Payment) 

Continuous Compounding 
Sinking Fund Factor 

Continuous Compounding 
Capital Recovery Factor 

Continuous Compounding 
Compound Amount Factor 
(Uniform Series) 

Continuous Compounding 
Present Worth Factor 
(Uniform Series) 


Use when; 


rN 

e - 

rN, 


rN 

1 a 

'rN 


e 

“rH 


r N 

Given 

P, 

to 

find p 

-rN 

i 

Given 

F, 

to 

find p 

‘1 

Given 

F, 

to 

find a 

•1 

Given 

P, 

to 

find a 

■i 

■i 

Given 

A, 

to 

find F 

•l 

Li 

Given 

A. 

to 

find P 


(e r -l) 


Group III. Continuous/ uniform cash flows: continuous compounding 

(Payments during one period only) ~~ 


uiUlB 

2* 

i 

B 

2 

* * 

3 

B 

N_1 

N 

[ 3 

mm 


No. 


Name of Factor 


re 


Algebraic Form 

i(Ui)“ ! 
£n(l+i) 


r . ...... —N 

e -1 _ l (lj-i) 


rN 


Use when: 

Given F , to find P 


Continuous Compounding 
Present Worth Factor 
(Single Continuous 
Payment) - 

irN 3T * * X vow 

Continuous Compounding a (e -l)_i(l+i) _ Given P, to find F 


Compound Amount Factor 
(Single Continuous 
Payment) 


re 


£n(l+i) 
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1 

I 

I 

I 

I 

1 


Prnun IV. Continuous, uniform cash flows 
~ (Payments during a continuous series of periods) 

Cash flow diagram for factors 17 through 20 


1 

2 

3 

• • • 

N -1 

N 

m 

IP 



BMP 

W 

l A l 


P 


T 

dJ 


gn 




% 


Jiaa. 
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Name of Factor 

Continuous Compounding 
Sinking Fund Factor 
(Continuous, Uniform 
Payments) 

Continuous Compounding 
Capital Recovery Factor 
(Continuous, Uniform 
Payments) 


Algebraic Form 

r _ &n(l+i) 
e rN -l (1+1) N -1 


Use whetn 

Given F, to find A 


„ * Given P, to find A 

re rN _ (l+i) N &n (1+i) 


e rN -l 


(l+i) N -l 


Continuous Compounding 
Compound Amount Factor * e rN_. 
(Continuous, Uniform — - — 

Payments) 


Given a, to find F 

(l+i) N -l 

£n(l+i) 


Continuous Compounding 
Present Worth Factor e r ^-l 

(Continuous, Uniform — — ^ 

Payments) re 


Given A, to find P 
(l+i) N -l 
(l+i) N in(l+i) 


% 
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APPENDIX C 


Interest Formula Values 
Discrete Compounding 
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Single Payment Compound Amount Factor 
Discrete Compounding 
(F/P, r%, N) 



i% 

2% 

4% 

6% 

8% 

10% 

12% 

15% 

20% 

25% 



To rindF 

To find F 

To find F 

To find F 

To find F 

To find F 

To fir i F 

To find F 

To find F 

To find F 



Given P 

Given P 

Given P 

Given P 

Given P 

Given P 

Given P 

Gii cn P 

C : -ven P 

Given P 


* 

F/P 

FtP 

F/P 

F/r 

NIP 

F/P 

F/P 

FfP 

F/P 

F/P 

— - * 


i 

I.OIGO 

1.0200 

1.0400 

1.0600 

1,0800 

1.1000 

1.1200 

1*1500 

1*2000 

I 1.2500 

1 

2 

1,0201 

1.0404 

1.0816 

1.1236 

1.1664 

1.2100 

1.2544 

1,3225 

1.4400 

1.5625 

2 

3 

1.0301 

1*0612 

1.1249 

1.1910 

1.2597 

* 1.3310 

1.4049 

1.5209 

1.7280 

1.9531 

3 

4 

1 .0400 

1.0824 

1.1699 

1.2625 

1.3605 

1.4641 

1.5735 

1.7490 

2.0736 

2.4414 

4 

„5 

1.0510 

1.1041 

L2162^_ 

1.3332 ... 

1.4693 

1.6105 

1*7623 

2.0114 

2.4883 

3.05 IS 

5 

6 

1,061 S 

1.1262 

1.2653 

1.4185 

1.5869 

1.7716 

1.9738 

‘ 2.3131 

2,9860 

— mvr - 

6 ! 

7 

1.0721 

1.1487 

1.3159 

1.5036 

1.7138 

1.9487 

2.2107 

2.6600 

3.5832 

4.7684 

7 1 

8 

1,0829 

1.1717 

1.3686 

1.5938 

1.8509 

' 2.1436 

2.4760 

3.0590 

4.2998 

5.9605 

8 

9 

1JW37 

1.1951 

1.4233 

1.6895 

1.9990 

2.3579 

2.7731 

3.5179 

5.1598 

7.4506 

9 

J£L_ 

-UStt - 

i-2>2a 

— 1-480?, . . 
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2.5937 

3.1058 

4.0456 

6.1917 

9.3132 

in 

u 

1.1157 

1.2434 

1.5395 

1.8933 

2.3316 

2.8531 

3.4785 

4.6524 

7.4301 

11.6415 

U 

12 

1.1268 

1.2682 

1.6010 

2.0122 

2.5182 

3.1384 

3.8960 

5.3502 

8.9161 

14.5519 

12 

13 

1.1381 

1,2936 

1.0651 

2.1329 

2.7196 

3.4523 

4.3635 

6.1528 

10.6993 

18.1899 

13 

14 

1.1495 

1.3195 

1.7317 

2.2609 

2.9372 

3.7975 

4.8871 

7.0757 

12.8392 

22.7374 

14 

13 

1.1610 

1.3459 

1.8009 

2.3966 

3.1722 

4,1772 

5.4736 

8.1371 

15.4070 

28.4217 

15 

16 

1.1726 

1,3728 

1.8730 

15404 

3.4259 

4.5950 

6.1304 

9.3576 

♦ 18.4884 

35.5271 

16 

17 

1.1843 

1.4002 

1.9479 

2.6928 

3,7000 

5.0545 

6.8660 

10.7613 

22.1H6’, 

44.4089 

17 

18 

1,1961 

1.4282 

2.0258 

2.8543 

3.9960 

5.5599 

7.6900 

12.3755 

26.6233 

55.5112 

18 

19 

1.2081 

1.4568 

2.1068 

3.0256 

4,3157 

6.1 ! 59 

8.6128 

14,2318 

31.9480 

69.3889 

19 

20 

1.2202 

1.4859 

--2.1911 1T .__ 

3.2071 

4.6610 

6.7275 

9.G463 

16,3665 

38.3376 

86.7362 

20 

2i 

1.2J24 

1.5157 

2.2783 

3,3996 

5.0338 

7.4002 

10,8038 

" 18.1215 

46.0051 

— mm — 

tC 

22 

1.2447 

1.5460 

2.3699 

3.6035 

5.4365 

8,1493 

12,1003 

21.6447 

55.206! 

135.525 

22 , 

23 

1.2572 

1.5769 

2.4647 

3.8197 

5.8715 

8,9543 

13.5523 

24,8915 

66.2474 

169.407 

• 23 • 

24 

1,2697 

1.60*4 

2.5633 

4.0489 

6.3412 

9.8497 

15.1786 

28.6252 

79,4968 

211.758 

i 

24 

25 

1.2824 

1.6406 

— 2.6658 4.2919 

- G-84B5 

10.8347 

17.0001 

32.9189 

95,3962 

264.698 

25 

26 

1.2953 

1.6734 

2.7 725 

4.5494 

7.3964 

11.9182 

19.0401 

37.8568 

i 14,475 

1 330.552 "*■ 

26 

27 

1.3082 

1.7069 

2.8834 

4,8223 

7.9881 

13.1100 

21.3249 

43.5353 

137*371 

413,590 

27 

28 

U213 

1.7410 

2.9987 

5.1117 

8.6271 

14.4210 

23.8839 

50.0656 

164.845 

516.988 

28 

29 

1.3345 

1.7758 * 

3.1187 

5.4184 

9.3173 

15.8631 

26.7499 

57.5754 

197.814 

646.235 

29 

30 

1.3478 

1.81 14 

- 3.2434 . 

5,7435 

10.0627 

17.4494 

29.9599 

66.2118 

237 376 

807.794 

30 

35 

1.4166 

1.9999 

3.9461 

7.6861 

14.7853 

28.1024 

52,7996 

133.176 

590,668 

"“SSTSCIv™ — 

35 

40 

1.4889 

2.2080 

4.8010 

10.2857 

21.7245 

45,2592 

93.0509 

267.863 

1469.77 

7523.16 1 

40 

45 

1.5648 

2.4379 

5,8412 

13.7646 

31.9204 

72,8904 

163.988 

538,769 

3657*26 

22958.9 j 

45 

50 

1.6446 

2,6916 

7.1067 

18.4201 

46.9016 

117.391 

289.002 

1083.66 

9100.43 

70064.9 1 

50 

55 

1.7285 

2.9717 

8-6464 , 

24.6503 

68.9138 

189.059 

509.320 

2179,62 

22644,8 


55 

60 

1.8167 n 

3,2810 

10.5196 

32.9876 

101.257 

5 304.481 

897*596 

4384,00 

56347.5 


60 

65 

1.9094 

3,6225 

12.7987 

44.1449 

148.780 

490.370 

1581*87 

8817*78 



65 

70 

2.0068 

3.9996 

15.5716 

59.0758 

218,606 

789.746 

2787.80 

17735.7 



70 < 

75 

2.1091 

4.4158 

18.9452 

79.0568 

321.204 

1271.89 

4913.05 

35672*8 



75 ! 

80 

2.2167 

4.8754 

23.0498 . _ 

105.796 

- 471.955 

2048.40 

8658*47 

71750*8 



80 

85 

2.3298 " 

5.3829 

28.0436 

141.579 

693.456 

3298.97 





85 

90 

2.4486 

5,9431 

34.1193 

189.464 

1018.92 ) 

5313.02 
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90 

95 

2.5735 

6.5617 

41.5113 

253.546 

1497.12 

8556.67 





95 

100 

2.7048 

7.2446 

50.5049 

339.301 

2199.76 

13780*6 
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ft 

1% 

To find T 
Given F 

P/F 

2% 

V 

To Tind P 
Given F 
P/F 

4% 

To find P 
Given F 

P/F 

Single Payment Present Worth Factor 
Discrete Compounding 
(P/F, r%, N) 

6% 8% 10% 12% 

To find P ■ To find P To find P To find P 

Given F Giv*r. F Given F Given F 

P/F ; P/F P/F W 

15% 

To find P 
Given F 

P/F 

20% 25% 

To find P f To find F 

Given F Given F 

P/F ! FfF \ 

N 

1 

0.9901 

0.9804 

0.9615 

0.9434 

0.9259 

0.9091 

0.8929 

■ 0.8696 

0.8333 

0.8000 r 

i 

2 

0.9803 

0.9612 

0.9246 

9.8900 

G-8573 

0.8264 

0.7972 

0.7561 

0,6944 

| 0.6400 

2 

3 

0.9706 

0,9423 

0.8390 

0.8396 

0.7938 

0.7513 

0.71 18 

0.6575 

0.5787 

0.5120 . 

3 

4 

0,9610 

0.9238 

0.8548 

0.7921 

0.7350 

0.6830 

0.6355 

0.5718 

0.4823 

1 0,4096. 

4 

5 

0.9515 

0.9057 

. 0.8219 

I _0.7473_ gl , 

0-6806 

0.6209 

0.5674 

0.4972 

0.4019 

0.3277 ; 

5 

6 

0.9420 

0.8 8H0 

0.7903 

0.7050 

0.6302 

0.5645 

0.5066 

”””074323“ — 

0.3349 

0,2621 

— T 

7 

0.9327 

0,8706 

0,7599 

0.6651 

0.5335 

0,5132 

0.4523 

0.3759 

0.2791 

0.2097 

7 

8 

0,9235 

0.8535 

0.7307 

0.6274 

0.5403 

0.4665 

0.4039 

0,3269 

0.2326 

0.1678 

8 

9 

0.9143 

0.8368 

0.7026 

0.5919 

0.5002 

0.4241 

0.3606 

0.2843 

0.1938 

0,1342 

9 

10 

0,9053 

„ O-Hffl-I ii - 




■ 0^855. 0.3220 

0.2472 

0,1615 ^ 

Qlf)74 

10 

11 

0,8963 

0..504J 

0.6496 

0.5268 

0.4289 

0.3505 

0.2875 

0.2149 

0.1346 

0.0859 

ii 

12 

0.8374 

0.7885 

0,6246 

0.4970 

0.3971 

0,3186 

0.2567 

0.1869 

0.1122 

0.0687 

12 

13 

0.8787 

0.7730 

0,6006 

0,4688 

0.3677 

0.2897 

0.2292 

0.1625 

0.0935 

0.0550 

13 

1+ 

0.8700 

0.7579 

0.5775 

0,4423 

0,3405 

0*2633 

0.2046 

0.1413 

0.0779 

' 0.0440 

14 

15 

0.8613 

0.7430 _____ 

- - 0.5553 _ 

^ Q.4173 


0.2394 

- 0.1827 

0.1229 

Q.L1649 


15_ 

16 

0.8528 ^ 

0.7284 

0.5339 

““ 0,3936 

0,2919 

1)2176 

0.1631 

“^0.106$“*“'™ 

0.0541 

0.0281 

trr HT 

17 

0.8444 

0.7142 

0.5134 

0.3714 

0.2703 

0.1978 

0.1456 

0,0929 

0.0451 

0.0225 

17 

13 

0.8360 

0.7002 

0.4936 

0.3503 

0.2502 

0.1799 

0.1300 

0.0808 

0.0376 

0.0180 

18 

19 

0.8277 

0.6864 

0,4746 

0.3305 

0.2317 

0.1635 

0,1161 

0,0703 

0.0313 

, 0.0144 . 

19 


0.8195 

__ 

-^0,4564^ 


H!SI Em 

0.1486 


0.0611 

mmESSmi 

■■iiuiiHsn 

20 

21 


0.6598 

0.4383 

0.2942 

0.1987 

0.1351 

0.0926 

0.053 i " 

0.0217 


21 

22 




0.2775 


0.1228 


0.0462 



22 

23 




0.2618 

0.1703 

0.1117 





23 

24 





0.1577 

0,1015 





24 

25 

0.7798 


HRS STH 

SHBJSJfRBE 

■Kiff < tM 

C.0923 

HBjGSJHI 

0,0304 

BBGJwt&Bi 

HRfljjijtHH 

25 

26 

0.7720 

0.5976 

0-J607 

0.2198 

bkimW ! vm 


0.0525 

wmmTmm 


IHfluSQBBi 

26 

27 




(3.2074 







27 

28 

0.7568 

0.5744 

SEH 

0.1956 

0.1159 

0.0693 
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0,3207 

0,1846 

0.1073 

0.0630 


Tm 



29 


0.7419 

MBfswIH 


0-1741 

mCmSim 

0.0573 
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■KYlffirHi 

30 

35 

6.7059 

o.5mm 

0.2534 

0.1301 

0,0676 

0.0356 

0,0189 

HBififiOBH 


BHSH 

3T s 

40 

0-67' 7 

0.4529 


0.0972 


0.0221 
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45 

0.63 >1 

0.4102 


0.0727 


0.0137 





45 

50 

0.6080 

0.3715 

0.1407 

0.0543 

0.0213 

0.0085 

0,0035 

• ■T . 

MESSuM 
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0,5785 

MmtJkkfawm 
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55 
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• 0,5505 

0.3048 

0.0951 

0.0303 


0.0033 

0.0011 

0.oo02 
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i 

60 

65 

0.5237 

0,2761 
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70 


0.2500 

0.0642 
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0.2265 
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80 
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0.4292 

0.(858 

0.0357 
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85 

90 

0.4084 

0,1683 
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0.0198 



\ 




( . i 















c~* 


: m mm 


azm n 






a*sn3 

Uniform 


ma P-ra mm « 

Series Compound Amount Factor 
Discrete Compounding 
(F/A, rZ, N) 


gCXg«g^j 

cr^™5 L 1 " 1 iu j LV U, ™J BH 


N 

To find F 
Give a A 
F/A 

' To find F 
* Given A 
F/A 

’ To find F 
■ Given A 
| F/A 

To find F 
Given A 
FfA 

To find F 
Given A 
FfA 

To find F 
Given A 
FfA 

\ TofindF 
Given A 
F/A 

“ To find F 
* Given 
, FfA 

To find F 
Given A 
F/A 

. To find F 
GIyctM 
F/A 

N 

1 

1.1)000 

1*0000 

1.0000 

1.0000 

l.oono 

1*0000 

1.0000 

1.0000 

1*0000’ 

1,0000 

1 

2 

2.0100 

2.0200 

2.0400 

2.0600 

2.0800 

2.1000 

2.1200 

2.1500 

2.20U0 

2.2500 

2 

n 3 

3.0301 

3.0604 

3.1216 

3.1836 

3.24 64 

3.3100 

3*3744 

3,4725 

3.6400 

3.8125 

3 

? 4 

4.0604 

4,1216 

4.2465 

4.3746 

4.5061 

4.6410 

4.7793 

4.9934 

5.3630 

5.7656 

4 

5 

5,1010 

5.2040 

5.4163 

5.6371 

5.8666 

6.1051 

6.3528 

6.7424 

7.4416 

8.2070 

5 

6 

6.1520 

6.3061 

6.6330 

6.9753 

7.3359 


8.1152 

8,7537 

9.9299 

11.2588 

6 

7 

7.2135 

7.434J 

7,8983 

8.3938 

8.9228 

9,4872 

10.0890 

11.0668 

12.9159 

15.0735 

7 

8 

8.2857 

8.5830 

9.2142 

9.8975 

10.6366 

11.4359 

12.2997 

13,7268 

16.4991 

19.8419 

8 

9 

9.3685 

9.7546 

10.5828 

11.4913 

12.4876 

13,5795 

14.7757 

16.7858 

20.79B9 

25.8023 

9 

30 

10.4622 

10.9497 

12.0061 

13.1808 

14.4866 

15.9374 

17.5487 

20.3037 

25.9587 

33.2529 

10 

Tl 

1 i.56f.B 

12.1687 

13.4863 

14.9716 

16.6455 

18.5312 ' 

20,6546 

24.3493 

32.1504 

42.5651 

u 

12 

12.6825 

13,4121 

15.0258 

16.8699 

18.9771 

21.3843 

24.1331 

29.0017 

39.5805 

54.2077 

12 

13 

13.8093 

14.6803 

16.6268 

18,8821 

21.4953 

24.5227 

28.0291 

34,3519 

48.4966 

68.7596 

13 

14 

14.9474 

15.9739 

18.2919 

21.0151 

24.2149 

27,9750 

32.3926 

40.5047 

59.1959 

86.9495 

14 

IS 

1 6.0969 

17.2934 

20.0236 

23.2760 

27.1521 

31.7725 

37,2797 

47.5804 

72.0351 

109.687 

15 

16 

17.2578 

18.6393 

21.8245 

25,6725 

30.03 

35.9497 

42.7533 

55.7175 

87.4421 

138.109 

16 

17 . 

18.4304 

20.0121 

23.6975 

28.2129 

33.750*. 

40.5447 

48.8837 

65.0751 

105.931 

173.636 

17 

18 

19.6147 

21.4123 

25.6454 

30.9056 

37.4502 

45.5992 

55.7497 

75.8363 

128.117 

218.045 

16 

19 

20.8109 

22.8405 

27.6712 

33.7600 

41.4463 

51.1591 

63.4397 

88.2118 

154.740 

273.556 

19 

20 

22.0190 

_ ^24.2974 

29.7781 

36.7356 

45.7620 

57.2750 

72.0524 

102.444 

186.683 

342.945 

20 

21 

23.2391 

25.7833 

31,96**2 

39.9927 

50.4229 


81,6987 

118 810 

225.026 

429.681 

21 

22 

24.4715 

27.2990 

34.2480 

43.3923 

55.4567 

71.4027 

92.5026 

137.632 

271.031 

538.101 

22 

22 

25.7162 

28.8449 

36.6179 

46.9958 

60.8933 

79,5430 

104.603 

159.276 

326.237 

673.626 

23 

24 

26,9734 

30.4218 

* 39,0326 

50,8155 

66.7647 

88.4973 

118.155 

184.168 

392.484 

843,033 

24 

25 

28.2431 

-120503 

41.6459 

54.8645 

73.1059 

98.3470 

133.334 

212.793 

471.981 

1054.79 

25 

26 

29.5256 

33.0709 

44.fl 17 

59.1563 

79.9544 

109.182 

150.334 

245.712 

567.377 

1319.49 

26 

27 

30.H208 

35.3443 

47.0842 

63.7057 

87.3507 

121.100 

169.374 

283.569 

681.853 

1650.36 

27 

28 

32.1290 

37.0512 

49.9676 

68.52C1 

95.3388 

134.210 

190.699 

327.304 

819.223 

2063.95 

28 

29 

33.4503 

38.7922 

52,9663 

73.6397 

103.966 

148.631 

214.583 

377.170 

984.068 

2580.94 

29 

30 

34.7848 

40.5081 ... 

56.0349 

79.0581 

113.283 

164.494 

241.333 

434.745 

1181.88 

3227.17 

30 

35 

™ 4 i 66(12 " 

49.9944 

73.6522 

111.435 

172.317 

- — TfTSl 3 “ 

431.663 

881.1/$ 

2948.34 

9856.76 

35 

40 

48.8863 

60.4019 

95.0255 

154.762 

259.0S6 

442.592 

767.091 

1779.09 

7343.85 

30088.7 

40 

45 

56.4809 

71.8927 

121.029 

212.743 

386.506 

718,905 

1358.23 

3585.13 

18281.3 

91831.5 

45 

50 

64.4630 

84.5793 

152.667 

290.336 

573.770 

U 63.91 

2400.02 

7217.71 

45497.2 

280256 

50 

55 

72.8523 

98.5864_ ir __ 

191.159 

394.172 

848.923 

1880.59 

4236.00 

.... 14524.1 _ H3219 _ _ 


55 

60 

" “ 81.6695 

1 14.051 

237.991 

533.128 

i 1253.21 

30K8T” 

7471.63 

29220.0 281732 

r * 

00 

65 

90,9364 

131.126 

294,968 

719.082 

1847.25 

4893*71 

13173.9 

58778.5 



65 

70 

100.676 

149.978 

364.290 

967.931 

2720.08 

7887*47 

23223.3 

118231 1 

• r 


70 

75 ■ 

110.912 

170.792 

448.631 

1300,95 

4002.55 

12708.9 

40933.8 

237812 



75 

80 

121.671 

193.772 

551.245 

1746.60 

5886.93 

20474.0 

72145.6 

478332 



80 

85 

~ D'2.979 

219.144 

676.090 

2342.98 

8655.71 

329^9^ "■*" 


\ 1 



85 

90 

144.86 

247.16 

827*98 

3141*07 

12723.9 

53120.2 

U- , t # - ..<* * 

H * * 1 



90 

95 

157.35 

278.08 

1012.78 

4209*10 

18071.5 

85556.7 
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100 

170.4B 

312.23 

1237,62 

5638.36 

27484.5 

137796 
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Uniform Series Sinking Fund Factor 
Discrete Compounding 
(A/F, r %, N) 


N 

To fiml/1 
Given F 
A/F 

To find A 
Given F 

A/F 

To find4 
Given F 
A/F 

To find/1 
Given F 
A/F 

To findi4 
Given F 
A/F 

To find A 
Given F 
A/F 

To find A 
Given F 
A/F 

To find A 
Given F 
A/F 

To find A 
Given F 
A/F 

To find/4 
Given F 
A/F 

tt 

i 

1X000 

1.0(100 

ToooT 

1X000 

i.onoo 

1.0000 

1.0000 

woST 

1.0000 

1.0000 

i 

2 

0.4975 

0.4950 

0.4902 

0.4854 

0.4808 

0.4762 

0.4717 

0.4651 

0.4545 

0.4444 

2 

3 

0.3300 

0.3268 

0.3203 

0.3141 

0.3080 

0.3021 

0.2963 

0.2880 

0.2747 

0.2623 

3 

4 

0.2463 

0.2426 

0.2355 

0.2286 

0.2219 

0.2155 

0.2092 

0.2003 

0.1 863 

0.1734 

4 

5 

0.1960 

0.1922 

0.1846 

0.1774 

0.1705 

0.1638 

0.1574, 

0.1 483 

0.1344^ 

0.121R_ 

5 

6 

0.1625 

0.1535 

0.1508 

0.1434 

0J3G3 

0.1296 

0.1 12Z 

0.1142 

0.1007 

0.0888 

6 

7 

0.1386 

0.1345 

0.1266 

0.1191 

0.1121 

0.1054 

0X991 

0X904 

0X774 

0.0663 

7 

8 

0.1207 

0.1165 

0.1085 

o.ioio 

0X940 

0X874 

0.0813 

0X729 

0X606 

0X504 

8 

9 

0.1067 

0.1025 

0X945 

0X870 

0.0801 

0X736 

0X677 

0X596 

0.0481 

0.0338 

9 

10 

0.095*6 

0.0913 

0.0333 

0X759 

0X690 

0X627 

0.0570 

0X493 

0X385 

0.0301 

10 

li 

ll.iinf>5 

0,0.322 

r 9.07 4 i 

0X668 

0.0601 

0X540 

0X484 

03)411 

0X311 

0X235 

n 

12 

0.0788 

0.0746 

0.0666 

0.0593 

0.0527 

0.046B 

0X414 

0.0345 

0.0253 

0.0184 

12 

13 

0.0724 

0.0681 

0 X 601 

0.0530 

0.0465 

0X408 

0X357 

0.0291 

0.0206 

0.0145 

13 

14 

0.0669 

0.0626 

0.0547 

0.0476 

0.0413 

0X357 

0.0309 

0.0247 

0X169 

0.0115 

14 

15 

o.nr.2t 

04)578 

0.0499 

0X430 

0.0368 

0.0315 

0X268 

0.0210 

0X139 

0.0091 

15 

16 

0X579 

0X537 


0X390 

0X330 

0.0278 

0X234 

0X179 * 

0X114 

0X072 

16 

17 

0.0543 

0.0500 

0.0422 

0X354 

0X296 

0X247 

0.0205 

0X154 

0X094 

0X058 

17 

18 

0.0510 

0,0467 

0.0390 

0.0324 

0.0267 

0.0219 

0X179 

0.0132 

0X078 

0.0046 

18 

19 

0.0481 

0.0433 

0.0361 

0X296 

0X241 

0X195 

0X158 

0.0313 

0.01)65 

0.0037 

19 

20 

0.0454 

0.0412 

0X336 

0.0272 

0.0219 

0X175 

0X139 

0X098 

0.0054 

0.0029 

20 

21 

0X4 JiT 

0.0333 

nsmr*** 


0.0198 

0.0156 

(1X122 

0X084 

0X044 

0X023 

r ^2* 

22 

0.0409 

0.0366 

0X292 

0.0230 

0.0180 

0X140 

0X103 

0X073 

0X037 

0X019 

22 

23 

0.0389 

0.0347 

0X273 

0.0213 

0.0164 

0X126 

0X096 

0X063 

0.0031 

0X015 

23 

24 

0.0371 

0,0329 

0X256 

0,0197 

0.0150 

0X113 

0X085 

0X054 

0X025 

0X012 

24 

25 

0.0354 

0,0312 

0X240 

0.0182 

0.0137 

0X102 

0X075 

0.0047 

O.O021 

pxoo? __ 

25 

26 

0.0339 

0.U2V7 

0.0226 

0X169 

0X125 

0,0092 

0.0067 

0.0041 

0.0018 

"0.0008 

26 

27 

0.0324 

0.0283 

0X212 

0X157 

0.0114 

0X033 

0.0059 

0.0035 

0.0015 

0.0006 

27 

28 

0.0311 

0.0270 

0X200 

0X146 

0X105 

0X075 

0.0052 

0.0031 

0.0012 

0.0005 

28 

29 

0.0299 

0.0258 

0X189 

0X136 

0X096 

. 0X067 

0.0047 

0.0027 

0.0010 

0.0004 

29 

30 

0.0287 

0.0246 

0X173 

0X126 

0.O08H 

0X061 

_W 4 JL 

0.0023 

_0.0008 

0.0093 

30 

35 

0.0240 

0,0200 

0X136 

0XU9O 

O.OU58 

0.tJ(T37 

0.0023 

0.0011 

0X003 

OX0Q1. 

35 

40 

0.0205 

0.0166 

0X105 

0X065 

0X039 

OOC'23 

0.0013 

0.0006 

0,0001 

* . 

40 

45 

0.0177 

0.0139 

0X033 

0.0047 

0X026 

0.00? * 

0.0007 

0.0003 


• 

45 

50 

0.0155 

0.0118 

0X066 

0X034 

0X017 

0.0009 

0.0004 

0.0001 


* ; 

50 

55 

0.0137 

0.0101 

0X052 

0X025 

0.0012 

0.0005 

0.0002 



i 

55 

61) 

0.0122 

0X088 

0X042 

0X019 

0.0003 

0X003 

0.0001 . 

-k. 

'■ 

; K V 1 . 1 , 1 A 

60 

65 

0.01 10 

0.0076 

0X034 

0X014 

0.0005 

0X002 

* 

' 



65 

70 

C, X )99 

0.0067 

0.0027 

o.ooio 

0.0004 

0X001 


>, 

^ l 


70 

75 

0.0090 

0,0059 

0X022 

0,0008 

0.0002 

* 


i 



75 

M _ 

0.0032 

0.0052 

0X018 

0X006 

0 0002 - 


« 

** 



SO 

85 

0.0075 

0.0046 

0.0015 

0.0004 

0*0001 

* 


! *4 



85 

90 

0.0069 

0X040 

0.0012 

0.0003 

* 

% 


, * - - 



90 

95 

0.0064 

0X036 

0.0010 

0.0002 

* 






95 

100 

0.0059 

0.0032 

0.000B 

.0.0002 

■- 1 

- 





1(H) 


§-&i 


Hrr-nsr* itrtrrr^/» 



To find P 

To find/ 1 

To find P 

To find P 

To find P 


Given A 

Given A 

Given A 

Given/1 

Given A 

N 

P/A 

P/A 

P/A 

P/A 

P/A 

1 

0.9901 

0.9804 

0.%1S 

0.9434 

0,9259 

2 

1,9704 

1.9416 

1.B861 

1.8334 

1,7833 

3 

2.9+10 

2.8839 

2.7751 

2,6730 

2,5771 

4 

3.9020 

3.3077 

3.6299 

3.4651 

3,3121 

5 

4.8534 

4-7135 ' 4.95 IB 

^>•2124 

3,9927 


6.7282 

7.6517 

8.5660 

0.4713 

10 . 3(376 

11.2551 

12.1337 

13.0037 

"74.7178 
15,5622 
16.3082 
17.2260 
18,0455 
18.8570 
19.6603 
20.4558 
21.2434 
iii r 22,023 1 n 
22.7952 
23,5596 
24.3164 
25.0657 
25 .8077 
29.4085 
32.8346 
36:0945 
39,1961 
42.1471 
44.9550 
47,6265 
50.1684 
52.5870 
_ 54.8881 
57.0776 
59.161 
61.143 
63*029 
100.000 


5.r>i)14 

6,4720 

7.3255 

8.1622 

«M82(L 

9.7868 

10.5753 

11.3484 

12.1062 

13.5777 

14.2919 

14.9920 

15.6785 


5.2421 

6.0021 

6.7327 

7.4353 

3.1 109 

8.76(15 

9.3851 

9.9856 

10.5631 

JLIJ5JU 

1 1.6523 
12.1657 
12,6593 
13.1339 


Aum 

17.0112 14,0292 


17.6580 
A 8.2922 
18.9139 

E jfL5221. in ,- iaj 15.622L 

20,1210 15.9828 13.0032 

20.7069 16.3296 13.2105 

21.2313 16.6631 13.4062 

21.8+44 16.9837 13.5907 

MJ9/A i>7rl9?0 _! 3.16+3, 

24.9986 18.6646 14,4982 

27.3555 19.7928 15,0463 

29.4902 20,7200 15.4558 

31.4236 21,4822 15.7619 

^u7Ax^„. ia jzm<^ 

34.7609 22.6235 16.16 L4 

36.1975 23.0467 16,2891 

37,4936 23,3945 16.3845 

38.6771 23.0304 , 16.4558 

J&zm 23^15^ &502L 

40.7113 24.1035 16.5489 

41.5869 24,2673 16.5787 

42.3300 24.3978 16,6009 

43.0983 24.5050 16.6175 

50,0000 25.0000 18.182 


14,4511 

14.3568 

15.2470 


4.9173 
5.5824 
6.2098 
6.8017 
7.3601 _ 
7,8869 
8.3838 
8.8527 
9,2950 
Jj. 7 122_ 
10.1059 
10,4773 
10.8276 
11.1581 
J 1 .4699^ 
11.7641 
12.0416 
12.3034 
12,5504 


4,6229 

5.2064 

5,7466 

6.2469 

6.7101 

mwr 

7.5361 
7.9038 
8.2442 
t 8.55 95 ^ 
3.8514 
9.1216 
9,3719 
9.6036 
9 .8)81 
liwuT 

10.2007 
10.3711 
10.5288 
10.6743 ± 
10.8100 
10.9352 
11,0511 
11.1534 
11.2578 
11.6546 
11.9246 
12.1084 
12,2335 
12.3186 
12,3760 
12.4160 
12.4428 
12.4611 
12.4735 
12.4820 
12.4877 
12.4917 
114943 
12,5000 
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Uniform Series Capital Recovery Factor 
Discrete Compounding 
(A/P, r%, N) 


N 

1% 

To find A 
Given ? 
A/P ' 

2% 

To find /I. 
Given P 
A/P 

4% 

To fincM 
Given P 
A/P 

67 

To find A 
Given P 

A/P 

8% 

To find/4 
Given? 
A/P 

10% 

To Sind A 
Given P 
A/P 

12% 

To find A 
Given? 
A/P 

15% 

To find /l 
Given? 
A/P 

20% 

To find A 
Given ? 
A/P 

25% 

. 4 

TcHlnd A 
Given? 
A/P 

N 

i 

i.otoo 

1.0200 

1.0400 

1.0600 

1,0800 

1.1000 

1,1200 

1.1500 

1.2000 

2 , 1 : -.oo 

1 

2 

0.5075 

0.5150 

0.5302 

0,5454 

0.5608 

0.5762 

0.5917 

0.6151 

0,6545 

0,t *44 

2 

3 1 

0.3400 

0.34G8 

0.3603 

0.3741 

0,3880 

0.4021 

0,4163 

0,4380 

0.4747 

0,5123 

3 

4 

0.25G3 

0.2626 

0.2755 

0.2886 

0.3019 

0.3155 

0.3292 

0.3503 

0.3863 

0.42.H 

4 

jJL egn 

0.2060 

0.2122 

0.2246 

0.2374 

0.2505 

0.2638 

0.2774 

0.2983 

0.3344 

0.3718 

5 

6 

”{1.1725 

0.1785 

0,1908 

0.2034 

0.21 63 

0.2296 

0.2432 

0.2642 

0.3007 

' 0.33B8 

6 

7 

0.1486 

0.1545 

0,1666 

0.1791 

0.1921 

0,2054 

0.2191 

0.2404 

0.2774 

0.3163 

7 

8 

0.1307 

0.1365 

0.1485 

0.1610 

0.1740 

0.1874 

0.2013 

jf.2229 

0.2606 

0,3004 

8 

9 

0.1167 

0.1225 

0.1345 

0.1470 

0.160! 

0,1736 

0.1877 

0,2096 

0.24BI 

0.2838 

9 


0,1036 

0.1113 

0.1233 

0.1359 

0.1490 

0.1627 

0.1770 

0,1993 

0.2385 

0.2801 

10 

11 

UJIWo 

0.1022 

0.1141 

“ 0.1268 

0.1401 

0.1540 

0.1684 

0.1911 

0.2311 

0.2735 

11 

12 

0.0888 

0.0946 

0.1066 

0.1193 

0.1327 

0.1468 

0.1614 

0.1845 

0.2253 

0.2684 

12 

13 

0.0824 

O.0HK1 

0.1001 

0.1130 

0,1265 

0.1408 

0.1557 

0,1791 

0.2206 

0.2645 

13 

14 

0.0769 

0.0826 

0.0747 

0,1076 

0.1213 

0.1357 

0.1509 

0.1747 

0/2 1 m 

0.2615 

14 


0.0 } 2 1 

0.0778 

0,0899 

0.1030 

- 0.1168 

p.1315 

0.1468 

0.1710 

0.2139 

0.2591 

15 

16 

lW)b79 

O.U737 

0.0858 

0.0990 

0.1130 

0.1278 

0.1434 

0.1679 

0.21 14 

0.2572 

iS 

17 

0.0643 

0.0700 

0,0822 

0.0954 

0.1096 

0.1247 

0.1405 

0*1654 

0.2094 

0.2558 

!7 

18 

0.0610 

0.0667 

0.0790 

0.0924 

0.1067 

0.1219 

0.1379 

0.1632 

0.2078 

0.2546 

18 

19 

0.0581 

0.0638 

0.0761 

0.0896 

071041 

0.1195 

0.1353 

0.1613 

0.2065 

0.2537 

19 

20 

0.0554 

0.0612 

0.0736 

0.0872 

0.1019 

— 0J_LZ5 ' 

0.1339 

0.1598 

0 205+ 


20 

21 

(ITTyflJ 

" 6.0588 

0.0713 

04)850 

0.0998 ' 

0.1156 ’ 

0.1322 

6.1584^ 

0.2044 


‘ 21 

22 

0.0509 

0.0566 

0.0692 

0.0830 

0.0980 

fl.1140 

0.1308 

0.1573 

0.2037 

0.2519 

22 

23 

0.0489 

0.0547 

0.0673 

0.0813 

0.0964 

0.1126 

0,1296 

0.1563 

0.2031 

0.2515 

23 

24 

0.0471 

0.0529 

0.0656 

0.0797 

0,0950 

0,1113 

0,1285 

0.1554 

0.2025 

0.2512 

24 

25 

0.0454 

0.0512 

0.0640 

0.0782 

04937 

0.1102 

0.1275 

0.1547 

0,202 T 

0.2509 

25 

26 

0.0439 

0,0497 

0,0626 

0.0769 

9.0925 

0.1092 

0.1267 

0.J541 

o. 2 o n 

0.2508 

26 

27 

0.0424 

0.0483 

0.0612 

0.0757 

0,0914 

0.1033 

01.259 

0.1535 

0.2Q1J 

0.2506 

27 

28 

0,0411 

0.0470 

0.0600 

0,0746 

0.0905 

0.1075 

0.1252 

0.1531 

0.2012 

0.2505 

28 

29 

0.0399 

0,0458 

0.0589 

0,0736 

0.0896 

0.1067 

0.1247 

0.1527 

0.2010 

0.2504 

29 

Win i 

0.0387 

0.0446 

0.0578 

0.0726 

0,0888 

0.1061 

0.1241 

0.1523 

0.2008 

0*2503 

30 

35 

0.0340 

0.0400 

0.0536 

04)690 

. 0.0858 

0,1037 

0.1223 

-ojsTI— 

0.2003 

0*2501 

35 


0,0305 

0.0366 

0.0505 

0.0665 

0.0839 

0.1023 

0.1213 

0.1506 

0.2001 

0.2500 

40 

45 

0,0277 

0.0339 

0.0483 

0.0647 

0,0826 \ 

0,1014 

0*1207 

0.1503 

0,2001 

0*2500 

45 

so 

0.0255 

0.0318 

0,0466 

04)634 

0,0817 ; 

0,1009 

0.1204 

0.1501 

0.2000 

0*2500 i 

50 

55 

0.0237 

0.0301 

0.0452 

0.0625 

0.0812 

_ 0.1005 ri — 

, 0*120^ iiim 

0.1501 

( 0,2000 

0.2500 

55 


0.0222 

0.0288 

0,0442 

0.0619 

0,0808 N 

0.1003 

0,1201 

0.1500 

0*2000 


60 

65 

0.0210 

0.0276 

0.0434 

0,0614 

0.0805 */• 

0.1002 

0,1201 

0,1500 

0*2000 


65 

70 

0.0199 

0.0267 

0.0427 

0.0610 

0,0804 

0.1001 

0.1200 

0.1500 1 

i 

J 

70 

75 

0.0190 

0.0259 

0.0422 

0,0608 

0.0802 

0.1001 

i 0*1209 

0*1500 



75 


0.0182 

0.0252 

0.0418 

0.0606 

0,0802 

0.1000 

0.1200 

* : 0.1500 



SO 

85 


0.0246 

0,0415 

0//O04 

0.U8U1 

0.1000 

0.1200 

0.1500 



— ST 

90 

0.0169 

0.0240 

0.0412 

0.C603 

0.0B01 

0.1000 

. - .. .,i 




90 

95 

0,0164 

0.0236 

0,0410 

0.0602 

0.0001 

0*1000 





95 

100 

0.0159 

0.0232 

0.0408 

0.0602 

0.0800 

0.1000 





too 

CO 

0.0100 

0.0200 

BP>W jnrqiu 

0,0400 

O.OfWJ 

. nnnpO „ . 

LSgj- 

i i t 

, j ? * - | 

/ * 

4j : * 1 ’ 




APPENDIX D 


Gradient Series Formula Values 
Discrete Compounding 


Gradient Series To Present Worth 
Discrete Compounding 
(P/G, r%, N) 


B 

' 

2% 

47. 

67. 

8 % 

10% 

127. 

157. 

207. 

25% 

n 

1 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1 

2 

0.98 

0.96 

0.92 

0,89 

0.86 

0,83 

0,80 

0.76 

0,69 

0,64 

2 

3 

2.92 

2.85 

2,70 

2.57 

2.45 

2.33 

2.22 

2.07 

1.85 * 

1.66 

3 

4 

5.80 

5.62 

5.27 

4.95 

4.65 

4.38 

4.13 

3.79 

3,30 

2.89 

4 

5 

9.61 

9.24 

8.55 

7,93 

7.37 

6.86 

6.40 

5.78 

4.91 

4.20 

5 

6 

14.32 

1168 

12,50 

11.46 

10.52 

9.68 

8.93 

7.94 

6.58 

5.51 

6 

7 

19.92 

18.90 

17.07 

15.45 

14.02 

12.76 

11.64 

10,19 

8.26 

6.77 

7 

8 

26.38 

24.88 

22.18 

19.84 

17.81 

16.03 

14.47 

12.48 

9.88 

7*95 

8 

9 

33.69 

31.57 

27.80 

24,58 

21,81 

19.42 

17,36 

14.75 

11.43 

9.02 

9 

10 

41.84 

38.95 

33.88 

29,60 

25.98 

72.89 

20.25 

16.98 

12.89 

9*99 

10 

11 

50,80 

47.00 

40,38 

34.87 

30,27 

26.40 

23.13 

19.13 

14.23 

10.85 

11 

12 

60.57 

55,67 

47.25 

40.34 

34,63 

29.90 

25.95 

2U8 

15.47 

11,60 

12 

15 

94.48 

85.20 

69.74 

57.55 

47,89 

40.15 

33,92 

26.69 

18.5? 

13.33 

15 

20 

165.46 

144.60 

m.56 

87.23 

69.09 

55.41 

44.97 

33,58 

21.74 

14.89 

20 

25 

252.89 

214,26 

156.10 

1 15.97 

87.80 

67,70 

53.10 

38,03 

23.43 

15.56 

25 

30 

355.00 

291.72 

201,06 

142.36 

103,46 

77.08 

58.78 

40-75 

24.26 

15.83 

30 

35 

470.15 

374.88 

244.88 

165.74 

116.09 

83,99 

62.61 

42,36 

24.66 

15.94 

35 

40 

5V6.85 

461.99 

286.53 

185.96 

126.04 

88.95 

65.12 

43,28 

24.85 

15,98 

40 

45 

733.70 

551.56 

325.40 

203,11 

133.73 

92.45 

66.73 

43.81 

24.93 

15.99 

45 

50 

879.41 

642.36 

361,16 

217.46 

139.59 

94,89 

67.76 

44.10 

24.97 

16.00 

SO 

60 

1192.80 

823.70 

423.00 

239.04 

347.30 

97,70 

68.81 

44,34 

24*99 

— 

60 

70 

1528.64 

999.83 

472.48 

253.33 

151.53 

98.99 

69.21 

44.42 

— 

— 

70 

80 

1879.87 

1166.79 

511.12 

262.55 

153.80 

99.56 

69.36 

44.47 

— 

— 

80 

90 

2240.55 

1322.17 

540.77 

268,39 

154.99 

99.81 

— 

— 

— 

— 

90 

100 

2605.76 

1464,75 

563.12 

272.05 

155.61 

99.92 

— 

— 

— 

— 

100 


1 Gradient Series To Uniform Series 

Discrete Compounding 
(A/G, r%, N) 

K 

17. 

• 27c 

47o 

6% 

a% 

10% 

12% 

15% 

20% 

25% 

n 

I 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

O.OQOQ 

0.0000 

O.ODDO 

O.GOOO 

0,0000 

i 

2 

0.4974 

0.4950 

0.4902 

0.4854 

0.48O8 

0.4762 

0.4717 

0.4651 

0.4S4S 

0.4444 

2 

3 

0.9932 

0.9868 

0.9739 

0.9612 

0*9487 

0.9366 

0*9246 

0.9071 

0*8791 

0.8525 

3 

4 

1.4874 

1.4752 

1.4510 

1*4272 

1,4040 

1.3812 

1.3589 

1.3263 

1.2742 

1*2249 

4 

5 

1.9799 

1.9604 

1*9216 

1.8836 

1.8465 

1,8101 

1.7746 

1.7228 

1.6405 

1.5631 

5 

6 

2.4708 

2,4422 

2*3857 

2.3304 

2.2763 

2.2236 

2.1720 

2.0972 

1.9788 

1,8683 

6 

7 

2.9600 

2,9208 

2,8433 

2.7676 

2.6937 

2.6216 

2,5515 

2.4498 

2.2902 

2.1424 

7 

8 

3,4476 

3.3961 

3,2944 

3.1952 

3.0985 

3.0045 

2.9131 

2*7813 

2.5756 

2,3872 

s 

9 

3*9335 

3,8680 

3.7391 

3.6133 

3*4910 

3.3724 

3,2574 

3.0922 

2.8364 

2.6048 

9 

10 

4.4177 

4.3367 

4.1773 

4.0220 

3.8713 

3.7255 ‘ 

3.5847 

3.3832 

3*0739 

2.7971 

10 

n 

4*9003 

4.8021 

4.6090 

4.4213 

4.2395 

4.0641 

3.8953 

3.6549 

3.2893 

2.9663 

11 

12 

5.3313 

5.2642 

5.0343 

4.8113 

4.5957 

5.3884 

4.1897 

3.9082 

3.4841 

3*1145 

12 

15 

6.8141 

6.6309 

6.2721 

5.9260 

5.5945 

5.2789 

4.9803 

4.5650 

3,9588 

3*4530 

15 

20 

9*1692 

8.8433 

8.2091 

7.6051 

7.0369 

6.5081 

6.0202 

5*3651 

.4*4643 

3,7667 

20 

25 

11.4829 

10.9744 

9.9925 

9.0722 

8.2254 

7.4580 

6.7708 

5.BS34 

4.7352 

3*9052 

25 

30 

13.7555 

13,0251 

11.6274 

10*3422 

9.1897 

8.1762 

7.2974 

6,2066 

4*8731 

3.9628 

30 

35 

15*9869 

14.9961 

13*1198 

11.4319 

9*9611 

8.7086 

7*6577 

6,4019 

4.9406 

3.9858 

35 

40 

18.1774 

16.8885 

14.4765 

12.3590 

10,5699 

9.0962 

7,8988 

6,5168 

4*9728 

3.9947 

40 

45 

20,3271 

18.7033 

15.7047 

13*1413 

11*0447 

9.3740 

8,0572 

6.5830 

4.9877 

3.9980 

45 

50 

22.4362 

20.4420 

16,8122 

13.7964 

11*4107 

9*5704 

8.1597 

6,6205 

4.9945 

3.9993 

50 

60 

26*5331 

23.6961 

18.6972 

14.7909 

11.9015 

9,8023 

8*2664 

6*6530 

4.9989 

— 

. 60 

70 

30.4701 

26,6632 

20.1961 

15*4613 

12.1783 

9.9113 

8*3082 

6,6627 

— 

— 

70 

80 

34.2490 

29.3572 

21.3718 

15.9033 

12*3301 

9.9609 

8.3241 

6.6656 

— 



80 

90 

37.8723 

31.7929 

22,2826 

16.1891 

12.4116 

9.9831 

— 

— 

— 

— * 

90 

100 

41*3424 

33*9863 

22.9800 

16.3711 

12.4545 

9,9927 

— 

— 

— 

— 

100 
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APPENDIX E 


Interest Formula Values 
Continuous Compounding - Discrete Flows 


Single Payment Compound Amount Factor 
Continuous Compounding - Discrete Flow 
(F/P, r %, N) 


CJl 

CO 



i% 

2% 

4% 

6% 

8% 

10% 

12% 

15% 

20% 

25% 



To find F 

To find F 

To find F 

To findF 

To find F 

To find F 

To find F 

To findF 

To llndF 

To find F 

| 


Given P 

Given P 

Given P 

Given P 

Given F 

Given P 

Given P 

Given P 

Given P 

CivcnP 


N 

F/P 

F/P 

F/P 

F/P 

F/P 

F/P 

F/P 

F/P 

F/P 

F/P 

J N 

i 

1.0101 

1.0202 

1,0408 

, 1.0618 

1.0833 

1.1052 

1.1275 

1.1618 

1.2214 

1.2840 

1 1 

2 

1.0202 

1.0408 

1.0833 

! 1.1275 

1.1735 

1.2214 

1.2712 

1.3499 

1.4918 

1.6487 

2 

3 

1.0305 

1.0618 

1.1275 

1.1972 

1.2712 

1.3499 

1.4333 

1.5683 

1.8221 

2.1170 

3 

4 

1.0408 

1.0833 

1.1735 

1.2712 

1.3771 

1.4918 

1.6161 

1.8221 

2*2255 

2.7183 

4 


1*05 13 

1.1052 

1.2214 

1*3499 

l im. 

1.6482 L«221 

2.1170 

2.7183 

3.4903 


6 

^1.0618 


“T27TT” 

^^"17433^ 

1.6161 

l,822l 

2.0544 

2.454i™' 

^“^3.3201 

""““"T33T7* 

6 

7 

1.0725 

1.1503 

1.3231 

1.5220 

1.7507 

2.0138 

2.3164 

2.8577 

4.0552 

5,7546 

7 

8 

1.0833 

1.1735 

1.3771 

1.6161 

1.8965 

2.2255 

2.6117 

3.3201 

4*9530 

7.3891 

S 

9 

1.0942 

l.:972 

1,4333 

1.7160 

2.0544 

2.4596 

2.9447 

3.8574 

6*0496 

9.4877 

9 

io , , 

.... 1.1052 

1.2214 

1.4918 

-JU2SS- 


3.3201 

4.4817 

7.3891 

— 12.1825,. _ 10 

11 

1.1 163 

”"T246r” 

" “*U527"" 


3.00+2 

3.7434 TS.2W0 

9.0250 

73.^26^11 

12 

1.1275 

1.2712 

1.6161 

2.0544 

2.6117 

3.3201 

4.2207 

6.0496 

11.0232 

20*0855 

12 

13 

1.1388 

1.2969 

1,6820 

2.1815 

2.8292 

3.6693 

4.7588 

7.0287 

13.4637 

25*7903 

13 

14 

L1503 

1.3231 

1,7507 

2.3164 

3.0649 

4.0552 

5.3656 

8.1662 

16.4446 - 

33.1155 

14 

IS 

- 1.1618 

1.3499 

1.8221 

2.45% 

3,3211] 


6.0496 

9,4877 

20.0855 

42.5211 

15 

16 

1.1 7 J5 

1.3771 

I.K‘Jf.5 

2.6117 

33966^ 

“4.9530’ 

6*8210 

It. 0232 

24.5325 

54.59S2 

16 

17 

1.1853 

1.4049 

1.9739 

2.7732 

3,8962 

5*4739 

7.6906 

12.8071 

29.9641 

70.1054 

17 

18 

1.1972 

1.4333 

2.0544 

2.9447 

4.2207 

6.0496 

8.6711 

14.8797 

36.5982 

90.0171 

18 

19 

1.2092 

1.4623 

2.1383 

3.1268 

4.5722 

6.6859 

9*7767 

17.2878 

44.7012 

115.584 

19 

20 * 

1.2214 

1.4918 

2.2255 

3.3201 

125M> Tjaai 

11.0232 

20.0855 

54.5981 

148.413 

20 

21 

1,2337 


’“’TTToT*’ 

3333T" 

5.3656 

8.1662 

“””Ti42sr 

23.3 $ & 1 ^ 

66.6863 

T iyb37!&“” 

21 

22 

1.2461 

1.5527 

2.4109 

3.7434 

5.8124 

9.0250 

14.0132 

27. S 126 

81.4509 

244.692 

22 

23 

1.2586 

1.5841 

2.5093 

3.9749 

6.2965 

9.9742 

15.7998 

31.5004 

99.4843 

314.191 

23 

24 

1.2712 

1.6161 

2.6117 

4.2207 

6.8120 

11.0232 

17.8143 

36,5982 . 

121.510 

403.429 

24 

25 

1.2840 

1.6487 

2.7183 

""" 

4.4817 

- - 7.3891 

■ „ 1 17 187 , ;„| 

20.0855 

Aj _ mi 42.5211 

148,413 

518.013 

25 

26 

1.2969 

1.6820 


—TfsT 

8.0045 

13.4637 

216464" 

^ 49.4024”^ 



26 

27 

utoo 

1.7160 

2.9447 

5.0531 

.8.6711 

14.8797 

25.5337 

57.3975 

221*406 

854.059 

27 

28 

1.3231 

1.7507 

3.0649 

5.3656 

9.3933 

16.4446 

28.7892 

66.6863 

270.426 

1096.63 

28 

29 

1.3364 

1.7860 

3. 1899 

5.6973 

10*1757 

18.1741 

32.4597 

77,4785 

330.299 

1408.10 

29 

30 

1.3499 

3.8221 

3.3201 

6.0496 

11 0232 


36.5982 

90.0171 

403.429 

1808.04 


35 

1.4191 

n 2,0 "" 


””““fe'.l662“" 

16.4446 

33.1155 ' 

66*6863 

TJiTIcS' """ 

10%*63 

n '6310.(j9" M 

35 

40 

1.4918 

2.2255 

4,9530 

11.0232 

24.5325 

54.5981 

121.510 

403.429 

2980.96 

22026.5 

40 

45 

1,5683 

2.4596 

6.04% 

14.8797 

36.5982 

90.0171 

221.406 

854.059 

8103*08 

76879.9 

i 45 

50 

1,6487 

2.7183 

7.3891 

20.0855 

54.5982 

148.413 

403.429 

1808,04 

22026.5 

268337 

so 

55 __ 

1.733U — 

r r 3.0042 

9.0250 

27.1126 


J4AfiOT 

— 73S.095 — 

3827.63 

.52874,1- _ 


_5 l 

60 

1.8221 

3,3201 

J 1,0232 

36,5982 

121.510 

403.429 

1339.43 

8103.08 162755 


60 

65 

1.9155 

3.6693 

13-4637 

49.4024 

181.272 

665.142 

2440.60 

17154.2 



65 

70 

2.0138 

4.0552 

16.4446 

66.6863 

270.426 

1096.63 

4447.07 

36315.5 



70 

75 

2.1170 

4.4SI7 

20.0855 

90.0171 

403.429 

1808.04 

8103.08 

76879,$ 



75 

92 * 

... 2.2255 

4.9530 

24,5325 

121.510 


UMimIiiii 

14764.8 

162755 



SiL 


90 

95 

100 


2.4596 

2.5857 

2.7183 


6.0496 

6.6859 

7.3891 


36.5982 
44.7012 

54.5982 


221.406 

298.867 

403.429 

ws ssssL 


1339.43 

1998.20 

2980.96 


4914.77 
8103.0B 
13359.7 


90 

95 

♦ on 



} t? a tr3 C — 5 

esent Worth Factor 
ding - Discrete Flow 
r%, N) 

10 % 12 % 15 % 20 % 25 % 


To find P 
Given F 


To find P 
Given F 


To find P 
Given F 


To rind P 
Given F 


To find P 
Given F 























Uniform Series Compound Amount Factor 
Continuous Compounding - Discrete Flow 
(F/A, r%, N) 


cn 

O 


12 

13 

14 

JiL 


16 

17 

IS 

19 

JSL 


21 

22 

23 

24 

25 


90 

95 

LOO 


12.6860 

13.8135 

14.9524 

16.1026 


17.2645 

18.4380 

19.6233 

20.8205 


23.25 J 2 
24.4349 
25.7309 
26.9895 
28.2608 


145.232 

157.780 

170.971' 



1% 

2% 

4% 

6% 

8% 

10 % 

12% 

15% 

20% 

25% 



To find F 

To find F 

To find F 

To find F 

To find F 

To flnJ F 

To find F 

To Hnd F 

To find F 

To Hnd F J 

. .. 


Given A 

Given A 

Given A 

Given A 

Given A 

Given / 

Given A 

Given A 

Given A 

Given A % 


jV 

F/A 

FIA 

F/A 

F/A 

F/A 

F/A 

F/A 

F/A 

FIA 

F/A 

N 

1 

LOOOO 

LOOOO 

1.0000 

1.0000 

LOOOO 

1.0000 

1.0000 

LOOOO 

1,0000 

1.0000 

1 

2 

2*0101 

2.0202 

2.0408 

2.0618 

2.0833 

2.1052 

2.1275 

2,1618 

2.2214 

2.2840 

2 

3 

- 3.0303 

3.0610 

3.1241 

3.1893 

3.2568 

3.3266 

3.3987 

3.5117 

3.7132 

3.9327 

3 

4 

4,0607 

4.1228 

4.2516 

4.3866 

4.5280 

4,6764 

4.8321 

5.0800 

5.5353 

6.0497 

- 4 1 

J w 

— 5.101S, 

5.2061 

5-4251 5,6578 5*)i!52 6.1683 

6.4481 

6.9021 

7.7609 

■iri . 1 nj_. m 

5 1 

v 

6.1528 

oJiii 

6.6465 

7.0077 

7.3970 

7.8170 

“Ti/iir* 

TWtW 10.4792 

12.2584 


7 

7.2146 

7.4388 

7.9178 

8.4410 

9.0131 

9.6391 

10.3247 

11.4787 

13.7993 

16.7401 

7 

8 

8.2871 

8.5891 

9.24N9 

9.9630 

10,7637 

11.6528 

12.6411 

14,3364 

17.8545 

22.4947 

8 

9 

9.3704 

9.7626 

10.6180 

11.5790 

12.6602 

13.8784 

15.2528 

17.6565 

22.8075 

29.8837 

9 

10 

10.4646 

10.9598 

12.0513 

13.2950 




18.1974 

21.5139 

28,8572 

— M 

10 


13.4273 

14.6985 

15.9955 

17.3186 


18.6685 

20.0456 

21.4505 

22.8839 

25.8380 

27.3599 

28.9126 

30.4967 

32.1128 


15.0959 
16.7119 
18.3940 
■ 2(U4 46_, 


21.9668 
23.8632 
25.8371 
27.8916 
0291 
32.2554 
34.5717 
36.9826 
39.4919 


16.9402 

17.0519 19.3511 22.0604 

19.1064 21.962B 25.3806 

21.287B 24.7920 29.0499 

23.6042 77.HS6Q 31.1051 

26.0638 31.1770 37- 

28.6755 34.7736 42.5398 

31.4487 3 8.6698 48.0137 

34.3934 42.8905 54.0634 


40.8403 

44.3657 

48.1091 

52.0840 

56.3047 


52.4158 

57.7813 

63.5938 

69.8903 


68.1383 

76.3045 

85.3295 

95.3037 


25.2610 

29.4817 

34.2405 

39.6061 

?S.65?r 

52.4767 

60.1673 

68.8384 

78 .6151 

89.6iir 

102.067 

116.080 

131.880 

149,694 


31.2026 

37.2522 

44.2809 

52.4471 

Vi.ttW 

72.9580 

85.7651 

100.645 

117,933 

-MoIT" 

161.354 

18B.467 

219.967 

256.565 


45.2712 

56.2944 

69.7581 


. 67.1966 
87.2821 - 
113.073 


86.2028 146.188 

106.288 188.709 

130.821 243.307 

160.785 313.413 

197.383 403.430 


296.682 

363.369 

444.820 

544.304 


667.427 

857.993 

1102.69 

1416.88 


249.966 

281.461 

316.269 


872.275 

1070.82 

1313.33 




3564.34 

4817.01 

6507.94 


16070.1 

23979.7 

35779.3 


77037.3 

127019 

209425 

a? 


12 

13 

14 

15 


tr 

!7 

18 

19 

20 


21 

22 

23 

24 


26 

Z9.5448 

33.70 15 

44.3219 



HDEES3 


tmsmsa 



27 

30.8417 

35.4435 

47.6511 

65.5452 



192.426 

348.489 

995.500 

3003.46 

27 

28 

32.1517 

37.1595 

50.5958 



146.853 

217.960 

405.886 

1216.91 

3857.52 

28 

29 

33.4748 


53.6607 

75.9639 

110.169 

163.298 

246.749 

472.573 

1487.33 

4954.15 

29 

■■■ 

34.3113 

40.6962 

56.33(16 

81.6612 

_ 120,345 

181.472 

. 279.209 

550.051 

1817.63 

6362.26 

30 

35 

41.6976 


74.8626 

1 15.889 

185,439 

305.364 



MMi 

22215.2 

ST" 

40 

48.9370 

60.6663 

96,8625 

162.092 

282.547 

509.629 

945.203 

24B6.67 

13459.4 

77547.5 

50 

45 

56,5476 

72.2528 

123,733 

224.453 

427.416 

846.404 

1728.72 

527J.19 

36594.3 

2)0676 

45 

so 

64,5483 

85.0578 

156.553 

308.645 

643.535 

1401.65 

3156.38 

11166.0 

99481.4 

944762 

50 

■■$91 

— miuj 





23645.3 


. , 


60 

81.3015 

114.350 

245.601 

575.683 

1446,93 3B26.43 10497,8 

50064.1 

735103 


60 

65 

91.0971 

132.135 

305.403 

781748 

2164.47 

6314.88 

19134.6 

105993 



65 

70 

100.869 

151.238 

378,445 

1062.26 

3234.91 

10417.6 

34872.0 

224393 



70 

75 

111.143 

172.349 

467.659 

1439.56 

4831.83 

17182.0 

63547.3 

475047 ■ 



75 

-SO. 

121.942 



_ 57.6.6 25^,^ 

1948.85 

221415 28314.4 

11S797 

1005680 



- .... 80 

85 

133.296 

22L468 

709.717 

2636.34 

10768”! 

46721.7 

, — ■ 

- fc — 



85 


91 

95 

100 





Uniform Series Present Worth Factor 
Continuous Compounding - Discrete Flow 
(P/A, r%, N) 


To find P 
Given A 
P/A 

0.9802 
1*9410 
2,8828 
3.8059 
4.7107 
5. 5976 
6.4070 
7.3191 
8.1544 
_8.973L 
9.7750 
10.5023 
11,3333 
12,0891 
. 12,8299 
13,5561 
14.2678 
14.9655 
15.6494 
_ 16.3197 
16.9708 
17.6208 
18.2521 
18.8709 
T ,19.4774 
20.0719 
20,6547 
21.2259 
21.7853 
22.3346 
24.51 99 
27.2591 
29.3758 
31.2910 

34.5921 

36.0109 

37.2947 

38.4564 

* 22 * 5075 ^ 

40.4585 

41.3191 

42,0978 

42.8023 


To HndP 
Given >1 
P/A' 

0.9608 

1.8839 

2.7708 

3,6230 

4.4417 

"TSSF" 

5.9841 
6.7103 
7.4079 
,,,8.0783 
8.7223 
9,3411 
9.9356 
10.5068 
1 1 ,0 556 
TL5SIT 
12.0895 
12.5763 
13.0439 
13.4933 
13.9250 
14.3398 
14.7383 
15.1212 
t 15.4891 
15.8425 
16.1821 
16,5084 
16.8219 
17,1231^ 
18.4609 
19.5562 
20.4530 
21.1872 

Jim. 

22.2804 
22.6834 
23.0133 
23.2834 
_ 23.5045^ 
23.6856 
23.8338 
23.9552 
24.0545 


To find P 
Given A 
P/A 

0.9418 
1.8287 
2.6640 
3.4506 
m J. 1914 |ui 
4.B891 
5.5461 
6.1649 
6.7477 
_ 7.2965 
7.8133 
8.3001 
8.7585 
9.1902 
,,9.5968^ 
"^9797“ 

10.3402 
10.6798 
10.9997 
11.3009 
11.5845 
11.8516 
12.1032 

12.3402 
12,5633 
12.7734 
12.9713 
13.1577 
13.3332 

^13.4985 
14.1913 
14.7046 
15.0848 
15.3665 
15.5752 | L 
1 5.7298 
15.8443 
15.9292 
15.9920 
16.03B6 
T l6.07ir 
16.0986 
16.1176 
16.1316 


To HndP 
Given A 
P/A 

0.9231 
1*7753 
2.5619 
12880 
3.9584 
^4.5771 
11483 
5.6756 
6.1624 
„6.6U7 
7.0265 
7.4094 
7.7629 
8.0B91 
_ 8,3903 
8.6684 
8*9250 
9.1620 
9.3807 

9.7689 
9,9410 
10.0998 
10.2464 
,10.3817 a 
10.5067 
10,6220 
10,7285 
10.8267 
,10.9174 
11.2765 
11,5172 
11.6786 
U.786B 

jl im~ 

11.9079 

11,9404 

11.9623 

11.9769 

11.9867 

11.9933 

11.9977 

12,0007 

12.0026 


To fmd P 
Given A 
P/A 

0.9048 

1.7236 

2.4644 

3.1347 

4,2900 

4*7866 

5.2360 

5.6425 

JkOi&U 

6.3433 

6.6445 

6.9170 

7.1636 

-7.38 67^ 

7.5886 

7.7713 

7.9366 

8,0862 

-82215^ 

8.3440 

8.4548 

8.5550 

8.6458 

^127^ 

8.8021 

8,8693 

8.9301 

8.9852 

_9,0349 

9.2212 

9.3342 

9.4027 

9.4443 

^4695 

9.4848 

9.4940 

9.4997 

9.5031 

jyosi 

*~9.5064 ,1L 

9.5072 

9.5076 

9.5079 


To findP 
Given A 
P/A 

c? * " *ua 

0.8869 

1.6735 

2.3712 

2.9900 

3.5388 

ToES” 

4.4573 

4,8402 

5.1798 

5.4810 

T.ftar 11 

5.9850 

6.1952 

6.3815 

6,5468 

-rzmr* 

6.8235 

6.9388 

7.0411 

7,13 18 

TTnr 

7,2836 

7.3469 

7,4030 

7.4528 

7.5362 

7,5709 

7.6017 

7.6290 

"Tv35 7" 

7.7788 

7.8079 

7,8239 

7,8327 

7^8375 ' 

7.8401 

7.8416 

7*8424 

7.8428 


ToflndP 
Given A 

PIA 

SSSfSSSSSS 

0.8607 

1.6015 

2.2302 

2.7880 

3.2603 

TOT* 

4.0168 

4.3180 

4.5773 

4.sn 04 

"TWIT 

5.1578 

5.3000 

5.4225 

5.5279 

teist 

5.6967 
5.7639 
5.8217 
, 5.8715 
5.9144 
5.9513 
5.9830 
6.0103 
6.0338 

"To! TT 

6.0715 

6.0865 

6.0994 

6.1105 

■w 

6.1638 

6.1719 

6.1757 

3.1775 

6.178+' 

6.1788 

64790 

64791 
64791 


To find P 
Given A 
P/A 

H- ii. 

0.8187 

1,4891 

2.0379 

2,4872 

2.8551 

3.1563 

3.4029 

3*6048 

3.7701 

3.9054 

"To 162” 

4,1069 

4*1812 

4.2420 

4.2918 

4.3325"“ 

4.3659 

4*3932 

4.4156 

4,4489 
4.4612 
4.4713 
4*4795 
_ 4,4862 
4*4917 
4.4963 
4,5000 
4.5030 
_ 4.5055 
4.5125 
4.5151 
4/5161 

4.5165 
_4,5166 

4.5166 


To find P 
Given A 
P/A 

0.7788 

1,3853 

1*8577 

2.2256 

JU&21-T 

2.7352 

2.9090 

3.0443 

3.1497 

3.2957 
3. 3455 
3*3843 
3.4145 

-iiaaa- 

3.4563 

3.4706 

3.4817 

3.4904 

J.i2ZJU 

3.5023 

3.5064 

3.5096 

3.5121 

JJLW5L- 

3.5155 

3.5167 

3.5176 

3.5183 

M. 

3.5203 

3.5207 

3.5203 

3.5208 


APPENDIX F 


Interest Factor Values 

Continuous Compounding - Continuous Uniform Flow 



1^4 tel ^.J M~3 tkrd c-rrrr3 fe^S M gM 


Uniform Series Compound Amount Factor 
Continuous Compounding - Continuous Uniform Flow 
(F/A, r%, N) 



i% 

2% 

• 

4% 

6% 

8% 

10 % 

12% 

15% 

20% 

25% 



* To flndF 

To fUidJ? 

TottndF 

To flndi^ 

To find F 

To And F 

To find F 

To And F 

To find F 

To find F 



Given A 

Given A 

Given A 

Given 

Given A 

Given A 

1 Given 3 

Given A 

Given % 

Given A 



FfA 

FfA 

F/A 

FfA 

F/A ^ 

F/A 

F(A 

FfA 

f/a 

F/A 

N 

1 

1.0050 

i.oioi 

1.0203 

1.0306 

1.0411 

1.0517 

1*0625 

1.0789 

1.1070 

1.1361 

1 

2 

2.0201 

2.040S 

2.0822 

2,1249 

2,1689 

2.2140 

2.2604 

2.332+ 

2.4591 

2,5949 

2 

3 

3.0455 

3.0918 

3.1874 

3.2870 

3.3906 

3.4986 

3.6111 

3.7887 

4.1106 

4.4680 

3 

4 

4.0811 

4,1644 

4.3378 

' 4.5208 

* 4.7141 

4*9182 

5.1340 

5.4808 

6.1277 

6.8731 

4 

5 

5.1271 

£.2585 

5.5351 _ 

5,8310 

6.1478 

■^n„ m ,6-4872 w 

6.8510 

7.4467 

8.5914 

9.9614 

5 

T" 

~ £1837 

'"0748™*" 

6.7812 

7.2222" 

7*7009 

8.2212 

8.7869'” 

*” *).7j07” 

1 1.6006 

13.926B 

J 6 

T 

7.2508 

7.5137 

8,0782 

8.6994 

9.3834 

10.1375 

10.9697 

12.3843 

15.2760 

19.0184 

7 

B 

8.3287 

8,6755 

9.4282 

10.2679 

11.2060 

12.255+ 

13.4308 

15*4674 

19.7652 

25.5562 

8 

9 

9.4174 

9.8609 

10.8332 

11.9334 

13*1804 

14.5960 

16.2057 

19.049$ 

25.2482 

33.9509 

9 

10 

10.5171 

11.0701 

12.2956 _ 

13,702!) — 


17.1828 19.3343 

23.2113 

31.9453 

44*7300 

10 

11 

1 1.6278 

' "Uioja 

13,8177 

15.5799 

17.6362 

20.0417 

22.86)8 

^”*20465" 

40.1251 

58*5705 


12 

12.7497 

13.562S 

15.4019 

17.5739 

20.1462 

23.2012 

26.8391 

33.6643 

50*1159 

76.3421 

12 

13 

13.8828 

14.8465 

17.0507 

19.6912 

22.8652 

26.6930 

31.3235 

40.1913 

62.3187 

99.161+ 

13 

14 

15.0274 

16.1565 

18.7668 

21.9394 

25.8107 

30.5520 

36.3796 

47.7745 

77.2232 

128.462 

14 

15 

16.1834 

17.4929 

- 20.5530 

24,3267 — 

-29.0015 _ 

,34.fll6^ rn . 

42.0804 

56,5849 

95^277 

166.084 

15 

16 

17,3511 

*”TOs54 

22.4120 

26.8616 

32.4580 

39.5303 

48.5030 

66.82)2™" 

*‘ i "ll7l(.63'" 

214.393 

16 

17 

18.5305 

20.2474 

24.3469 

29.5532 

36.2024 

44.7395 

55.7551 

78.7140 

144.820 

276,422 

17 

18 

19.7217 

‘ 21.6665 

26.3608 

32.4113 

40.2587 

50.4965 

63.9261 

92.5315 

177*991 

356.068 

18 

19 

20.9250 

23.1142 

28.4569 . 

35.4461 

44.6528 

56.8589 

73.1390 

108.5S5 

218.506 

458.337 

19 

20 

22.1403 

24.5912 

30.6385 

38.6686 

, 49 . 4 j. 2 s_ 

63.8906. 

83.5265 

127.237 

267.901 

589,653 

20 

*21“ 

ziAtfl 

^TTJ5ST^ 

32.9092 

42,0904 

54.5694 

71.6617 

tail nr " 

148.967 

328*432 

750,265 

21 

22 

24.6077 

27.6354 

35.2725 

45.7237 

60.1555 

80.2501 

10B.443 

174.084 

402.254 

974*768 

22 

23 

25.8600 

29,2037 

37.7323 

49.5817 

66.2067 

89.7418 

123.332 

203.336 

492*422 

1252.76 

23 

24 

27.1249 

30.8037 

40.2924 

53,6783 

72.7620 

100.232 

140.119 

237.322 

602552 

1609.72 

24 

iL- 

28.4025 

32.4361 

- 42 9570 


1 7Q.K632 

— H1825 

159.046 

. 276.807 

737*066 

2068.05 

25 

26 

29.6930 

34.1014 

45.7304 

62.6470 

87.5559 

124.637 

180,386 

322,683 

9111.561 

^ 2656.57 

26 

27 

30.9964 

. 35.8003 

48.6170 

67,5515 

95.8892 

138.797 

204.448 

375.983 

1102,03 

3412.23 

27 

28 

32.3130 

37.5336 

51.6214 

72.7593 

104.917 

154.446 

231.577 

437.909 

1347.13 

4382*53 

28 

29 

33.6427 

39.3019 

54.7483 ' 

* 78.2891 

114.696 

171.741 

262.164 

509.856 

1646.50 

5628*42 

29 

30 

34.9859 

41.1059 

S8.O029 

R4.1608 

125.290 

mmur 1 90.8S5 ^ 

296.652 

593,448 

2012.14 

7228.17 

30 

35 

41.9068 

50.6874 ' 

• 7O.3SU0 

119,436 

193.058 

321.154 

" ' S^.5it6 

""05373“*" 

5478. l r 7 

iaiibji 

111 j -sr 

40 

49.1825 

61.2770 

98.8258 

167.053 

294,157 

535.982 

1004.25 

2682.86 

14899.8 

88101.9 

40 

45 

56.8312 

72.9802 

126.241 

231,329 

444.978 

890.171 

1836.72 

5687.06 

40510.4 

307516 

45 

50 

‘ 64,8721 

85.9141 

159.726 

318.092 

669.977 

1474.13 

3353.57 

12046.9 

110127 

1073350 

50 

55 

73,3253 

100.208 

200.625 

435.211 

1005.64 

2436.92 

6117.46 

2S510.8 

-299366 


55 

60 

82.2119 

116.006 

250.579 

593.304 

1506.38 

4024.29 


1 54013,9™” 

813769 


^ " 60 

65 

91.5541 

133,465 

311.593 

806.708 

2253,40 

6641.42 

20330,0 

114355 



65 

70 

101,375 

152,760 

336.116 

1094,772 

3367,83 

10956.3 

37050*6 

242097 



70 

75 

il 1.700 

174,084 

477.138 

1483.619 

5030.36 

18070*7 

67517.4 

512526 

. ' 


75 

80 

122.554 

197,652 

588.313 __ 

2008.507 

75 ^ 5f> 


123032 

1085030 



80 

85* 

133.965 

223.697 

724.102 

2717.032 

11210.6 

49137.7 


— * - - 



as" 

90 

145.960 

252.482 

889.956 

3673*440 

16730.4 

81020.8 





90 

95 

158.571 

284.295 

1092.530 

4964.457 

24964.9 

133537 





95 

100 

171.828 

319.453 

1339.954 

6707,146 

37249.5 

220255 





100 


ifi 


1 

( 


f 

. i, '■ 


r 

?! 
*i • 

■ i 

'i 
; .1 

ij - 

\\ 

i s 

if 

If 








Uniform Series Present Worth Factor 
Continuous Compound ing_- Continuous Uniform Flow 

(P/A, r%, N) 



1% 

2% 

4% 

6% 

8% 

10% 

12% 

15% 

20% 

25% 



To find P 

To find P 

To find P 

To find F* 

ToTmdP 

To find P 

To find P 

To rind P 

To find P 

To find P 



Given A 

Given A 

Given A 

Given 3 

Given A 

Given A 

Given 

Given A 

Given?! 

Givenvi 


jV 

P/A 

PlA 

PlA 

PlA 

PlA 

P/A 

P/A 

P/A 

pit : 

P/A 

jV 

1 

0.9950 

0.9901 

0.9803 

0.9706 

0.9610 

0*9516 

0.9423 

0*9286 

0.9063 

0.8848 

i 

2 

1.9801 

1.9605 

1.9221 

1.8847 

1.8482 

1.8127 

1.7781 

1*7279 

1.6484 

1.5739 

2 

3 

2.9554 

2.9118 

2.8270 

2.7455 

2*6672 

2,5918 

2.5194 

2.4158 

2.2559 

2*1105 

3 

4 

3.9211 

3.8442 

3.6964 

3.5562 

3,4231 

3.2968 

3*1768 

3.0079 

2.753+ 

2.5285 

4 

5 

4.8771. _ 

4.7581 

4*5317 

4*3197 

4.1210 

3.9347 

3.7599 

3,5176 

3.1606 

2.8540 

5 

6 

5.8235 

5*6540 

5*3343 

5.0387 

4*7652 

4,5119 

4.2771 

3.9562 

3*4940 

3.1075 “ 

6 

7 

6.7606 

6.5321 

6.1054 

5.7159 

5.3599 

5.0341 

4.7357 

4.3337 

3.7670 

3.3049 

7 

8 

7.6884 

7.3928 

6.i’463 

6*3536 * 

5*9088 

5*5067 

5.1426 

4.6587 

3.9905 

3.4587 

8 

9 

8.6069 

8.2365 

7*5581 

6*9542 

6*4156 

5.9343 

5*5034 

4.9384 

4.1735 

3.5784 

9 

10 

9.5163 

9.0635 

8.2420 

7.5198 

6.8834 

6.3212 

5.8234 

5.179! 

4.3233 

3.6717 

10 

"ii 

10.4166 

9.8741 

8.8991 

ii.0^5 ~ 

— rcrcr— 

6.6713 

'"'6.1073 “ 

" r 5.38 6 J 

4.4460 7 

3.7443 

11 

12 

11.3080 

10.6686 

9.5304 

8.5541 

7.7138 

6.9881 

6.3589 

5.5647 

4.5464 

3.8009 

12 

13 

17.1905 

11.4474 

10.1370 

9.0266 

B.0B18 

7.2747 

6.5822 

5.7182 

4.6286 

3.8449 

13 

14 

13.0642 

12.2108 

10.7198 

9.4715 

8.4215 . 

7.5340 

6.7802 

5.8503 

4.6959 

3.8792 

14 

15 _ 

119292 

12*9591 

11.2797 

9.8905 

8.7351 

7.7687 

6*9558 

5,9640 

4*751! 

3.9059 

15 

* 16 

14.7856 

13*6925 

11.3177 

10.2851 

' ' "^0245 ”” 

7.9810 

tJ * r 7.lil6 

6.0619 r' 

4.7962 

3.9267 

16 

17 

15.6335 

14*4115 

12,3346 

10.6568 

9.2917 

8*1732 

7.2498 

6.1461 

4.8331 

3.9429 

17 

18 

16.4730 

15.1162 

12.8312 

11.0067 

9.5384 

8.3470 

7.3723 

6.2186 

4*8634 

3.9556 

18 

19 

17.3041 

15.8069 

13.3083 

11*3363 

9.7661 

8.5043 

7.4810 

6*2310 

4.8881 

3*9654 

19 

20 

18.1269 

16.4840 

13.7668 

11.6468 

9*9763 

W El! M Mil 

7.5774 

6.3 34S 

■■nm 


m 

21 

1 

17.1477 

14.2072 

11.9391 

10.1703 • " 

8*7754 

7.6628 




» 21 

'22 

' 22 

19.7481 

17.7982 

14.6304 

12.2144 



7.7387 


4.9386 

3*9837 

23 


18.4358 


12*4737 


. 8.9974 

7.8059 

S'?’ 4 

4.9497 

3.9873 

23 

24 

21*3372 


15.4277 

12*7179 

10.6674 


7.8655 

6*4845 

4.9589 


24 

25 

22.1199 


15.8030 

12.9478 

10,8083 

9.1791 

7.9184 


4.9663 

3.9923 

25 

' 26- 

22.8948 

20.2740 

16,1636 

" 13.1644 ' 


9,2573 

" 7,9654 "™" 

6*5317 

4.9724 


26 

27 

23.6621 

20*8626 

16.5101 

13.3684 

11*0584 

9.3279 



4*9774 

3.9953 

27 

2S 

24*4216 

21.4395 

16*8430 


11*1693 

‘9.3919 


6.5667 

4.9815 

3.9964 

28 

29 

25.1736 


17.1628 

13.7413 

11.2716 

9.4498 

8.0766 

6.5806 

4.9849 

3.9972 

29 


25.9182 

22.5594 

MlrAU'JIM 

13,9117 

11.3660 

9.5021 

8.1056 

6.5926 

4.9876 

3.997^ 

HI 

35 

29*5312 

25.1707 

18.8351 

14.6257 

"" 11.7399 ™ " 


8.20B4 

6.6317 

4.9954 

3.9994 

35 

40 

32.9GH0 

27.5336 

19.9526 

15.1547 

11*9905 

9.81 6B 



- 4*9983 

3.9998 


45 

36.2372 

29*6715 

20*8675 

15.5466 

12.1585 

9.8389 

8.2957 

6.6589 

4.9994 

3.9999 

45 

50 

39.3469 

31.6060 

21.6166 

15.8369 

12.2711 

9.9326 

8.3127 


4.9998 



55 

mfWrnM 1 

33.3564 

22.2299 

16.0519 

12*3465 

9*9591 

8.3220 

6.6649 

4.9999 

! .. . 

j_55 

60 

45.1188 

34.9403 

22.7321 




■‘fmwrni 


5.0000 


UKU 

65 

47.7954 

36.3734 

23.1432 

16.3293 

12,4310 


8.3299 

6.6663 



65 

70 

50*3415 

37.6702 

23.4797 

16.4167 

12.4538 


8*3315 

6*6665 



70 

75 

52*7633 

38.8435 

23.7553 

16.4815 

12.4690 

9*9945 
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Shaping Mustang I Life Cycle 
By An Interplay Of Product Innovation Rates 
Versus Forces Of Product Mortality 


65 




SUMMARY 


The Mustang I life cycle is portrayed as a dynamic resultant of an 
interplay of two countervailing forces: the force of an innovating per- 

formance to keep the product alive against a force of mortality generated 
by the introduction of internal and external competing effects to re- 
place the product by other similar or equivalent products. A generalized 
two state Markov process is used to caricature this dynamic interplay of 
innovative forces against forces of mortality. The actual Mustang 1 
product life cycle is first smoothed out to dampen random, seasonal and 
cyclic effects. Thereafter, a "best fit" is obtained to match the Markov 
model results against the smoothed out empirical results. This matching 
process synthesizes a force of innovation and a force of mortality, which 
are then interpreted in terms of pertinent observable phenomena. This 
approach can provide a method of accounting for performance of an 'nter- 
prise in a specific environment and in relation to a specific product 
beyond conventional measures of business success. A discussion is pro- 
vided on the applicability of the concept for marketing and product 
planning needs of an enterprise. 


INTRODUCTION 


Life cycles of products, enterprises, societies and cultures are 
considered as a more or less established empirical fact. The use of the 
life cycle concept in new product and marketing planning is a far less 
obvious fact. We may divorce ourselves from the heroic attempts to 
interpret the rise and fall of Roman Empire; nevertheless, we remain 
immensely curious about the life cycles of very familiar entities over 
our rather short life spans of experience. This is particularly true 
in the case of practical business planning of new products in a setting 
of competition and uncertainty. The rather well recognised life cycle 
concept appearing through historical, political, social and business 
literature closely parallels that of the biological life cycle processes 
of living organisms: there is a birth, growth, maturation, decline and 

eventual termination of the process. Actuarial practices utilize this 
concept in insurance business. In the practice of new products planning 
and marketing the practical and operational usefulness of the life cycle 
concept is not so obvious before the fact. Usually, the life cycle con- 
cept appears as a post mortem judgment after the fact. An important 
aspect of planning is prediction and forecasting of the faith of a 
product in the future market settings. A life cycle of a product should 
be understood in such relevant terms that are subject to some meaningful 
forecasting processes. The main purpose of this article is to investi- 
gate this possibility in the context of a specific example. 

The literature on product life cycle concepts is considerable. A 
small sample of references are provided here [360-372]. A general impres- 
sion obtained from this literature is that the product life cycle 
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concept is not particularly useful in planning a new product and pre- 
dicting its faith in a marketplace. The purpose of this article is not 
to review in detail this literature, but merely point out that there is 
an agreement that a process of innovation is needed to support a product 
against a variety of forces tending to kill this product in the market- 
place . 

It is possible to account for innovative and renovative performance 
of an enterprise in terms of the life cycles of its product.-, one by one 
or in an appropriate composition. Such an accounting could be done in 
financial terms, in terms of real output, or in terms of some relevant 
measure of the "worth" of the activity generated by this enterprise. 

With a dimensional consistency of a measure of performance, the same 
could be done with the effects that promote the mortality of this enter- 
prise within and without in a setting of a particular environment. At 
least in some relative terms one can introduce a measure for a force of 
renovation and innovation on one h'tid, and a force of mortality on the 
other hand. Then is is possible to develop a simple model in which 
these two opposing forces produce in a dvnauic interplay a bounded prod- 
uct life cycle [370], With a meaningful interpretation such a life cycle 
model can be applied in post mortem to interpret such phenomena as the 
railroad passenger service life cycle in the United States [371], 


A TWO-STATE LIFE CYCLE MARKOV MODEL 


References [370,371] provide a derivation of a simple generalized two- 
state Markov process life cycle model repeated here. The two states for 
the model are as follows: 

State #0: A product is produced by the enterprise and is desired by 

the society or the marketplace. 

State #1: The product cannot be produced by the enterprise or it is 

not desired by the society or the marketplace. 

Let A(t) be the probability that the product is "alive", i.e. in the 
state #0, at the time t. Let D(t) be the probability that the product 
is dead, i.e. in the state #1, at the time t. Further, let h(t) be the 
hazard or force of product mortality (events of death per unit time) 
pushing the life of the product from the state #0 toward the state #1. 

Let v(t) be the innovative recovery rate or force of innovation or force 
of renovation forcing the product from the state #1 toward the state #0. 
Then 

(1) dA(t)/dt = -h(t)A(t) + v(t)D(t) 

A(t) + D(t) = 1 
A(0) = A 0 ; 0 < A q < 1 

As shown in references [370,371], the appropriate solution to this differ- 
ential equation with the indicated constraints is as follows: 
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t 


(2) 

A(t) = S(t)U(t)[ / [v(x)/S(x)U(x)]dx + A q ] 


0 

where 

t 

(3) 

SCt) = exp[- / h(x) dx] 


0 

and 

t 

(4) 

U(t) - exp[- / vCx) dx] 


0 


S(t) is the survival probability of the product without any supportive 
innovation, and U(t) is the probability that the innovation will not 
occur within time t; S(0) = U(0) = 1. Therefore, knowing h(t) and v(t) 
will suffice to shape the probability A(t) that the product is alive. 

Let K be an appropriate dimensional scaling constant. Then K A(t) can 
be treated as the measure of the output for a particular product by the 
enterprise. In the spirit of references [370,371] Figure G-l illustrates 
the flow-diagram of the two-state Markov process life cycle generator. 
h(t) and v(t) represent a pair of relative forces which with A Q and K 
suffice to provide a match to a bounded life cycle of a "reasonably well 
behaved type". The idea is now to match an empirically observed life 
cycle in an appropriate manner and then extract out h(t) and v(t) for 
further interpretation. If h(t) and v(t) can be explained in some rel- 
evant manner subject to possible forecasting, then it would be possible 
to predict also the shape of an expected product life cycle. Such a 
possibility would be operationally helpful to new products and market- 
ing planning processes. A particular product life cycle history could 
be helpful in testing out such a possibility. In a previous study the 
product life cycle of the US passenger railroad service was investigated 
[371], The results had some encouraging correlation with other empirical 
investigations [372]. This preliminary study motivated the investigation 
described subsequently. 


THE MUSTANG I PRODUCT LIFE CYCLE 


Specific lines of automobile products provide good examples of 
product life cycles. The Mustang I product life cycle is a particular 
example. It could have been just as well a life cycle of some other 
specific brand and line of automobiles. Figure G-l illustrates the flow 
diagram of the two state Markov process discussed previously [370,371]. 
The quarter!; data for Mustang I actual production and the respective 
smoothed data by a four quarter moving average process from the year 
1964 to the year 1972 is given in Table 1. Figures G-2 and G-3 illus- 
trate this data. The purpose of this investigation was to match the 
life cycle generated by the Markov Life Cycle Model to the smoothed data 
for Mustang I, and thereby extract out the respective force of innovation 
v(t) and the force of mortality h(t) that generated this seemingly "best 
fit". This is a trial and error process indicated by the flow diagram 
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TABLE 1 


Calendar 

Year 

Quarter 

Actual 

Production 

4 Quarter 
Movi ng 
Averaqe 

Calendar 

Year 

Quarter 

Actual 

Production 

4 Quarter 
Moving 
Average 

1964 

I . 

6,954 


1968 

I 

119,671 

97,536 


II 

77,315 

i 


II. 

78,084 

87,719 


III 

112,120 



Ill 

48,556 

81 ,918 


IV 

107,019 

75,852 


IV 

98,906 

86,554 

1965 

I 

141,702 

109,539 

1969 

I 

50,323 

68,968 


II 

168,796 

132,409 


II 

76,064 

68,462 


III 

115,034 

133,138 


III 

43,420 

67,178 


IV 

154,655 ’ 

145,047 


IV 

73,384 

60,798 

1966 

I 

170,592 

152,269 

1970 

I 

41,256 

58,531 


II 

160,638 

150,255 


II 

43,135 

50,299 


III 

108,114 

148,499 


III 

27,750 

46,631 


IV 

141,423 

145,192 


IV 

53,184 

41 ,331 

1967 

I 

124,009 

133,546 

1971 

I 

31,983 

39,013 


II 

117,354 

122,722 


II 

31,577 ‘ 

36,123 


III 

71,758 

113,636 


III 

34,028 

37,793 


IV 

81,361 

98,620 


IV 

33,104 

32,673 





1972 

I 

27,968 

31,674 


Sources: 1. Automotive Industries, March 15, 1965-1970 Issues. 

2. Economics, Research and Statistics Dept., Automobile Mfrs. Asso., Inc., Detroit, Michigan 48202. 


of Figure 6-1. Figures G-4 and G-5 illustrate ■Hie result of such a 
trial and error search process. With some insight to this matching 
process the number of trials is not very large. Corresponding to this 
seemingly "best fit" life cycle modeling, the following forces were 
found: 


The force of innovation v(t) - 1.1 innovation events/year 
The force of mortality h(t) = 0.25 failure events/year 

It is readily noted that the force of innovation or renovation is 
essentially constant over time while the force of product mortality in- 
creases proportional to the 3/2 power of time from the starting moment 
of the product life cycle at the time t = 0 (at about third quarters of 
the year 1964 for the smoothed data). 

The interpretation of the force of innovation or renovation seems 
not so hard: the automobile business by tradition introduces new models 

on an annual basis on the top of which there are soma technological 
i Approvements (such as learning curve effects). Therefore, each year the 
potential automobile buyer is exposed to a "new" model . In addition to 
the new model, there are also new optional features including such items 
as window defoggers, sun roofs, "Fire-Frost" paints, ride and handling 
packages, and so on. It is reasonable to assume that the consumers' 
perception of a "new model" represents a nearly constant innovation 
event per year. Such an innovative event is generated by total business 
planning. This innovation per year can be related to new options, minor 
technical improvements, new consumer product perceptions created by 
advertisement, improvements in sales promotion, and so on. It should 
be noted that the concept and perception of a "New Model" is an aspect 
of American consumer expectations . 

The force of mortality extracted out by fitting the model into the 
smoothed data, h(t), is not so obvious to interpret at the outset. 

There can be a multitude of factors affecting product survival in the 
marketplace, such as: 

Decline in the number of buyers in a specific age group, e.g. 

18-35 year old males, or females, etc. 

Decline of number of buyers in a specific income bracket, e.g. 

$10,000 - $15,000 per year, etc. 

Introduction of new "product substitutions" by competitors or 

by the firm producing the product subject to forces of mortality. 

Increase in vehicle sales prices. 

Increase and change in vehicle size, weight, characteristics, 

etc. 

Decrease in vehicle operating cost efficiencies. 
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Many of these specific forces affecting the product life cycle might be 
aggregated into some relevant single measure such as the total competi- 
tors' market share of effective "product substitutes" that could affect 
a specific Darticular product perception in the minds of customers 
associated with, for example. Mustang I. 

In order to compare Mustang I to "Mustang-like" automobiles certain 
assumptions must be made. First, the number of people interested in this 
type of automobile is assumed to be slowly growing and with some relative- 
ly fixed characteristics such as a certain age bracket (say, 18-35 years) 
and income. With a slowly growing total market for the "Mustang-like" 
cars the Mustang I sales deteriorated while that for all the substitute 
vehicles increased. This would suggest that the force of mortality for 
Mustang I is somehow related to the total market share of the competition. 
Table 2 lists the automobiles that the authors of this article considered 
reasonably "Mustang-like". The table gives the four quarter moving 
average of the actual production divided by 1,000 for these various 
"substitute-products" of Mustang I. Since there was some question whether 
VW and Toyota would really be sufficiently "Mus^ang-like" , totals were 
obtained with and without these two types of cart' . The final results 
did not seam to be affected strongly either way. 

Figure G-6 illustrates the total market share of competition in- 
cluding VW and Toyota. This markest share is called the actual mortality 
market share, and it exhibits a cyclical trend with a period of some 
three years (corresponding, perhaps, to a three year automobile replace- 
ment period of customers?). A Theoretical Mortality Market Share (TMMS) 
was developed to conform to the force of mortality, h(t), obtained from 
the "best fit" life cycle results discussed previously. The assumed 
form of TMMS was as follows: 

(5) TMMS = A + B t 1,5 

Figure G-6 illustrates the "best fit" match of this theoretical mortality 
market snare to the data. The relative cyclical trend, AMMS/TMMS, is 
given in Table 2 and illustrated in Figure G-6. 

The numerical expression for the "best fit" TMMS is as follows: 

(6) TMMS = 150,000 + 65,530 (.25 t 1 * 5 ) 

The initial condition of the Mustang I life cycle occurs at t = 0 years 
(1964 3 I). Therefore 5 noting the initial value of TMMS, the difference 
impact of the competing market on Mustang is from Equation 5 s 

(7) TMMS - A = B t 1 * 5 

This represents the competing market share over and above the initial 
competing market share that Mustang I faced at the beginning of the 
life-cycle. Therefore, the force of mortality, h(t) is calculated from 
the equation 
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( 8 ) 


h(t ) = 0 . 25 t Xm0 


= <0.25/B)[TMMS(t) - A] 

= (0 *25/65530 )[TMMS(t) - 150000] 


While the cyclic trend is present, its effects tend to cancel and smooth 
out in the integration process indicated in Equations 2 and 3. 

The heuristic method used here in connection with Mustang I life 
cycle yields a result that can be interpreted: For Mustang I the in- 

novation rate remains essentially constant over its life cycle, while 
the force of mortality is proportional to the competitors total market 
share minus their initial market share at the start of the Mustang I 
life cycle. 


CONCLUSIONS 


' The heuristic shaping of product life cycles presented here generate 
a relative pair of countervailing forces: the force of product innovation 

and renovation, v(t), and the force of product mortality, h(t). These 
are actuarial concepts that can be related to the actual performance of 
a firm in a specified competitive setting. Ttey are sort of generalized 
performance measures potentially useful for ?jew products and strategic 
planning of a corporation. The illustration here was kept intentionally 
at a simplified level sufficient, however, for meaningful interpretations. 
The method can be useful for planning purposes provided that the analyst 
has a proper insight into the actual processes of a firm and understands 
the basic nature of its competitive environment. 






A major concern in the implementation of the total assessment pro- 
file relates to the identification of the regional scale at which a study 
would be carried out and the identification of key impact regions such as 
the total impact area, fringe impact areas and filter areas. It will be 
'the purpose of this Appendix to review the basic methodology geographers 
utilize in the delineation and analysis of spatial /temporal phenomena. 

The Appe- lix will represent a compendium of pertinent geographic method- 
ologie' ?n d their applications to the operationalization of the Total 
Assessment Profile. Through the implementation of these methodologies 
the researcher will be able to introduce and evaluate the spatial dimen- 
sions and impacts of new technological innovation. 


THE ANALYSIS OF NODAL DISTRIBUTIONS 


A major methodological thrust in geographic methodology is associated 
with the analysis of phenomena which exhibits a point or nodal character 
over space. Much of the characteristics of human settlement and activity 
patterns will reveal either an abstracted point configuration or an orien- 
tation to specific points in the space economy. An example of this, would 
be the location of human settlements, employment opportunities, and move- 
ment pattern orientations. Areas that have a concentration of such phenom- 
ena would be those likely to be hardest hit with the unslaught o s a natural 
disaster. This necessitates the development of methodology that can eval- 
uate the probable impact area and the extent of damage to the region . 

The geographic analysis of nodal patterns in space is normally struc- 
tured on a two-dimensional coordinate system. Within the coordinates each 
node is located in respect to an x and y coordinate position. The pre- 
ceding provides the basic information required to locate the node, define 
its relationship or spatial situation to other nodes in the system, and 
to study the aggregate relationships between subsets of nodes. Geographers 
have noted that every spatial arrangement of nodes possesses the following 
three characteristics: 

1. Shape - This quality was defined by Bunge [373] and is 
described as a closed curve which eliminates the collec- 
tion of objects and provides an areal measure of the 
distribution in two-dimensional space. 

2. Pattern - This is a zero -dimensional characteristic of 
~ spatial arrangement which describes the spacing of a 
sjI of objects with respect to one another [374]. 

3. Dispersion - This quality may be viewed as a one-dimen- 
sional characteristic of a spatial arrangement which 
measures the spacing of a set of objects in relation to 
one particular shape of a given area [375]. 
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THE ANALYSIS OF SHAPE IN A GEOGRAPHIC DISTRIBUTION 


The study of shape in geography has currently only limited methodol- 
ogy associated with it. The implications of shape, however, may be of 
substantial interest in technology assessment. Elongated political and 
cultural regions for example, provide significant difficulties in the 
provision and accessibility to essential services [376] . Service bottle- 
necks and administrative complication- may arise in such a regional con- 
figuration , Normative location theories for the allocation of industrial 
and service center locations over space have provided a comprehensive and 
convincing literature on the optimality of circular service regions. Hier 
archical models and the restrictions associated with spatial packing of 
service areas on triangular network lattices has forced a general compres- 
sion of these circular shaped service regions into a generally hexagonal 
network system. The economic optimality of such arrangements in conjunc- 
tion with the minimization of movement costs has been verified by an ex- 
tensive theoretical and empirically based literature [377]. The following 
figures show respectively a theoretical hexagonal configuration for settle 
ments and an adjusted empirically based service distribution system: 






An emergency logistics system structured on fast and error free services 
may be seriously hampered by an administrative regionalization that is 
inefficient in geographic structure. The full implications of shape, in 
particular the shape of administrative and functional regions is of concern 
to the geographic researcher. The measurement of shape, however, presents 
many difficulties. 


Most measurements of shape employ various indices that are ratios 
involving parameters of a distribution. The range for most indices fall 
between the value of zero and one. Departures from a line or circle tend 
to be the normative shapes utilized in these ratios : 





S 1 = 

(A/0.282 )P 

s 2 - 

(A/0.866L) 

S3 = 

R 1 ! r 2 

s 4 = 

(A/ir(0.5L) 

*5 - 

(1.27A/L 2 ) 


Where : 

A = Area of shape being measured 
L = Length of Longest Axis 
Rq = Radius of largest inscribing circle 
R 2 = Radius of smallest circumscribing circle 
p = Length of perimeter 


In all of the above ratios a value approaching zero would suggest a linear 
elongated shape, while a value approaching one would imply a circular shape 
to the region in question [378]. 


Two other approaches have been commonly used in geography to define 
the quality of shape. The Boyce-Clark [379] method is based on a formula 
that calculates an index which varies between 0 and 175 regardless of size 
of the area concerned. The index is represented below: 


Where : 

SBC - Boyce Clark Shape Index 

SBC = £ |(rj/£r£)l00 - 100/n | r^ = the length of the ith radial 

n = number of radials 


A major problem with the above shape index is the determination of the 
point within the shape from which the radials are to be drawn. Most re- 
searchers use the center of gravity as this orientation point. Interpreta- 
tion of the ratio does provide some problems with the identification of 
ground truth. An approximate shape correlation with the Boyce-Clark scale 
is provided in the following: 

0 12 18 25 28 175 


Square Star Rectangles 

Cruciform Rectangle (x>y) 

(2x+y) 

Bunge provides an alternative measure of shape that is much more com- 
plex than the previous ratios and indices [373]. The method Is based on 
two theorems : 

1. Any simply connected shape can be matched by a polygon of 
any number of sides in which the lengths of each side are 
equal, but the lengths can of course vary from shape to 
shape . 

2. If the distance between all vertices of the polygon lag 1 
are summed and then squared and summed, lag 2 are summed 
and then squared and summed, and so on, there will exist a 
unique set of sums that define the shape of the polygon. 

Thus, the first step, in the Bunge method is to define for any given shape 
an abstract polygon which adequately represents the shape to be defined: 
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Actual Shape 


Abstract Polygon 


The unique set of sums which describe that polygon can then be calculated 
with respect to its vertices. All of the above methods have computational 
problems associated with them as well as interpretation problems with de- 
fining the ground truths of the index values. In a cost benefit model the 
shape indices would be employed with "hard data" to define the degree of 
relationship and its respective implications. 


THE ANALYSIS OF SPATIAL PATTERNS 


The analysis of spatial patterns , is of perhaps s greater significance 
in a cost benefit application than the function of shape. It Is through 
the use of patterns that spatial relationships may be defined. The dif- 
fusion network , its areal extent s and its configuration would be defined 
by the concept of pattern. Point pattern analysis has been of major 
methodological interest to geographers. The general orientation of these 
methods is to establish an observed set of points and compare these points 
to a theoretical set of points that would be generated by a number of prob- 
ability processes. A basic division of pattern classification may be 
structured from the utilization of probability distributions. Point pat- 
terns will either reflect a random spatial arrangement or will exhibit a 
bias spatial pattern that would reveal a vacancy or occupancy bias; 


QUADRAT 




This figure demonstrates a classification structure for a point distri- 
bution analysis. A major division in the taxonomic chart is reflected in 
the determination of a pattern reflecting a random arrangement opposed to 
the bias arrangements. To define the implications of randomness in space 
the following conditions are assumed: 

1. Any point on the surface has equal probability of occur- 
ring at any position on the surface. 

2. The position of a point at the pJbove place is independent 
of any other point. 

A major issue in the analysis of point distributions is based on 
identifying whether the point pattern reflects a regular,, random s or clus- 
tered spatial distribution. Two basic methodologies may be employed to 
accomplish the above identification: quadrat analysis* and point to point 

analysis. In quadrat analysis a grid is laid over the study region with 
a frequency count of points contained in each cell being calculated: 
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The observed frequency distribution is then compared to a theoretical one. 
For example, a Poisson distribution may be used in the case of randomness, 
or a negative binomial function for a clustered distribution. 

Point to Point distributions, on the other hand, have been developed 
from nearest neighbor statistics. This method involves computation of 
the mean distance and associated variances for each order of nearest neigh- 
bor points. The preceding are th^n used in a comparison to define the 
parameters that would character the expected distances : 

Where : 

« R = Nearest Neighbor Statistic 

R = (Er/n)/(l/2(/n/A) r = Observed Nearest Neighbor Distance 

n = Total Number of Nodes 
A = Area of the Study Region 

Thus , under the assumption that the first order distances are drawn from a 
normal population, a density dependent expected mean would be defined. 



Randomness can then be tested using the standard normal curve. Alterna- 
tively a ratio of observed and expected mean distances can be computed 
(R, the nearest neighbor statistic). Employing the nearest neighbor 
statistic the following range in values may be derived [380,3813: 

0 1.0 2.14 

Absolute Random Absolute 

Clustering Regularity 

The lslationship of the spatial pattern to a cost benefit analysis 
of an early warning disaster system rests with the extent, orientation, 
and duration of the impacted region. If the distribution of settlements, 
for example, indicate a clustered pattern the expected impact would be 
far greater than may be anticipated from a region in which the settlements 
reflect a regular or random spatial arrangement. It is also through the 
measurement of pattern that disaster agents and their expected impacts may 
be isolated. Disaster events that reflect a random spatial occurrence 
would prove to be more difficult to develop an adequate warning network 
compared to those that reveal a definite and predictable spatial pattern. 


THE ANALYSIS OF SPATIAL DISPERSION 


Closely related to the element of pattern, is the concept of spatial 
dispersion. Dispersion, is perhaps, the most important spatial concept 
to be considered in introducing the areal dimension to a cost benefit anal- 
ysis. It is by means of this dimension that the areal dispersion, and to 
a great extent, the necessary intensity of impact over the earth's surface 
may be isolated. The relationship of dispersion to the preceding concepts 
may be defined as an attribute of a pattern that is located within a par- 
ticular shape at a given distance. 

The analysis of spatial dispersion draws heavily on the utilisation 
of centrographic statistical procedures. The centroid or arithmetic mean 
center of an areal distribution {bivariate distribution) is analogous to 
the concept of the arithmetic mean of a linear distribution. Neft [3823 
indicates that virtually all arithmetic mean centers have been calculated 
as the point representing the arithmetic mean of the X values and the 
arithmetic mean of the Y values where X and Y are a pair of orthogonal 
axes : 

Z X£ _ 2 Yi 

x = n“ y = if 


The standard deviation of a linear distribution has for a counter- 
part in bivariate statistics the standard distance deviation. Hence, 
this statistic represents a dispersion measure over space, Lee [3833 
notes that the standard deviation is minimized about lines Xq and Y g 
drawn through the mean centers of the X and Y reference axes: 
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Standard Distance Deviation X 



Standard Distance Deviation Y 


Normally the standard distance deviation about each of the two orthogonal 
axes will not be of the same length except in the case where a distribu- 
tion is circular in shape. 


A measure of spatial dispersion computed along the axis of the dis- 
tribution is referred to as the standard distance or the standard radius 
of the distribution. The standard radius is defined as the hypotenuse of 
a right triangle formed from the sides of the distance variance of the X 
and Y axes; 


Standard Radius (SR) 



The combination of the arithmetic mean center with the standard distance 
deviation, distance variance, and the standard radius provide a quantita- 
tive expression of an areal distribution. 


The standard distance deviations will vary in their orientation to a 
set of perpendicular central axes which are parallel to the X and Y refer- 
ence axes. The central reference axes do not necessarily maximize or min- 
imize the standard distance deviation. Lee indicates that a set of central 
axes called the principal axes can be found by rotation which will maximize 
the standard distance deviations. The principal axis about which the 
standard distance deviation is minimized will be referred to as the major 
principal axis. Consequently the central axis perpendicular to the major 
principal axis that maximizes the standard distance deviation is called 
the minor principal axis. The above minimum and maximum standard distance 
deviations about the central reference axis can be used to provide an in- 
dex of circularity. The coefficient of circularity is simply a ratio of 
the standard distance deviations about the major principal axis to the 
standard distance deviation about the minor principal axis : 


Coefficient of Circularity (CC) 


q yp 

°xp 


Where: Oyp = Standard distance deviations about the 

major principal axis 

a X p = Standard distance deviations about the 
minor principal axis 

If the preceding ratio is equivalent to one this would imply both the 
maximum and minimum standard distance deviations about the principal axis 
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are of the same length, thus the distribtuion would be circular in shape. 

As the index approaches zero the standard distance deviation about the 
minor principal axis would become very small suggesting a linear pattern 
in the distribution. In most cases the distributions of spatial phenomena 
are not evenly dispersed; thus an elliptical shape will characterize the 
majority. The preceding is provided by a graphic device called the Standard 
Ellipse. The larger the ellipse the more dispersed will be the phenomena 
being studied. On the other hand, a small circular ellipse would indicate 
a concentration of the phenomena in space. The following illustrates the 
ellipse and the general structure of a centrographic analysis of a point 
distribution : 



THE ANALYSIS OF SPATIAL NETWORKS 


The preceding discussion has provided a means of identifying the 
nodal intensity and distributional characteristics for a regional system. 
Included within the above system, but not yet introduced would be the 
connecting networks between points and phenomena demonstrating a nodal 
character. The demand and level of interaction between nodes would ex- 
hibit over space a meshed fabric of intersecting lines and routes . The 
intensity and connectivity of these networks would indicate the level 
of potential interaction and the degree of redundancy in potential cir- 
culation circuits. Both of the above would be critical considerations 
within a cost/benefit analysis of technological innovation. Specifically, 
in the case of a disaster warning system, the nature of the existing net- 
work may be a critical factor in defining the temporal lag be Lween the 


warning sequence and the probable response time for evacuation and emer- 
gency aid. 


The geographic analysis of the spatial structure of networks has 
primarily utilized finite graph theory as the key methodology* The geo- 
metry of a graph network system is reduced to an abstracted finite set of 
nodes (V) and to a possible empty set of edges (E). Hence, graphs are 
defined as sets of points or nodes which may or may not be connected by 
edges or links to one another. It is through the use of graph theory that 
the basic topological structure of a network may be defined [384.385,386]* 
Three major types of networks appear to be relevant in a cost benefit anal- 
ysis centered on a disaster warning system: 

1* Branching Networks 


The study of branching networks has its applications focused primarily 
in the area of the earth sciences* In finite graph theory a branching net- 
work represents a minimally connected graph in which all nodes are connected 
however, no circuits or loops exist in the system* The major emphasis on 
the study of branching networks, or what are sometimes referred to as tree 
graphs, are in applications associated with the configuration and operation 
of river drainage systems. 

A stream network, for example, represents one of the simplest config- 
urations from a topological perspective* All bifurcating networks are, in 
an abstract sense, topologically identical regardless of their unique at- 
tributes. An example of a simple branching graph is provided below: 


Vf 



Where: - Internal Edge 

E e - External Edge 
V t - Terminal Node 
V-^ - Branch Node 
V r - Source Node 


The symbolization associated with the graph is derived frcm Melton’s [3873 
research on drainage system configurations are components . The geometry 
of a bifurcation graph permits the identification of the following rela- 
tionships between nodes and edges [387,3883: 

E = (V t + V b + V s ) - 1 
E = 2V t - 1 
V b + V r = V t 


89 


* 

Based cr the preceding Shreve [389] has developed a procedure to describe 
the topological structure of the tree network. The method simply involves 
tracing a systematic route along a graph while recording all interior and 
exterior links traversed. Each time an interior link is traversed a +1 
is recorded and each exterior link encountered is recorded as a -1: 



An important consideration in the research associated with branching 
networks has been oriented to the differentiation in edges based on a hier- 
archical order. This concept is of particular importance in the study of 
a drainage system * since each segment may exhibit different morphometric 
and hydrologic features. Most of the hierarchical models assign a level 
of relative importance to each segment in a branching network. The approach 
is based on a sequential .arrangement of tributary and major stream segments. 
An example of a hierarchical ordering system is illustrated in the work of 
Horton [390] and in the modifications provided by Woldenberg [388], Hor- 
ton^ approach defines the first order stream as one receiving no tribu- 
taries with the second order stream as one receiving two first order streams. 
Subsequent orders to stream segments are then based on the junctions of the 
previous order segments: 



This diagram is a modification of the Horton system by ffoldenberg [388]. 
Several other systems of hierarchical ordering methods are available and 
it is suggested that the interested reader review the work of Strahler 
[391] , Melton [387] 3 Wamtz [413] s and Woldenberg [388]. 

A review of this literature will define two major research implica- 
tions: (l) the significance of the relationship between different orders 
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of stream segments, and (2) the relationship of the various orders of 
stream edges to the entire drainage network. Haggett [392] defined the re- 
lationship between the different orders of the network edges as the "Law 
of Path Numbers". This concept suggests that the number of edge links for 
different orders will approximate an inverse geometric series in which the 
first term equals unity, while the ratio between each order in the sequence 
would represent the bifurcation ratio. Individual orders of bifurcation 
ratios may then be used to define a mean bifurcation ratio (the aggregate 
ratio may be defined by a simple arithmetic mean, a weighted mean, or a 
geometric mean) for the entire network: 



The study of bifurcation ratios has susggested a relatively high sta- 
bility between regional drainage systems. The minimum bifurcation ratio 
between 3,00 and 5.00 is seldom found in nature for it requires a basin 
topography of relative homogeneity. It is through the utilization of the 
bifurcation ratio and its relationship to the order of the network that 
descriptive relationships may be defined. 

The application of branching networks to a disaster warning system 
analysis are primarily related to hydrologic effects during the impact 
phase. The simulation of water flow through a hierarchically ordered 
drainage system determines warning time and necessary lag times for an 
effective response. The applications of the study of branching networks 
has several implications to the human interaction system. Evacuation 
and recovery efforts may reflect over space a tree like collection or 
dispersal system. The identification of minimal travel cost spanning 
paths would be critical in defining major emergency escape routes. It 
would also be essential to define high risk links within the graph that 
could function as potential bottlenecks or breaks which would inhibit or 
prevent movement over the system. Most human interaction systems, how- 
ever, are more complex and the branching network would merely represent 
a subgraph or partial graph of the more complex system. 

2. Circuit Networks 


A circuit network would represent a more complete graph that con- 
tains closed loops and circuits. Attempts at measuring the structural 
character of circuit networks have resulted in a basic set of metrics. 
The cyclomatic index, for example, defines the number of circuits within 
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the graph. The alpha index is a ratio measure of the number of actual 
circuits to the maximum number possible in a given network. The third 
measure is the Gamma index which is also a ratio which compares the num- 
ber of existing linkages to the maximum possible. Thus, the cyclomatic 
number, the alpha index, and the gamma index provide a set of measures 
of the degree of connectivity and circuitry which characterizes a par- 
ticular network: 

Cyclomatic Number: 

C = E - V +’ G 

Alpha Index: 

a = '(E - V + G)/(2V - 5) 100 For Planar Graphs 
a - ( (E - V + G)/(V( V - l)/2 - (V - l)) 100 For Non-Planar Graphs 


Gamma 

Index: 


Y = 

(E/(V(V - l)/2) 100 

For Planar Graphs 

Y = 

(E/(3V - 2)) 100 

For Non-Planar Graphs 


Where : 



C - Cyclomatic Number 

E = Number of Edges 

V - Number of Nodes 

G = Number of Subgraphs 

The combinations of the Gamma and Alpha indices may be used to define the 
degree of complexity characterizing the network. Taaffe and Gauthier f393j 
indicate three major stages in network structure: spinal networks, grid 

networks, and delta networks. The spinal network represents the previous- 
ly discussed minimally connected graph or tree, while the grid and delta 
networks characterize intermediate and maximally connected graphs respec- 
tively. According to Taaffe and Gauthier T393*] the following breaking 
points may be defined for the relationship of the Gamma and Alpha indices 
to the level of the network: 

Gamma Index Alpha Index 

Spinal Network .33 y .50 Where v >_ 4 a - 0 Where v = e + 1 

Grid Network .50 y .66 Where v > 4 0 a .50 Where v 3 

Delta Network .67 y 1,00 Where v > 3 .51 a 1.00 Where v > 3 
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Hence, by using the above measures in association with one another the 
degree of connectivity and circuit redundancy may be defined for a net- 
work. 


The implications of the level of network development as measured by 
the above provides a framework to define the configuration and complexity 
of an existing transport or communications network. Individual roads and 
their relative locations to the graph are measured through the utiliza- 
tion of incidence and connectivity matrices. The simplest connectivity 
matrices are binary in structure and symmetrical. The elements of the 
matrix defines the existence or lack of existence of a linkage (unity is 
used to indicate a link while zero indicates the lack of a linkage). By 
powering the connectivity matrix both the direct and indirect connections 
between nodes may be defined. 

An analysis of the incidence matrix for the networks will provide an 
initial measure of accessibility. The incidence matrix is an n by m rec- 
tangular array with the columns of the matrix representing the edges and 
the rows associated with the network nodes: 



The elements of the preceding matrix are equivalent to one if an edge is 
incident to the node (a^. = 1 when v^ is linked to e) or zero if a lack 
of incidence exists. Multiplication of the incidence matrix (A) by its 
transposed matrix (A^) will give a connectivity matrix (C) which is sym- 
metrical: 

T 

n^m • tiAi - n^Ti 

Garrison [415] suggests a removal of the positive diagonal elements of 
the connection matrix replacing them with 0 ' s based on the rationale that 
a reflexive incidence structure at a node is meaningless. 

Powering the connection matrix to the diameter of the network and 
then summing the rows of tlie powered matrices will provide a measure of 
node accessibility for both direct and indirect linkages: 

nCn + nCn^ + nCn^ + + nlhi^ = nTn 

Where C = The Diameter of the Network (£ = max. d(v^vj) 
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The above nodal accessibility metric has significant shortcomings 
because it includes redundant paths generated by powering the connection 
matrix. Removal of the redundant paths in matrix T is mathematically pro- 
hibitive, thus creating an additional complication to the problem. A 
potential solution to the above is provided by a technique developed by 
Shimbel [394], Shimbel* s approach consists of recording the power, in a 
distance matrix (D), at which a one enters the cells of the powered con- 
nectivity matrices. Thus, if a one enters element C 44 on the third power- 
ing of the connection matrix a three is recorded for cell C44 of the 
distance matrix: 



Where p = the power at which a one enters the 
connectivity matrix 

It may be noted that the diagonal remains zero with the off diagonal el- 
ements indicating minimal path structures in the network. 

Summing the rows of the Shimbel distance matrix provides a measure 
of node accessibility. A high row sum would indicate that linkages to 
other nodes in the system are indirect and do not occur until the higher 
powers of the connection matrices are reached, thus the node may be con- 
sidered inaccessible to the system. 

The preceding techniques define nodal accessibility on the basis of 
equal length linkages between nodes. All linkages in the preceding were 
assumed to be of equal distance or travel time. The preceding employed 
a shortest path matrix procedure in which the distance metric was purely 
topological in structure. Distance was defined in the above simply as 
the number of linkages between nodes. Hence, all linkages were assumed 
to be of the same value. If additional information regarding linkage 
structure is available the network may be studied as a valued graph. In 
the latter, linkage distances may be weighted by actual mileage or travel 
time. A similar procedure to define nodal accessibility may be used in 
a weighted or valued graph approach [393], The initial step is based 
on the construction of a connection matrix in which the cell entries are 
not binary data but rather the weighted values of the linkages. Where 
linkages are absent infinity is used to define the value of the linkage, 
while direct links are recorded as their actual distance value. Diagon- 
al cells remain as zero cells in the connection matrix: 
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This matrix provides the direct linkages between nodes in the network. 
Indirect linkages are determined by matrix powering procedures similar to 
the one used to define accessibility in a binary connection matrix. Un- 
like the previous multiplication procedures Boolean mathematical rules 
are employed. Instead of element by element multiplication of the rows 
times the columns of the connection matrix we utilize an element by ele- 
ment addition. Rather than summing the results, we use the minimum value 
and insert this minimum value in the cells of a new matrix [393]. Hence, 
the cell ij value in the new matrix is not the sum of the products of all 
links but is rather the minimum values of the sums of links from the ori- 
gin to the destination nodes. An example of this is provided in the study 
by Taaffe and Gauthier [393] on the following page. The matrix L repre- 
sents a valued connection matrix with each power of the matrix defining 
the two, three, and four step indirect linkages respectively. Summing 
the rows of the fourth power linkage matrix (L^) provides a measure of 
nodal accessibility of Individual nodes to the entire network. 

This study represents a static treatment of network structure. An 
important consideration in the spatial organization of a network may be 
related to the dynamics of flows through the system. These flows may 
represent population migration and Interaction, the communication of ideas 
or messages, or the transmission of goods and services over space. The 
basic typological elements discussed earlier are important considerations 
in defining a networks flow capabilities. The static typological treat- 
ment of graphs must be modified so that dynamic concepts such as direc- 
tion and velocity may be included in addition to the analysis of connec- 
tions between nodes. A digraph is constructed where the linkages are rep- 
resented as having a direction orientation. In the binary case, a +1 in 
the incidence matrix would represent a destination node with a -1 Indicat- 
ing an origin node. Due to the limitations in length of this paper an 
extensive treatment of flow problems will not be undertaken. The general 
structure of these problems, however, will be briefly discussed. 

An operational format for a typical network flow problem is presented 
in the following diagram: 

Network Intermediate Network 

Source Feeder Network Collector Sink 

Node Nodes Nodes Nodes Node 
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In the' diagram the source and sink nodes are used to interject a flow in- 
to the feeder nodes and to recover it from the collector nodes of the 
network. The digraph may or may not be symmetrical with feeder, collec- 
tor, and intermediate nodes having one-way or two-way linkages. 

A study of the circulation through a network may have several objec- 
tives. A general application of a flow problem would be one oriented to 
defining a minimal travel spanning tree through the network. Several 
algorithms are available which determine the minimal paths within networks. 
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A second application is associated with problems of maximal flow in a net- 
work system. In this case, constraints of capacity are interjected on 
network linkages and sometimes on the nodes themselves. A third major 
application is associated with combinations of the above. A minimal cost 
maximal flow problem, for example, would be illustrative of this type. 
Several labeling algorithms are available for the solution of the minimal 
cost problem having maximal flow constraints [395], 

The study of circuit networks and the treatment of flows over these 
networks are critical elements in the evaluation of costs and benefits 
over-space. The potential capability of an emergency delivery system 
would be constrained by the network’s configuration and its linkage and 
nodal capacities. Changes in the emergency warning technology would re- 
quire an evaluation of its impact on the communications and transporta- 
tion networks. 

3 . Barrier Networks 


Barrier networks reveal a substantial difference in function from 
the previous two types . Branching and circuit networks consist of chan- 
nels that conduct communications, goods, and people over space while 
branching networks consist of links which either block or resist flows. 
Chorley and Hagget [386] observe that such barrier systems demonstrate 
aspects of closure by consisting of closed loops which are isolated or 
contiguous. We may regard barrier networks as the ’’duals of flow networks 
in both a functional sense and in a topological sense”. To illustrate the 
above point the following diagrams were taken from a study by Warntz [396 



The cellular nets are a form of graph theory that is associated with 
planar, graphs (two-dimensional graphs). They are what Ore [397] calls 
polygonal graphs in which the edges form a set of adjoining polygons on 
the plane. An example of such a polygonal plane vrould be a map of U.S. 
counties. Euler’s formulation of the relationship between edges and 
vertices in a simple polyhedron establishes the foundation for the geo- 
metric structure of polygonal graphs: 
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V - E + F = 2 Where: V = Number of Nodes 

E = Number of Edges 
F = Number of Faces 

A polyhedron is defined as a solid whose surface consists of a number of 
polygonal faces. A simple polyhedron has no "holes" in its surface, so 
that its surface is continuous. An example is provided below: 


V - E + F = 2 

8 - 12 - 5 - 6 = 2 


Euler's formula permits one to establish the conditions for the distri- 
bution of regular polygons over space. This condition is defined when 
the following relationship is fulfilled: 

P P 

1 + = — Where: P = Number of edges at each node 

^ P* = Number of edges bounding each 

cell 

Chorley and Fngget observe that this relationship will only work when 
both the number of edges at each vertex and the number of edges bounding 
each cell are equal. These conditions would be demonstrated in regular 
polygons such as triangles, quadrangles, and hexagons. In a regular hex- 
agonal pattern, for example, the contact number would equal six since 
each cell would be contiguous with six neighboring cells. 

The implications of the barrier networks or polygonal graphs rests 
with problems associated with the partitioning of space. Efficient ser- 
vice and administrative districting plans would be dependent on an optimal 
spatial allocation system. The preceding permits the researcher to evalu- 
ate an existing or projected areal arrangement of service and adminis- 
trative districts according to basic geometric concepts of efficient space 
packing and partitioning. Inefficiencies in the configurations of service 
hinterlands would reflect a cost to the delivery system. An emergency 
warning technology that would be dependent upon reaching a maximum num- 
ber of people in the shortest period of time requires the consideration 
of geometric and areal dimensions. 

The preceding provides the basis for the analysis .of the spatial 
infrastructure. The methodologies noted earlier define the regional geo- 
metry and its applications to the distribution and interaction of loca- 
tions in space. Critical to the application in this Appendix is the 
identification of aggregate patterns in the construction of regions. It 
is through the concept of a region that the key impacts and interactions 
of an emergency disaster warning system may be evaluated. The region 
provides a spatial benchmark through which pertinent information may be 
categorized and later evaluated. In addition, it gives the researcher 
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an areal unit in which accounts of benefits and costs for the new system 
may be measured. It is critical that the identification of meaningful 
planning regions be a central issue of the temporal/spatial dimensions 
of a cost benefit analysis. 


REGIONAL IDENTIFICATION 


Regions appear to be one of the most logical ways of recording and 
evaluating geographic information. Geographers have defined several major 
categories of regions. One of the most popular typologies was developed 
by Whettlesey [398] and is presented in the following diagram: 



Single Feature 
Region 


Areal Classes 


Multiple Feature 
Region 


Modal Regions Uniform Regions 


The first level of regional classification is based on the number of 
features or phenomena considered in the definition of the region. A re- 
gion is defined as a human construct in which an area of space is found 
to be uniform or homogeneous in relation to a feature or group of fea- 
tures. The notion of a human construct is important for regions are 
mental images and are real only to the extent of their applications. In 
a cost benefit application we are dealing with a potential impact area 
f in which a multiplicity of phenomena would be of interest. This applica- 
tion requires that a multiple feature impact zone be defined in the in- 
itial stages of the cost benefit analysis. The degree of homogeneity 
and its distributional characteristics are basic considerations at the 
second level of the regional typology. Uniform regions are those in 
which the phenomena is spread over the area at approximately the same in- 
tensity. This differs fipom a nodal region where homogeneity is not based 
on intensity, but rather on orientation to a node or point. 

Most human interaction systems reflect a nodal regional character 
in which activities are focused on a core or central area. In the nodal 
region intensity of activity or orientation appears to weaken with in- 
creasing distance from the core. The majority of applications in the 
cost benefit analysis will deal with the interface of both nodal and un- 
iform multiple feature regions. 

Uniform regions reflect the degree of overall impact and would be 
utilized in defining a threshold level. All areas that fall above the 
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threshold would represent the impact region to be considered in the anal- 
ysis. Critical to the utilization of the uniform regional concept is the 
definition of the intensity level or threshold. The definition of risk 
versus non-risk areas, impact versus non-impact, fringe impact versus 
non-fringe impact, and filter areas versus non-filter zones are all prob- 
lems of regionalization which require the identification of all contiguous 
locations experiencing a specific intensity of a phenomena. 

Nodal regional systems relate to the orientation of activities and 
are probably more relevant when viewing the delivery and interaction as- 
pects of technology. Both perceived and objectively determined regional 
orientations would be significant considerations in the evaluation of an 
early warning emergency disaster network. It may be noted that a stoch- 
astic element may exist within the regional identification problem. This 
is demonstrated in the probabilities and uncertain it ies associated with 
the identification of risk regions and with the orientation of interactions 
within nodal regions. 


APPROACHES TO REGIONALIZATION 

I. REGIONAL IDENTIFICATION BY BOUNDARY DEFINITION: 


The procedures employed in regional boundary definitions is a problem 
in the delineation of the limits of a phenomena and its spatial configur- 
ation. Fletcher [399] suggested that set theory could be employed in the 
problem of regional delineation. An example of Fletcher's approach is 
illustrated below: 




A jf B • 


The Venn diagrams A and B reflects a collection of counties associated 
with two states: A and B. In the above regional format a simple loca- 
tional discrimination problem between two mutually exclusive groups is 
defined. In the following example a set of hierarchical relationships 
are demonstrated: 


In the diagram x represents a set of counties which are included in state 
B which is part of a nation A, Hagget [392] notes that this concept of 
set theory demonstrates that maps are a special type of Venn diagrams. 

Problems of regional delineation may be solved by two maior approach- 
es; qualitative and quantitative. Qualitative approaches to regional 
delineation are based on intuitive decision processes. The Midwest, for 
example, is a regional construct that intuitively may demonstrate sub- 
statial variability in its perceived location. Several problems in the 
perception literature relate to qualitative delineations of regions. The 
map below represents & study by Cox [400] in which students were asked 
to define the location of specific cultural regions in the U.S. 



The most serious problem of qualitative approaches to regionalization 
occurs in the identification of boundaries between regions. The over- 
laps in regions shows confusion among the respondents and indicates 
their failure to generate mutually exclusive regional concepts. By the 
superimposition of boundaries , core areas of regions are usually defin- 
able, however, a serious element of uncertainty remains in the location 
of boundaries . 

Quantitative analysis for regional delineation appears to be more 
promising than qualitative attempts. Several approaches to regional- 
ization may be found in the current geographic literature: 
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1 


Distance - Minimization Functions 


This approach utilizes methodology from the operations research lit- 
erature associated with the solution of the transportation problem. Demands 
and supplies are allocated in a manner which minimizes the interaction 
costs of the following objective function: 


n 

ill 


m 




minimum Where : 

d^j = distance 
X£_: = units being 
allocated 


This type of analysis is useful in defining optimal spatial allocations of 
nodal regional systems. King, Muraco, and Vezner [4013 provide an example 
of the distance minimization model in a problem associated with the region- 
al allocation for mental health care service centers. The diagram on page 
103 shows a minimal distance solution to the regionalization problem. On 
page 104 a modified version of the model is demonstrated with weighted 
linkages based on a predicted case load over space. The regional differ- 
ences in the proposed catchment areas (service regions) comparing the 
weighted versus the distance minimization model is provided on page 

2. Discriminant Analysis 

Discriminant analysis may also be used in the. regional identification 
problem. Discriminant analysis is used to classify a set of observations 
which have already been classified according to some categories. The 
techniques associated with this form of analysis are drawn from the work 
of statisticians concerned with biological and anthropological data. The 
analysis attempts to define or classify observations which could conceiv- 
ably be associated with several classes. 

The approach rests on taking observations which may have arisen from 
several populations and differentiating them in relation to one another. 

The geographic analogy to this approach would be the identification of 
two core regions and an intermediate area which may he associated with 
one or the other of the two core areas. 

A solution to this problem is provided by the generation of a single 
linear discriminate function for the p variables so as to ensure maximum 
discrimination between the two populations while minimizing the probabil- 
ity of error in the assignment of new cases to one of the two populations 
(core areas). The more general problem may involve more than two regions 
and require the derivation of a multiple linear discriminant function. 

King [4023 notes that two major considerations underlie the utiliza- 
tion of discriminant, analysis: 

- a pre assumption that some form of classes are given. 
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- the assumption that individuals or objects to be classi- 
fied belong to one or the other of the classification 
categories. 

Fisher [414] states the problem as one of derivation of linear functions 

for a set of observed variables: (X^ 5 X 2 ,X 3 X~) where Y = X^X^, 

X 2 X 2 ,-.. XpXp. The objective of the analysis is to maximize the 

difference of the ratio between the sample means to the standard error 
within the two samples. 

King [402] provides an operational solution to the above problem 
where N observations are divided into two groups or samples of ra a and ng 
members respectively. Each observation in sample a is given a Y score 
of (n a /n), whereas the numbers in sample 3 are given the value (~n a /N) . 
According to King, these Scores are arbitrary and assure that the mean 
of the Y scores are zero. 

The discriminant function: Y = b^X-^ + b 2 X 2 + + bpXp represents 

the multiple regression equation with the constant omitted. The b values 
in the above equation are obtained by the use of the least squares tech- 
nique. The mean predicted Y score may be obtained by the utilization of 

(Xi,X 2 , X ) for each group. New individuals are then classified 

on the basis of their Y scores with reference to the mean values. 


The frequency misclassification is defined by the use of a t-statistic 
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King suggests that the significance of the discrimination may be tested 
using an analysis of variance model: 

Where : 

F = (SS T /p)/(SS w /N-p-l) SS-. = Within Sum of Square 

SS^ - Total Sum of Square 
p = Number of Variables 
N = Number of Obs 

Discriminant analysis is similar to principal components analysis in 
that it allows the researcher the capability of identifying the basis for 
the groupings. Casetti ' s [403] work in discriminant analysis has pro- 
vided much of the basic programming for regional classification applica- 
tions. His contributions are associated with multiple discriminant iter- 
ations and the means by which they are used to force classification of 
an optimal solution. The normative state is defined when within group 
variance is minimised and between group variance is maximized. 


3. 


The Utilisation of Graph Theory 


The growing importance rf graph theory in the study of geographic 
structure has resulted in the emergence of several graph oriented methods 
of regionalization. The majority of the graph approaches to regionaliza- 
tion are based on hierarchical models of interaction networks. Nystuen 
and Dacey's [4041 study of trade flows between cities in Washington re- 
flects the graph approach to regionalization. The following table in- 
dicates a hypothetical flow arrangement: 
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Within the table four dominant trade centers are identified. Associated 
with each of these four are a set of satellite communities. These obser- 
vations are based on the size and direction of the flows within the 
transaction matrix. A recent extension of the preceding approach may be 
found in the research of William Black [405 1. Black utilizes a dyadic 
factor analysis to extract the basic dimensions of a set of interregional 
commodity flows. 

The preceding approaches define regions by taking defined locations 
and clustering them so as to maximize or minimise basic distributional 
components. Distance minimization, for example, may be used to define 
the boundary separating two distinctive nodal regions. Through the use 
of discriminant analysis, components associated with the variance of 
clusters of observation may be utilized to obtain an optimal structure. 
More recently, the graph theoritic approach has provided a means of 
regionalization based on functional interactions between locations and 
the assumption that a hierarchical structure exists. The traditional 
non-quant itative approaches also suggest several new directions to the 
problem of regionalization. As the paradigm of the behavioral environ- 
ment becomes a dominant research theme the measurement of cognitive 
regional structure becomes a central issue. Social and cultural affini- 
ties demonstrate a spatial component [406] that can only be isolated 
through measurement of cognitive spatial dimensions. 


The problem of regional identification as treated in the preceding 
methodologies assumed that the researcher had complete control over the 
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basic areal units from which the regions would be derived. Often this 
is not the case and the geographer must utilize data that is already 
aggregated. When the data is already defined by areal units the problem 
of isolating regions becomes an assignment task. The researcher is 
forced to cluster these pre-defined areal units so as to create meaingful 
classifications that would represent true regions. The following section 
examines the approaches that may be utilised in assigning areal units to 
regions . 


II. THE REGIONAL ASSIGNMENT PROBLEM: 


In the assignment problem, regionalization becomes primarily a process 
of classification. Berry [4071 defines four basic principles which under- 
lie the regional assignment problem: 

- every location must ultimately be placed in a region 

- only under special conditions will the assignment problem 
not be mutually exclusive. 

- regions may be defined by successive splitting of bigger 
• regions into the smaller regions. 

- larger regions will exhibit greater internal variability 
than the smaller regions. 

Berry observes that the methodology of regionalization is based on defin- 
ing the degree of similarity for each pair of places . These places are 
then grouped with other pairs with the objective of maximizing the similar- 
ity. The above is subjected to the constraint that the groupings of places 
should be contiguous. 

Several methodological approaches may be utilized in the assignment 
problem. The major regionalization procedure utilize the following com- 
binations of techniques: 

- Factor Analysis - A multivariate approach which permits 
the researcher to define major components of the corre- 
lation matrix and to establish the areal structure of 
these components. It is from the areal distribution of 
the factors or dimensions that regions may he defined. 

- Dimensional Analysis - This- approach utilizes the pre- 
ceding factor output or may work with empirically 
derived data to define the degree of similarity of each 
location by its attributes. The multidimensional scaling 
approaches apply this basic concept and extend it into 
complex and sophisticated research designs [408 ] . 
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- Grouping or Cluster Analysis - It is through this technique 
that the observations are clustered into regions of con- 
tiguous spatial sets. 

The following section of the paper will briefly examine each of the above 
methodologies and indicate how they may be used in a regionalization prob- 
lem. 

1, Factor Analysis 

A variety of research designs may be defined under the general title 
of factor analysis. All of these models utilize the relationships that 
exist between sets of variates to define their underlying roots. The 
models take a large number of variables and attempt to define a composite 
set of variables that are fewer in number while minimizing the variability 
of the original relationships. 


THE BASIC MODEL 


The mathematical structure of the factor analysis problem applied to 
a spatial case may be demonstrated by the following data matrix: 

\ 12345 6 a 


Data Matrix 


In the above, matrix A represents an n by a matrix m which a set of 
attributes (variables) a and a set of observations (areal units) n^ and 
their associated elements (magnitude of a at n) are presented. The at- 
atribute set, a, may represent physical, environmental, or cultural qual 
ities which are the characteristics to be regionalized. The observation 
set n, on the other hand, are the areal locations which may be represent 
ed as townships , counties, states, census tracts, nations, etc. 

A principal axis factor model would then take the above matrix A 
and transform it to a matrix Z in which the original variables would be 
normalized : 


n A a 


n Z a 




The above normalized matrix (standardized data matrix V would then be sub- 
jected to a correlation analysis to generate matrix R: 

n Z a Correlation Analysis ^ a R a 

The symmetrical matrix a R a would represent the zero-order correlations 
among the standardized variable set. The correlation matrix would then 
be subjected to a principal axis factor analysis of (R - tr) in which 
would represent a diagonal matrix which would contain the unique variance 
of each of the original a_ variables, The result of the above produces a 
matrix F of order a by s in which (R - U ) = FF^ and F F - !1, Matrix F 
is a matrix of factor loadings and may be interpreted as the loading or 
relation of an attribute a on a generated set of factor components or di- 
mensions s. Multiplying the transformed factor loadings matrix by F would 
generate £2 which is the diagonal matrix containing the eigenvalues associ- 
ated with each factor. The above expresses that portion of the total 
common variance accounted for by each of the underlying dimensions of var- 
iation. The extracted dimensions in the orthogonal factor model are not 
correlated;, so each dimension expresses an additive share of the original 
variable set. In addition, to computation of the factor matrix, a set of 
factor scores may be generated. The factor scores matrix n S s is computed 
by the following: 

S = Z F S2 _1 

The factor scores provide the weighting of the observations on the factor 
dimensions. The spatial distribution and areal variations of these scores 
may then be mapped or used to generate multivariate regions. 

The preceding approach reflects an R mode factor analytic model. It 
is possible, however, to carry out the analysis for the rows of matrix A 
in which case a Q-mode analysis applies. The Q-mode analysis defines an 
immediate regionalization scheme in which the generated dimensions are 
composites of areal units. The majority of regionalizations employing 
factor analysis utilize the R-mode analysis and then attempt to cluster 
the generated factor scores. This approach permits a stronger definition 
of the variable structure and the rationale that would provide the base 
for a regionalization effort. 

2. Dimensional Analysis 

Dimensional analysis identifies the degree of similarity for each 
pair of observations by utilizing the concept of a taxonomic distance be- 
tween observations. The concept of an n-dimensional space is illustrated 
in the following diagram: 
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The diagram illustrates the taxonomic space for a three-dimensional area. 
The observations a, 6 , and X which are located in that space relative to 
their positions on the orthogonal vectors Vj_ 5 V 2 , and V 3 . Distance be- 
tween points in n-dimensional space follows from the utilization of the 
Pythagorean theorm: 


Where : 


x i 6 Vi 


~ Yi) ; 


is the distance (similarity) between points 
1,2,3,. .....n Dimensions 

are the values of characteristics for 
observation i 


The object of classification is simply to place in one group observations 
that are closer together in n dimensional space and to separate groups 
which are farther apart. The following two-dimensional diagram represents 
the above separation concept : 
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Berry [407] suggests that if the straight line distances are generated 
between observations it would be possible to compute a similarity matrix 
This similarity matrix would then be used as an input to a cluster anal- 
ysis. 

3. Grouping or Cluster Analysis 

Grouping or cluster analysis would take the preceding pairwise sim- 
ilarity matrix and through a stepwise procedure provide an optimal class 
ification of observations. There are several grouping algorithms which 
may be employed to obtain an optimal classification: 

- Centroid Groupings - This procedure groups with the objec- 
tive of minimizing squared distances between groups. The 
centroid grouping algorithm is the simplest and proceeds 
by identifying that pair of observations in which d|.= is at 
a minimum. It then takes the distance matrix (D) ana com- 
bines the row and column vectors of D into a single row and 
column vector representing the new group. The elements of 
these new vectors are the squared distances from the group 
centroid to all other points. The preceding is an iter- 
ative process with each successive stepwise procedure re- 
ducing the elements of the distance matrix until all linkages 
are defined. The following linkage tree demonstrates a cen- 
troid grouping technique for count: es: 

SPECIFIC 


GENERAL 

It may be observed that at each step those minimal dis- 
tances representing the greatest taxonomic similarity are 
linked together. Thus, the linkage tree represents an 
orientation from the specific to the generalized. Through 
the breakdown process we gain generality at the cost of 
individual detail. 

In the above example,; the determination of the breaking 
points would represent levels of a regional hierarchy. 

Most researchers define the optimal breaks in the branch- 
ing tree by noting where substantial jumps occur in the 
ratio of within to between group distances. 


Counties 



- Total Increment Grouping - This approach minimizes the in- 
crement to within group distances. At every step in the 
linkage analysis a check is made of total within-group 
distances so that each step minimizes the increment to 
within group distance. This approach is generally superior 
to the centroid algorithm for it provides a better control 
for the size of the final groups. It does require, however, 
a greater input of computer time. 

- Gravity Grouping Algorithm - This approach starts in the 
opposite direction from the preceding approaches. It begins 
at the highest general level and proceeds to the specific. 

The method permits larger groups to annex smaller groups 
over longer distances. The algorithm employs the basic 
gravity formulation: 

p 3 / 

Where : 

P£ = Measure of Mass for i 

P. = Measure of Mass for j 
3 

2 

d^j - Distance between i and j 

All of the above methods are subjected to the same problems of selection 
of significant cut-points to define the regional hierarchical levels. 

It may be noted that a minor modification should be introduced into 
these algorithms. The methods just discussed, did not incorporate a 
linkage structure having a contiguity constraint. These approaches simply 
linked areal units by similarity, the end result is the derivation of 
clusters of areal units of high similarity. A basic construct of a re- 
gion, however, is that the areal units assigned to it must be spatially 
contiguous. The linkage contiguity constraint may easily be introduced 
through the use of a binary connection matrix. The elements of the matrix 
that have a zero would indicate non-contiguous areas, while a one would 
imply that the places are joined over space. With the contiguity con- 
straint a linkage is made only with clusters that are linked in space. 

When clustering proceeds under the contiguity restraint, each clustering 
step reduces the number of remaining parcels by one. 

The preceding contiguity concept introduces an. important geographic 
component to the problem of regionalization. Haggett [392] notes that 
a large share of the regionalization literature is not truly geographic 
in that it merely classifies geographical data not locations. The con- 
cept of contiguity forces the researcher to carry on a "progressive ' 
comparison" of each unit with its neighbors. Haggett defines two major 
approaches that have been reasonably successful in accomplishing the 
above : 
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Variance Analysis - Drawing from the work of Zobler [409 3 re- 
gions may be constructed from smaller units by defining their 
components of variance. Zobler suggests that there is vari- 
ation among areal units within a region (within-region vari- 
ation) and that there is variation among the regions (between- 
region variation). Zobler uses variance analysis to assign 
West Virginia to three major regions with the following results: 

Variance Ratio 


Variance (V^) 

Variat^ I on ( V w ) 

V V w 

Alternative Assignments : 






To Mid Atlantic 

46.09 


8.91 


5.17 

To South Atlantic 

71.55 


4.66 


15.35 

To East- South Central 

72.13 


4.57 


15.78 


In the table the between -regional variance demonstrates the var- 
iation of regions around the grand mean for all regions. It may 
be noted that substantial differences are observed between the 
Mid- Atlantic Group and the two southern regions. The within- 
regional variance, on the other hand, shows the variations of 
the states around their respective regional means . The F-ratio 
generated from the preceding reflects the interregional dif- 
ferential, or in other words, how successful the grouping pro- 
cedure has functioned. In this case, the optimum allocation for 
West Virgina would be in the East-South Central region. The 
high variance ratio indicates maximal between group variation 
compared to within regional variability. 

Correlation Analysis - A Second approach is represented by a 
correlation analysis which denotes a measure of spatial con- 
tiguity by defining the strength of functional ties between 
areal units. Hagood [410,411] utilises this approach in di- 
viding the United States into contiguous groups of states with 
the objective that states in each group show a high degree of 
homogeneity in relation to agricultural and population attrib- 
utes. Based on this data agricultural and population profiles 
were generated and then standardized. Correlations were then 
computed between the profiles of adjacent states. The cor- 
relations exhibited the nature of "regional bonds" and were 
portrayed by lines of varying width on a map: 


BOND VALUES 
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Hagood later used the correlation linkages in a factor analyt- 
ic model to generate a single regional index. The preceding 
map indicates those states that are strongly linked as opposed 
to those that exhibit relatively weak bonds . 

In summary, the problem of regionalization is a basic issue in a 
cost-benefit analysis. Several methods of approaching the problem from 
both an identification and an assignment orientation have been illustrated. 
The choice of methods is dependent on the purpose, type of application, 
size and configuration of the areal base units , and the attributes or 
variables to be employed in the regionalization. When the above are uti- 
lized in conjunction with the earlier discussed methodologies the struc- 
ture of the spatial dimensions of a cost-benefit analysis may be derived. 


CONCLUSIONS 


The function of this paper has been to present a compendium of 
methodologies that measure the spatial component of events. In general, 
the paper has been devoid of direct applications. In this final section, 
the interactions of the previously mentioned methodologies will be placed 
in the context of a cost -benefit application for an early warning satel- 
lite system. It will not be the purpose of the appendix to actually 
measure the spatial components of the above technology. Its function 
is directed toward demonstrating which methodologies should be utilized 
and at what stage they should be introduced into the analysis. 

The utilization of a spatial/temporal orientation to events will 
provide the infrastructure for the methodological review. In an earlier 
paper Muraco [412] has shown the events leading to and following a nat- 
ural disaster may be placed in the context of a temporal frame over 
space. Multiple events characterized each stage of the pre and post 
disaster experience continuum. At each stage in the sequence various 
spatial components appeared. The accurate measurement of these compo- 
nents is critical both in the evaluation of the disaster itself and in 
the assessment of technology that may alter the event sequences. 

Pre -Event Phase - In the pre-event stages the major emphasis was placed 
on the assessment of cultural and physical adoptions, primary community 
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functions, risk and hazard perception and behavior, organizational tv ‘-ruc- 
ture, and the general characteristics of the existing communications 
system. The interrelationship of these various functional elements as 
noted by Muraco [412] results in a utility level that is specific in 
space and time. A change in the utility level, introduced by the adop- 
tion of new technology would imply the generation of new costs and 
benefits to the system. It is essential that the spatial components 
of these costs and benefits be defined. 

In the pre-event stages the cultural and physical setting must be 
isolated. In a spatial context this would require a definition of the 
settlement pattern and the nature of functional relationships between 
settlements. The methodologies that would be employed in this phase 
of the analysis would be associated with nodal dispersion and network 
analysis. The intensity of the settlement pattern and its degree of 
connectivity are important components of the regional infrastructure. 

The integration of the preceding with the functional linkages would 
establish the nature of the settlement hierarchy and the general geo- 
graphic orientation of the area. This establishes the physical-func- 
tional setting for the population. 

An essential aspect of the pre-event sequence phase is associated 
with the identification of risks and the dissemination of emergency 
warnings. This phase of the analysis would utilize the previous method- 
ologies coupled with attempts at multidimensional regionalization. The 
identification and classification of cultural areas would employ a re- 
gional assignment methodology. It is through this approach that the 
phenotypic composition and distributional character for the area may be 
differentiated. Based on the spatial distributions of the phenotypes, 
new behavioral regions may be defined by using a regional identification 
methodology. These regions would then he employed in isolating per- 
ceived risk and non-risk areas , as well as , predicting the phenotypic 
implications of the perceptions. The following diagrammatically illus- 
trates the preceding: 
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The Event Phase - The event stage would also utilize a composite of the 
previously discussed methodologies. During the event, problems of cir- 
culation and the identification of major impact regions would be central 
spatial considerations. The identification of potential bottlenecks and 
barriers ^ would necessitate a network and flow analysis. Problems of 
regionalization would also be central issues during this stage. The 
identification of the location and extent of the event over space are 
essential elements in measuring the primary costs and the anticipated 
secondary costs. Much of the information generated during the event 
phase must be coupled with the earlier distributional analysis of the 
regional infrastructure. If the event, has a serious impact on that in- 
frastructure, its effect on total utility will have substantial conse- 
quences in both time and space: 


Communi 



Emergency Functions 

NETWORK\AND DISTRIBUTIONAL 
ANALYSIS 


Phenotypes 


Administrative Functions 


Post-Event Phase - The post-event stage would also require the inclusion 
of spatial components. During the early portions of this phase regional 
assessments of damage and emergency evacuation would be major activities. 
The derivation of damage regions to concentrate efforts for recovery 
would be illustrative of a multivariate regional identification problem. 
The interface of the damage regions against the cultural regions defined 
in the pre-event stages may represent a means of establishing the dis- 
tributional aspects of the disaster on various phenotypic groups. In the 
latter phases of this stage, adjustments of a normative nature may take 
place in coping with the reoccurrence of the event. The development of 
optimal interaction networks and administrative regions would be critical 
considerations in assessing the existing technologies efficiency: 
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The evaluation of the effectiveness of emergency and administrative func- 
tions requires the identification of optimal hinterlands that would min- 
imize harrier effects. The organizational delivery system would be 
assessed, in part, by its propensity to interact with the impact regions. 
To identify these interactions requires a comprehensive analysis of the 
functional networks over space. 

The preceding temporal sequences illustrated the methodologies re- 
viewed in this paper within the context of evaluating an event such as 
a natural disaster. In a technology assessment problem the changes that 
the technology introduces to the sequences over space would become the 
central issue. In the event and post -event phases, for example, sub- 
stantial changes in geographic orientation and interaction may be intro- 
duced by an early warning communications technology. Knowledge of the 
extent, and intensity of these changes over space would be required in 
assessing the true benefits and costs derived from the innovation. 


SUMMARY 


In the following table the general application of the methodologies 
discussed in this paper are displayed: 
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The table provides a comparison of the different methodologies and their 
potential applications. The circles in the body of the table suggests 
under which applications a methodology is suitable. In all of the sample 
applications multiple methodologies would be required to insure adequate 
spatial coverage. 

In summary, it is the purpose of this ppendix to provide a compen- 
dium of geographic methodologies that would have applicability in a cost- 
benefit problem. The ppendix is largely educational in structure and 
designed to provide the reader with a background in the techriques 
geographers employ in the analysis of events over space. The operax’onl- 
ization of these methodologies in a cost-benefit problem is demonstrated 
in an earlier section of this report. 
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This Appendix deals with multidimensional scaling. Multidimensional 
scaling (MDS) is a generic term that identifies several techniques for 
obtaining metric, cardinal information from a matrix of ordinal-level 
rankings. It has been described as a response to "the problem of repre- 
senting n objects geometrically by n points, so that the interpoint 
distances correspond in some sense to experimental dissimilarities between 
objects" [345] . Thus, MDS is concerned with obtaining a satisfactory 
spatial representation of relationships among stimuli or objects. Usually 
no tests of significance are computed, although such statistics as STRESS 
or SQUARIANCE are used to uncover the appropriate final configuration, 

The points are clustered in multidimensional space so that axes may be 
inserted. Measurements made along them can describe The relationship of 
specific points to the several latent dimensions [346], 

The purpose of this Appendix is to present the basics of the MDS 
methodology. Consideration is given to the basic principles of MDS, to 
the common properties of MDS computer algorithms, and to the methodological 
limitations of the methods. 


PRINCIPLES OF MULTIDIMENSIONAL SCALING 


One reason that has been advanced for the lack of empirical theory in 
the behavioral scienes is that they lack the capability for precise meas- 
urement. The importance of this assertion has diminished with the develop- 
ment of many non-parametric tests and measures during the past twenty years. 

Simultaneously with these advances in ordinal non-parametric statistics, 
the behavioral sciences have benefited from developments in the field of 
cumulative scaling. Table 1 divides attitudinal, behavioral and cognitive 
scales into direct and indirect types. In direct scaling, respondents are 
asked to "rank themselves" (ordinal scales) or to locate themselves along 
a continuum (some interval scale procedures). Ej eriments also have been 
made to establish self-rating at the ratio level of measurement. For a 
variety of reasons, direct ratings are problematic. Preferred alterna- 
tive scaling approaches involve some third person to do the rating. 'What 
this means is that it is deceptively easy to claim that a reliable rank- 
ordering (or a reliable cardinal measurement) has been reached. Actually, 
even ordinal measurements are more difficult to attain than it might 
appear. 

Cumulative scales can check upon the truth of ordinality, but they 
also can be used to generate ordinality from simple nominal information. 

A cumulative scale assigns magnitudes to entities in relation to the per- 
ceived amount of a common attribute that each possesses. That is, it 
assesses behavior by combining several similar variables into a single 
composite indicator. It assigns scale scores along a single continuum or 
"dimension", and these scores constitute an ordinal or higher level of 
measurement. A wide variety of scales have been constructed, and measure- 
ments have been attempted of preferences, of hostility, job prestige. 
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psychomotor dexterity, friendship, and of untold other concepts. At the 
ordinal level, the work of Louis Guttman has been very important; Likert 
and Thurstone scales are representative of work toward interval levels 
of measurement. Psychometricians have been laboring during the past few 
years at attitude measurement at even the ratio level — at developing 
measures with a true zero point. 

Table 2 depicts a typical Guttman-type scalogram situation. Simple 
"yes” and "no" responses (nominal dichotomies) have been entered in 
accordance with the types of storm damage suffered by the residents of 
forty-four counties. The categories of loss have been ordered according 
to their degree of "hardness" as represented by the four marginal totals. 
The number of counties experiencing each level of loss is specified under 
N, and the scale score is indicated at the right. Of course, as with all 
ordinal measures, the interval between scale scores 3 and 4 may bear no 
relationship to the one between 2 and 3. Only an interval (or ratio) 
level of measurement is concerned with inter-class magnitudes. 

Multidimensional scaling can be understood best in relation to uni- 
dimensional scale analysis. In the preceding example , the use of model- 
data for illustrative purposes allows statistical tests of unidimension- 
ality. The units that are scaled may be almost anything — survey re- 
spondents, governmental units, verbal concepts — and when found to be 
unidimensional they can be ranked in relation to the degree that -they ' : 
are perceived to possess an attribute, X. ' ' 
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TABLE 2 


MODEL GUTTMAN-TYPE SCALOGRAM 


Counties 
Reporting 
(N = 44) 

Types of Damage Reported 

Scale 

Scores 

Rain 

Damage | 

| Housing 
[ Losses 

Disease 

Outbreaks 

Deaths 

15 

No 

No 

No 

No 

1 

12 

Yes 

No 

No 

No 

2 

10 

Yes 

Yes 

No 

No 

3 

5 

Yes 

Yes 

Yes 

No 

4 

2 

Yes 

Yes 

Yes 

Yes 

5 


MULTIDIMENSIONAL SCALING 


Multidimensional scaling assumes that objects to be scaled may possess 
more than a single common attribute, or that Attribute A can be broken down 
into several sub-attributes. Multidimensional analysis transforms the 
threat of non-uni dimensionality that plagues Guttman-type scales from a 
debt into an asset. Using multidimensional scaling, for example, it can 
be demonstrated that color is a composite of three attributes: hue, 

brightness and saturation. Multidimensional scaling is concerned, then, 
with organizing the 0* objects so that they can be arranged along two di- 
mensions when a pair of £cales exist among data, along three dimensions 
when three underlying scales are present, or up to (n - 1) dimensionr „iien 
n_ objects are to be scaled multidimensionally . Figure 1-1 assumes the 
existence of three dimensions among a set of hypothetical data bearing 
upon the "Agnes" hurricane disaster of 1972. It should be pointed out : 
that not only may solutions be sought in even higher dimensionalities, but 
that they may be found id non-Euclidean space as well. 

Figure 1-1 shows that objects can be arranged in three-dimensional 
space. The "objects" indicated here might be units (such as persons who 
have responded to a series of survey questions) or they could be stimuli 
(the survey questions themselves). Again, they might be units or stimuli 
acquired from one of the many non-survey data gathering methods. The 
model depicted here, identifies the first dimension as financial loss in- 
curred from the Hurricane Agnes disaster. Those who suffered more heavily 
in a purely financial way would be located very close to this axis, and 
Would be scaled high on the concept. Those who actually benefited finan- 


cially from "Agnes" — from selling supplies or from profitting on new con- 
struction to replace damaged structures — also would be on this dimension 
but would approaching the value of -1. If all of the n_ objects were 
scaled unidimens ionally , such tests as coefficient of reproducability, 
coefficient of scalability and minimum marginal reproducability all would 
be low because it is apparent from the configuration in Figure 1-1 that 
more than a single dimension is present among the points. Dimension II, 
the vertical axis, might represent mental stress . The third axis could 
record the deliterous impact of the hurricane upon personal health . Be- 
cause a number of points remain at the centroid of the configuration, 
however, it may be that additional dimensions are required to extract ad- 
ditional explanations from the data. The expectation of higher dimension- 
ality also would come from a more general scattering of points within the 
cubic space or from any peculiar clusters that are extant within it. Un- 
acceptably high levels of multidimensional scalability (e.g. , STRESS, 
SQUARIANCE) also call for a search at higher dimensionality or in some 
other Minkowski-p_ geometric space. 

Although the roots of MDS extend into factor analytic concepts and 
include work in matrix analysis and psychophysics [355,356!! its independent 
development appears to owe more to the innovations of Torgerson and of 
Coombs during the 1950s. The classic formulation was detailed in two 1962 
papers by Shepard, however, which showed that tightly constrained geo- 
metric plots of points can be constructed from only rank-ordered informa- 
tion about relationships among a set of n objects. Shepard stressed the 
need for a monotone relationship between experimental dissimilarities 
(represented conventionally by the lower-case letter delta — 6 ) and the 
distances among points (or d.) in the configuration scattergram shown in 
Figure 1-4 ( c ) . Shepard demonstrated that a matrix expressing ordinal re- 
lationships among objects could be used to group the objects themselves, 
and that interval-scale axes could be inserted among the points so as to 
yield metric outcomes from only ordered inputs. 

The first few years after Shepard's paper saw a number of advances 
in multidimensional scaling. In 1964, Kruskal [345,347] improved upon 
Shepard's approach by using a measure of "badness of fit" to assist in 
coverging the distances among the points with the given matrix of dis- 
similarity orders. In Kruskal* s approach, a monotonic regression of 
distance is preformed on dissimilarity, and the normalized residual 
variance — STRESS — is used as a measure of how poorly the point config- 
urations reflect the rankings among the data. Other approaches to MDS 
followed in short order: Guttman-Lingoes ' [348] smallest space analysis 

(SSA); Torgerson and Young's TORSCA series of computer algorithms, Carroll 
and Chang's INDSCAL [357] and PROFIT [358] routines, and a number of 
others. Although these conceptualisations vary in approach, in numerical 
operations, in size capabilities and even in their ultimate goal, a 
common model and explanation of the operation of most MDS computing pro- 
cedures can be offered. 

Although the objective of all MDS routines is the spatial represen- 
tation in ratio-scale of simple ordered relationships among n_ objects, 
or the description of ordinal data at a cardinal level, quite a large 
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number of computing algorithms have been developed to accomplish this 
task. They vary with regard to (1) the types of associational matrices 
that may be input, (2) methods of identifying an initial configuration 
from which the ordinal scaling may begin, (3) types of statistical guides 
or tests that are used in reaching a solution, (4) methods of "jiggling" 
the points during the scaling process, (5) types of Minkowski-^ spaces 
in which findings can be reported, and (6) computing efficiency. 

All MDS programs include the steps identified in Figure 1-2. They 
begin from a set of input data that represent similarities or differences 
in attitudes, behaviors or cognitions. Sometimes extant as an integral 
part of the algorithm, sometimes supplied externally for it, often there 
is a routine for pre-processing the data that allows the special treatment 
of tied or missing cases or for reducing the data from a cardinal to a 
continuous-ordinal level. Third, some initial configuration of data points 
is proposed through either randomization or a purposeful assignment process 
so that the actual clustering might begin from this point. Fourth, the 
difference between the rank-order of interpoint distances and the extant 
point configuration is ascertained through computing a "loss function" . 
Whenever the configuration is found to reproduce acceptably the ordinal 
structure of the data, it could be decided to terminate further iterative 
"jigglings" of points and to output the present result. Failure to have 
passed the termination test or to have exhausted the number of iterations 
called for in the solution (e.g. , 25 iterations) again activates the con- 
vergence method for carefully moving the configuration of points so that 
they match with the given input matrix of rankings as nearly as possible. 
Consider each of these steps more closely. 

1. Input Data . 

Any type of raw data that satisfy the requirements of a data matrix 
concept and that can be measured at the ordinal level or higher can be 
used for input to any MDS algorithm. Many of the early applications of 
scaling were drawn from the fields of psychology and marketing. Perhaps 
because the extensive use of "pair comparison" techniques similar to 
those discussed in operationalizing T.A.P. as a basis of constructing 
rank -orderings has been common to these areas, many publications involving 
MDS use this approach. In addition, the focus on identifying dissimilar- 
ities among brands of automobiles, foods or colors might lead one to 
believe that the potential uses of MDS are limited. As Kruskal [345] 
pointed out a decade ago, however, the techniques can work with measures 
of similarity or dissimilarity, with correlation/association coefficients, 
and with information of a "most diverse kind". These data might be 
gathered in any of a number of ways. They might be "discrete" or "aggre- 
gate" sources; they may be exante or post-hoc indicators of attitude, 
behavior or cognition. Although there are conditions under which the 
data themselves can be analyzed directly, it is customary to correlate 
the columns of a raw data matrix (i.e., the stimuli or variables) with 
one another to obtain a matrix of stimuli (the R-technique of factor 
analysis). Where interest centers upon the actions of the persons or 
other units of analysis in the data matrix, however, a resulting matrix 
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of units can be constructed by associating the rows with one another (the 
Q-technique of factor analysis). Six types of resulting or secondary 
matrices ordinarily are encountered at inputs to MDS algorithms. They 
are shown in Figure 1-3. 

The most common input matrix is the symmetric triangular one shown 
in Figure 1-3 (a). In this case, raw data have been associated into n(n-l)/2 
unique combinations. Because the main diagonal records the association 
of variables with themselves and the reflexive ci efficient always is unity 
(1.0), it is usually omitted from further analysis. Symmetric triangular 
matrices result from correlating all possible pairs of stimuli with one 
another, or from associating all pairs of units with some sort of co- 
efficient of agreement. 

Figure 1-3 (b) is a non-symmetric matrix with the main diagonal again 
being missing. In this case, one or more of the entries in the lower 
left-hand segment of the matrix will be dissimilar from its analog in the 
upper right corner. If the coefficients entered in matrices such as these 
are indicated as r^j, then r^j ^ Eji* r s is a very realistic type of 
data representation, for it represents the myriad situations in which in- 
teractions are not perfectly reciprocal. One individual may think more 
highly of a second person in a social setting than the latter does of the 
former. The availability of a scarce raw material will have a stronger 
impact upon the behavior or a manufacturer in a purely competitive market 
than the activity of that manufacturer will have upon the raw material 
supply. The input-output models of Leonte if exemplify this sort of asym- 
metry, and the development of parametric and non-parametric coefficients 
of similarity or dissimilarity (e.g. , asymmetric uncertain coefficients 
from information theory; lambda^ ; Sorrier’ s^ x and Somer 1 s X( j ) enable these 
relationships to be summarized neatly. 

Figure 1-3 (c) depicts an asymmetric matrix in which the main diagonal 
does not contain unities . Such a square matrix can result when the con- 
text in which the stimulus is encountered varies, or because of some 
inability of respondents to identify a concept as being identical with 
itself. Again, the substitution of communality estimates in the main 
diagonal for principal components factoring exemplifies the use of a 
square matrix. 

Other types of input matrices are possible as well. Conditional 
matrices can be described by comparing the row or columnar entries (Fig- 
ures 1-3 (d) and 1-3 (e) — something that yields a rectangular effect. Again, 
Green and Carmone [349] have described partitioned matrices in which 
unconditional, symmetric relationships are obtained for a core set of 
units or stimuli, and the other interrelationships are obtained on the 
basis of their interaction with the respective core set entries. 

2 . Preprocessing . 

When data about to be scaled multidimensionally are expressed in 
interval or some other form than that of a strict rank-ordering, they 
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must be converted into an ordered matrix before the actual scaling can 
begin. During the early development of MDS this sometimes was done by 
an external preprocessing routine. One of these, TRICON [350], construct- 
ed an ordinal similarity matrix of up to only fifteen stimuli. Another, 
WAGS [351], was designed to accomplish this with data encoded into zero- 
one dichotomies. The use of such small quantities of objects for scaling 
would impose serious limitations on research. Fortunately, such MDS 
packages as M-D-SCAL, TORSCA and SSA allow the input of as many as eighty 
stimuli and do the necessary rank -ordering automatically. 

The procedures for preparing non-ordered data into strictly ordered 
matrices have a common goal, although they do vary with regard to their 
treatment of intransitives, of missing or tied data, and of asymmetric 
information. Table 3 provides an example of how this ranking process is 
done. A matrix of indexes of association is supplied in this case. These 
happen to be phi coefficients, or coefficients that approximate Pearson's 
r when x and y are dichotomies. Matrix (b) shows that the information 
has been reduced to a ranking. Of course, there might be little or great 
similarity of appearance between the coefficients with which we begin and 
the ordered matrix with which we conclude. The degree of similarity will 
depend upon both the data with which we work and the coefficients used to 
summarize them. In this example, the median phi value among these 45 co- 
efficients is .09. Were these initial coefficients expressed as Yule’s 
Q’s, however, the twenty-third (or median) value would have been .25, 
given the same data configuration. 


TABLE 3 

CONVERSION TO RANK-ORDERINGS 


(a) Phi Coefficients 
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3. Initial Configuration. 


Spurious outcomes can result from MDS techniques under certain con- 
ditions. One problem is that of the degenerate solution — one in which 
the goodness of fit between the scaled points is unacceptable but where 
STRESS approaches zero. It can occur when the object points cumulate in 
clusters so that within-cluster dissimilarities are much less than those 
between clusters, or when all individuals polled have chosen one stimulus 
in preference to all other possibilities. Therefore, quite different 
final solutions are obtained when the same data are re-analyzed. (This 
is not the same as the local minima difficulty. There, the step value 
that regulates the distance that points are to be moved on the next 
iteration is too small. This means that points can be trapped in loca- 
tions that are not optimal solutions. The appropriate remedy in this 
case is to increase the step values). 

It has been found that the choice of an initial configuration of 
points that is as near to the likely final MDS solution as possible mini- 
mizes the possibility of degeneracy. If, when the "jiggling" of points 
begins, a correct solution is close at hand, the likelihood of a degener- 
ate result is greatly reduced. MDS algorithms exhibit a variety of 
techniques for selecting initial configurations. Although points might 
be input randomly, most users today prefer some type of non-random ini- 
tial configuration. M-D-SCAL allows an input configuration of one's 
choosing, such as a square, circle, or other shape. Metric inputs also 
are used, as with the TORSCA series. Here, a semi-metric configuration 
is prepared by converting the data to scalar products that are factor 
analyzed to yield a preselected set of r dimensions. In TORSCA-9, five 
cycles of this semi-metric iterative process normally are repeated prior 
to the actual non-metric analysis. Often the initial configuration is 
very close to the final outcome. 

4. Converging the Data and the Object Points . 

Given an initial configuration of points, the fundamental problem 
of MDS is to move them in such a way that the distances between them best 
reproduces the order of entries in the strictly-ordered matrix. Early in 
the development of multidimensional scaling, solutions were sought by 
placing pins in a board, measuring the distances among them, comparing 
the rankings of these distances with the ordered matrix, moving the pins 
some more, and repeating these steps until the fit between pins and data 
were compared with the. ranks of the inter-pin distances, and not with 
the actual lengths that separated them. 

Figure 1-4 shows the interrelationship among ranked data, a config- 
uration of points, and a scattergram that is used for calculating a 
measure of STRESS useable in re-configuring the scaled objects. Let the 
interpoint distances between objects i and j be d^ j , and the corresponding 
ordering in the matrix be • ■ . 


FIGURE 1-4 


(A) DIAGRAM, D 


DATA, POINTS AND STRESS DIAGRAM 

(B) MATRICES, DEEXA 



CASE X. PERFECT MATCH OF 
MATRIX AND DIAGRAM 

BCD 
A X 2 5 

B 3 6 

C 4 

(STRESS *= 0) 


CASE II. IMPERFECT MATCH 

BCD 
A 1 3 5 
B 2 6 


c h 

(STRESS > 0) 







Kruskal's description of a procedure for mapping points is typical 
in many ways of convergence methods. Given an initial input matrix of 
ordered relationships. Figure 1-4 (b), an initial clustering of objects can 
be done (Figure 1-4 (a). Next, a scattergram can be constructed, 1-4 (c) 
so that the data rankings are represented on the ordinate axis and the 
"measured" distances on the abscissa. What amounts to a best-fitting 
line then is fitted among the d^j and the distance between the actual 
objects and their expected location (denoted as "hat" or dij) is 
found. This is the distance between where the points presently are lo- 
cated and where they would be expected to be if there were a perfect 
monotonic relationship between points d^j and ranking <5^4 . This differ- 
ence is squared and summed over all cases. Obviously, when the scaled 
object is located where it "should" be — where a perfect ^monotone rela- 
tionship exists between data and point, the difference d^j - d^j = 0. 

An equation 


I(di j - <Hj ) 2 

can be interpreted as the "residual sum of squares" for the real and ex- 
pected interpoint distances. 

Kruskal's important contribution to MDS is his development of a "bad- 
ness of fit" statistic that can be used both to judge the adequacy of an 
intermediate or final solution and to make the minute adjustments or 
"jiggles" of the points in the iterative process that ensues. The statis- 
tic is a normalization of the residual sum of squares, and is called STRESS: 


STRESS 



When all of the clustered points fall directly along the monotonic best- 
fitting line, STRESS will be zero and a perfect fit will exist between 
points and data. While STRESS always will equal zero when the number of 
dimensions in which a solution is sought is one less than the number of 
points according to Bennett and Hays [352 ] , it also has been argued by 
Guttman [348] that perfect configurations will obtain in n-2 dimensions. 
Less satisfactory values usually are encountered at lower dimensionalities . 
Kruskal offered this evaluation of STRESS and of the corresponding good- 
ness of fit: 
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Values of STRESS and similar indices according to Klahr [359] are larger 
with greater numbers of scaled objects, with greater error in the data, 
and with lower-dimensional solutions. 

The values of d^j and d— that are found in computing STRESS also 
are used to guide the J direction and magnitude of point jigglings in re- 
ducing the discontinuities between points and input data. Kruskal notes 
[345] the process of accomplishing this can be understood intuitively as 
a process of successive approximation, theoretically as in iterative tech 
nique of numerical analysis called the method of steepest descent, or 
operationally as any of the applications designed to accomplish the jig- 
gling process. Alternative computational approaches also exist. One of 
them is Gleanson's alternating algorithm. This process switches between 
satisfying metric distance requirements first, then by satisfying the 
ordering requirement. A solution is sought whereby both interval and 
ordinal needs are fulfilled. 

It should be noticed in Figure 1-4 (c) that the calculations in the 
scattergram are executed on the distances (d^j) arranged on the abscissa, 
rather than on the 6^ • shown on the ordinate. To base the calculation in 
reverse would imply that the distances between the ordinal data were of 
equal intervals — that the data were interval in level — precisely the 
assumption about cardinality that MDS was devised to circumvent. As 
Kruskal [345] points out deviations between the curve and the plotted 
points are measured along the distance axis because 

". . .if we measure them along the dissimilarity axis, we shall 
find ourselves doing arithmetic with dissimilarities. This we 
must not do, because we are committed to using only the rank 
ordering of the dissimilarities .... We wish to measure good- 






ness of fit in such a way that monotone distortion of the dis- 
similarity axis will not have any effect. This clearly prevei. ;s 

us from measuring deviations along the dissimilarity axis" [345]. 

5. Dimens ionality . 

In most scaling situations involving more than a few scalable ob- 
jects, acceptable convergence of points and data can be attained only 
multidimensionally. Examples used above suggest correctly that acceptable 
two-dimensional solutions — and even outcomes in a single dimension — are 
possible. The development of computerized MDS techniques since 1964, how- 
ever, have made it possible to attain not only quick and accurate solutions, 
but also to reach configurations of points in more than two dimensions. 
Between the extremes of one or two-dimensional solutions on the one hand 
and n-1 or n-2 solutions on the other, the question of identifying an 
appropriate dimensionality for obtaining solutions is important. 

Identification of an optimal number of dimensions to be used in a 
solution can be made in at least two ways. First, the values of the 
"badness of fit" statistic (STRESS in M-D-SCAL, SQUARIANCE in TORSCA-9, 
COEFFICIENT OF ALIENATION in SSA, etc.) can be plotted for each dimension- 
ality of a given body of data. Figure 1-5 shows that a distinct bend 
appears at the fourth dimension for the Euclidean solution (rho = 2.0), 
from which it can be inferred that this would be an appropriate number of 
dimensions. In the case of rho = 1.0, the bend is found at three dimen- 
sions. Sometimes there may be no distinct "elbow" and the investigator 
is provided with no firm clue as to the appropriate dimensional .y . A 
second answer to the question of how many dimensions should be extracted 
is to obtain only as many as can be interpreted reasonably in a substan- 
tive sense. 

6. Solution Spaces. 

Although the implication to this point often has been that objects 
are being scaled in Euclidean multidimensional space, this need not be 
the case. STRESS, SQUARIANCE, and other indicators of badness of fit 
can operate in a wide variety of outcome spaces. Euclidean space is but 
a special case of L -Norms, or of Minkowski rho - metrics . Under more gen- 
eralized rules of spatial geometry, a Minkowski-^ = 2.0 yields the Euclidean 
solution. When rho = 1.0, however, a different spatial outcome ensues. 

This is the so-called City Bloc or Manhattan metric space. City bloc 
solutions are useful in mapping the real distances that persons must travel 
under certain conditions or the psychological space perceptions that they 
have internalized. Consider as an example that the c .'stance between two 
urban locations may not be thought of "as the crow flies", but rather is 
conceived of as "up to the comer, then turn to the right". Minkowski 
metrics range from zero through infinity. Theoretically, any Minkowski-p_ 
spaced can be treated multidimensionally. Translation of precisely what 
is meant by a rho- space of 0.8 or of 3.5 can be troublesome, however. 

Figure 1-5 also suggests that measures of stress can be plotted across 
varying Minkowski spaces to aid in the search for an appropriate dimension- 
ality. 
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7. Outputs . 


The results obtained from any MBS algorithm vary somewhat, but typi- 
cally they include a listing of the matrix of similarities or dissimilar- 
ities, a table of the cardinal coordinates or configurations of points 
with respect to each of the dimensions called for in the solution, and a 
series of scatter diagrams cv plots for the visual display of the arrange- 
ments of points in the appropriate dimensionalities. Also included may 
be a scattergram of distances by original rank order, as shown in Figure 
1-4 (c). 


METHODOLOGICAL ISSUES 


Before addressing the topic of the use of MDS in total assessment 
profiling, two methodological issues must be addressed: (1) the reasons 

that have been advanced for using non -metric MDS, and (2) the limits or 
caveats that must be understood by those who would use it. 

Consider the uses of MDS. First, its use is justifiable when data 
are ordinal and it is inappropriate to assume that the intervals between 
the ordered units are equal. Oddly enough, even since the development of 
MDS the literature of many fields continues to apply metric techniques to 
ordinal information. Second, it has been pointed out by Lingoes and 
Guttman [348] that a smaller space is required to reflect ordered relation- 
ships alone than is needed to represent both order and intervening distance 
intervals. Third, reducing stronger data to an ordered level can make 
clearer the dimensions that can be identified uhrough MDS than would have 
been the case with factor analysis. This is because ordinal coefficients 
set aside the requirement of linearity that constrains techniques based 
upon least-squares measures. However, recent developments have been en- 
couraging in the area of polynominal factor analysis. Finally, as Lingoes 
and Guttman, also have observed, "when some lawful structure or pattern 
is present in the data, e.g., a simplex, a circumplex, or a radex, a non- 
metric analysis will reveal the configuration whereas a metric approach 
will obscure the lawfulness" [353], 

What can be said about the limitations of this family of approaches? 

It has been twelve years since the publication of Shepard’s major article. 
The first six years after 1962 were characterized by innovative develop- 
ments of several sorts in the field of MDS; since 1968, however, more 
attention has been devoted to consolidating those gains by testing the 
performance of the several approaches. The application of multidimensional 
scaling in Total Assessment Profiling must be approached with these recent 
evaluations in mind. 

The best general summary of methodological limitations in MDS prob- 
ably is that given by Green and Carmone [349 ] although a growing body of 
studies consider the relative merits of two or three algorithms at a time. 
These authors have grouped their methodological caveats and suggestions 


for further research into three areas: computational problems, empirical 

considerations, and conceptual difficulties. 

Computational problems are best approached by comparing the perfor- 
mance of several algorithms under similar conditions. Involved here are 
questions of uniqueness and determinateness of solutions within a single 
approach (e.g. , does M-D-SCAL successfully replicate the first solution 
on successive runs?) the invariance of solutions across algorithms, and 
the impact of error upon algorithm efficiency. It may not be at all sur- 
prising to learn that comparative studies of performance reveal few dif- 
ferences among the leading three approaches to MDS. One investigation 
reported by Green and Carmone [349] attempted to recover a synthetic con- 
figuration of points (the capital letter "R" represented by 27 points in 
two dimensions) and found that TORSCA-8, M-D-SCAL IV and SSA-1 performed 
rather similarly" under both error free and "noisy" data conditions. In 
a comparative study of the Shepard-Kruskal and the Torger-^on approaches 
to multidimensional scaling, it was found that, in thirty-six simulations, 
both "models produced highly accurate solutions. ..." £354]. Again, in 
a comparison of M-D-SCAL, SSA-1 and TORSCA-9, it was found by Spence C346] 
that "differences between the solutions obtained by the algorithms were 
typically so small as to be of little practical importance". 

Since the computational considerations raised by Green and Carmone 
seem difficult to resolve — fortunately, because the leading MDS procedures ' 
all perform quite well — what is the best course to follow in any one case? 

A conventional approach when selecting the best alternative from among a 
small number of choices probably is to do the analysis each way. Some- 
times one solution will be quite superior to the others , and a knowledge 
of the contrasting techniques allows the investigator to tell why these 
contradictions have been found. 

Methodological matters relevant to the use of MDS in total assessment 
profiling include, as well as computational considerations, certain em- 
pirical and conceptual considerations. Empirically, one must consider 
the impact of data gathering techniques upon scalability. Although the 
questions raised by Green and Carmone are important ones, they are common 
to many critiques of multivariate statistical methods and need not be 
considered at this point. Conceptually, there Is considerable food for 
thought. Care will have to be given to the ecological fallacy, or to the 
problem of "analysis at the wrong level". When TAP findings are reported 
on the basis of survey data, for example, it must be understood that what 
is true for the individual is not necessarily valid for the group. If 
the recent concern with ecological problems can be used as one example, 
we can see that. the benefits and costs that accrue at the level of the 
firm are not identical with those that obtain at a higher level of ag- 
gregation. Externalities seldom are foreseen by isolated respondents. 

Other conceptual issues are significant. Just how well do the coefficients 
of similarity or preference reflect social realities? How quickly, and 
why, do the preferences mapped at any given time alter dischronicly? How 
do they change with distance? These and other important theoretical ques- 
tions are addressed in the larger context of the Total Assessment Profile. 
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APPENDIX J 


A Tentative Register Of American Values 


Source: Baier, Kurt and N. Rescher (eds.), Values 

and The Future, The Free Press, New York, 
1969 
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SELF-ORIENTED VALUES 


1. Personal "material" welfare (the right to life and the pursuit of 
happiness ) 

a. health (physical and mental well-being) 

b. economic security and well-being ("materialism" and the 
American way of life) 

c. personal security (stability of the conditions of life) 

2. Self-respect (the right to be treated as a person and as a member 
in good standing of the community ; honor, h on or ab 1 en e s si 

3. Self-reliance (self-sufficiency, rugged individualism and the 
pioneer tradition) 

4. Personal liberty (the right to endeavor to "shape one's own life", 
to work out major facets of one's own destiny and to go one's own 
way) 

a. freedom (from interference) 

b . privacy 

c. property rights 

5. Self-Advancement ("success", ambition, diligence) 

6. Self-Fulfillment (and "the pursuit of happiness") 

7. Skill and Prowess 

a. the intellectual virtues (intelligence, education, know-how, 
realism, practicality, versatility, etc.) 

b. the physical virtues (strength, dexterity, endurance, good 
appearance, cleanliness, etc.) 

c. the virtues of the will (strengths .of character) 

1. readiness for hard work (industriousness) 

2. toughness (fortitude, endurance, bravery, courage) 

3. initiative and activism (the "go getter" approach) 

4. self-control (temperateness, sobriety) 

5. perseverance and stedfastness 
d„ competence (pride of workmanship) 

e. inventiveness and innovativeness 

f. initiative (the "self-starter") 

g. well- informedness (access to information, being "in the 
know" ) 

h. faith ("believing in something" including "having a sense 
of values") 

’i. appreciation and appreciativeness (of "the good things of 
life") 


II. GROUP-ORIENTED VALUES 


1. Respectability (group acceptances avoidance of reproach, good 
repute, conformity, the ’’done thing" and the "herd instinct") 

2. Rectitude and personal morality (honesty, fairness, probity, 
reliability, truthfulness, trustworthiness — the "man of honor") 

3. Reasonableness and rationality (objectivity) 

4. The domestic virtues (love, pride in family role, providence, 
simplicity, thrift, prudence, etc.) 

5. The civic virtues (involvement, good citizenship, law-abidance, 
civic pride — the "greatest little town" syndome) 

6. Conscientiousness 

a. devotion to family, duty 

b. personal responsibility and accountability 

c. devotion to principle (especially of one's religion — 

"the godfearing man") 

7. Friendship and friendliness 

a. friendship proper 

b. loyalty (to friends, associates) 

c. friendliness, kindliness, helpfulness, cooperativeness, and 
courteousness (the good scout; "getting along with people") 

d. fellow-feeling (compassion, sympathy, and "love of one's 
fellows") 

e. gregariousness 

f. receptivity (openness, patience, "the good listener") 

g. personal tolerance ("live, and let live", "getting along 
with people") 

h . patience 

8. Service (devotion to the well-being of others) 

9. Generosity (charity, openhandedness) 

10. Idealism (hopefulness in human solutions to human problems) 

11. Recognition (getting due public credit for the good points 
scored in the game of life; success and status) 

12. Forthrightness (frankness, openness, sincerity, genuineness; 
keeping things "above board", the fair deal) 

13. Fair play (the "good sport") 


III. Society-Oriented Values 


1. Social welfare (indeed "social consciousness" as such) 

2 . Equality 

a. tolerance 

b. "fair play" . fairness 

c. civil rights 

3. Justice (including legality, proper procedure, recourse) 

4. Liberty (the "open society"; the various "freedoms") 

5. Order (public order, "law and order") 

6. Opportunity ("land of opportunity" concept; the square deal for 
all) 

7. Charity (help for the "underdog") 

8. Progressivism optimism (faith in the society's ability to solve 
its problems) 

9. Pride in "our culture" and "our way of life" 


IV. NATION-ORIENTED VALUES 

1. The patriotic virtues (love of country, devotion to country, 
national pride) 

a. national freedom and independence 

b. national prosperity and national achievement generally 

c. patriotism and national pride 

d. concern for the national welfare 

e. loyalty (to country) 

f. chauvinism (nationalism, pride in national power and 
preeminence ) 

2. Democracy and "the American way" 

3. "Public service" in the sense of service of country (the nation) 


V. MANKIND-ORIENTED VALUES 


1. The "welfare of mankind" 
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a. peace 

b. material achievement and progress 

c. cultural and intellectual achievement and progress 

2. Humanitarianism and the "brotherhood of man" 

3. Internationalism 

M-. Pride in the achievements of "the human community" 

5. Reverence for life 

6. Human dignity and the "worth of the individual" 


ENVIRONMENT-ORIENTED VALUES 


1. Aesthetic values (environmental beauty) 


2 . Novelty 



This Appendix contains the results of Zelinsky's findings. The 
tables are obtained from the article. Careful examination of the re- 
sults indicates substantial regional variations. 
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TABLE 0 


CoJuiEEJtnoH cocmciiHTit xagtoh loadings and jcomcs ron 
AT ATX S EXTRACTED XftOSt ROTATED MATntX OX 103 SPECIAL-INTEREST MAGAZINES 
and voluntary associations correlated wmi selected social* economic, 

AND DEMOOUArjIIC VARIABLES 
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Fig. 0. R-3: Migrant Factor (factor scores In quintiles). 


Fig. 5. R-2: Innovative West vs. Traditional South Factor (factor scores ir quintiles). 
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DERIVATIONS 


CLASS I: ASSERTION 

I-a. Assertions of facts, experimental or imaginary 

I-b. Assertions of sentiments 

I-c. Mixtures of fact and sentiment 


CLASS II: AUTHORITY 

Il-a. Of one individual or a number of individuals 

Il-b. Of tradition, usages, and customs 

II-c. Of divine beings, or personifications 


CLASS III: ACCORDS WITH SENTIMENTS OR PRINCIPLES 

Ill-a. Accord with sentiments 

Ill-b. Accord with individual interest 

III-c. Accord with collective interest 

Ill-d. Accord with juridical entities 

Ill-e. Accord with metaphysical entities 

Ill-f . Accord with supernatural entities 


CLASS IV: VERBAL PROOFS 

IV-a. Indefinite terms designating real things; indefinite things 
corresponding to terms 

IV-b. Terms designating things and arousing incidental sentiments, 
or incidental sentiments determining choice of terms. 

IV-c. Terms with numbers of meanings, and different things desig- 
nated by single terms 
Metaphors, allegories, analogies 

Vague, indefinite terms corresponding to nothing concrete 


IV-d. 

rv-e. 
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